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FOEEWORD 


The Martin Marietta Thermal Radiation Analyzer System ClRASTS) 
program marks the first instance that thermal radiation analysis 
has been put on the same basis as thermal analysis using program 
systems such as MITAS and SINDA. As with these thermal analyzer 
programs, the user is provided the powerful options of writing 
his own executive or dri'wer .logic and choosing, among several 
available options, the most desirable solution techniques for the 
problem at hand. In addition, TEASYS provides many features never 
before available in a single radiation analysis program. Among 
the more important are: 

- A 1000-node problem size capability with shadowing by inter- 
vening opaque or semitransparent surfaces; 

- A choice of diffuse, specular, or diffuse/ specular radiant 
interchange solutions; 

~ A capability for time-variant geometry in orbit; 

- A choice of analytically determined or externally supplied 
shadow data for environmental flux calculations; 

- Form factors and environmental fluxes computed using an 
internally-optimized number of surface grid elements, selected 
on the basis of user-supplied accuracy criteria; 

- A general editing capability for updating thermal radiation 
model data stored on tape; 

- A plot package that provides a pictorial representation of 
the user’s geometry. 

TRASYS is indebted to a number of predecessor programs in the 
thermal radiation analysis field. The major contributors v?ere 
HEATRATE, MTRAP version 2.0, RAD FAC, and the MRI computer pro- 
gram for detenrining external radiation absorbed by the Apollo 
spacecraft. 

This progratomers ' manual represents an effort to provide scientific 
programming personnel with Che descriptive material necessary to 
reach an understanding of the various program segments. Due to 
the highly modularized design of TRASYS, there are 166 preprocessor 
subroutines and over 300 processor library subroutines described 
herein. Although this is a rather large number in total, the 
individual subroutines have a more moderate size, so that the user 
can develop a working understanding of any routine after devoting 
a reasonable amount of effort to reading the Fortran code and 
consulting Che appropriate material herein. 
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iH This manual was generated for the National Aeronautics and Space 

J Administration's Lyndon B. Johnson Spacecraft Center under NASA 

Contracts NAS9-13033 and NAS9-14318, Development of a Thermal 
Radiation Analysis/Heat Rate Computer Program System. The technical 
monitoring for this program was provided hy Mr. Robert A. Vogt of 
the Thermal Technology Branch of the Structures and Mechanics Divi- 
sion, NAS Lyndon B. Johnson Spacecraft Center. His helpful sugges- 
tions during the development of TEASYS are gratefully acknowledged. 

TRASYS would not exist without the superb design and programming 
efforts of Messrs. Ronald E. Paulson and Robert J. Connor, who 
were responsible for generating the majority of the TRASYS code. 
Their efforts are gratefully acknowledged. Extensive thanks are 
also due to Mr, G. M. Holmstead for his efforts in developing the 
direct irradiation program segment and for the valuable consulting 
effort he performed during the course of program development. Mr. 
Richard G. Goble is also recognized for his praiseworthy efforts 
in developing the specular-diffuse radiation interchange segment, 
the orbit plotter segment, and for his solution.'^ of many knotty 
problems that cropped up during program checkout. 
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INTRODUCTION 


^ffiAT IS TRASYS? 

TRASYS, the Martin Marietta Theriiial Radiaticin Analysis System, is 
a digital computer software system with a generalized capability 
to solve the radiatipn-related aspects of thermal analysis prob- 
lems . VJheh used in conjunction v^ith a generalized thermal analy- 
sis program Siidh aS the Systems Improved 'Numericdi Differencing 
Analyzer (SINDA)' program, any thermal probleni. that can be e:cpressed 
in terms of a lumped -parameter, radiation-conductor thermal net- 
work can be solved. 

The function of TRASYS is twofold. It provides: 

1) Intemode radiation interchange data; 

2 ) Incident and absorbed heat rate data from environmental 
radiant heat sources.' 

Data of both types are provided in a format directly usable by 
the thermal analyzer programs, 

A primary :feature Of TRASYS is that it allows the user to write 
his ohn executive or driver program, which organizes and directs 
the program libraiy routines to solve each specific problem in the 
most expeditious manner. The user also may write his own output 
routines;, thus , the _ data output can directly interface with any 
thermal analyzer using the R-C network concept . 

Other outstanding features of TRASYS include: 

1) A 1000-node allowable problem size; 

2) The ability to accommodate time-variable problem geometry; 

3) An editing capability that allows the combination or separa- 
tion of multiple thermal radiation models; 

4) A plot package that provides pictorial plots of input 
geometry and orbit data, as well as output data. 

The TRASYS system .consists of two majo components: the pre- 

processor, and the processor library i The preprocessor has two 
THajor functions.- First; it reads and converts the user’s geometry 
input data into the fprin used by the processor library routines. 
Second, it accepts the user’s driving logic (inritten in the TRASYS- 
modified FORTRAN iariguage) that directs tiser-prbvided and/or 



library routines in the solution of the problem. The processor 

library consists of FORTRAN language routines that perform the , 

functions cominonly needed by the user. The user has, in so^e J 

cases, a choice of solution techniques for performing the same 
function . 


SYSTEM STRUCTURE 


In the usual engineering environment, a programmer is commissioned 
to prepare an applications program that is subsequently made avail- 
able to the engineer on a production basis , The engineer supplies 
input data and receives output data, as shotm in Figure I-l, 
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Basic Flow in Using an Appli 'ations Program 


In roost cases, changes to the logic and equations are difficult 
for the program user to implement conveniently since they must 
be written in a computer-oriented language and may be submitted 
through a formal programming organization. When TRASYS is used, 
however, the engineer need only call on the programmer to supply 
a standard deck of computer-oriented "control cards" that will 
call the various elements of the system into action in the proper 
sequence. The engineer then fo3cmulates his problem in the 
engineering-oriented TRASYS language, assembling both data and 
solution techniques (i.e., logic and equations) into this card 
deck., which then serves as the complete input to the TRASYS sys- 
tem, Programmer support has been minimized since the bulk of 
the programming effort is already built into the TRASYS pre- 
processor and processor library. The engineering user need only 
specify the data and the order and type of "program building 
blocks" he deems necessary to solve his problem (see Fig. 1-2) . 




Figure 1-2 Basic Flow in Using TRASYS 

It should be evident that TRASYS is much more than an applications 
program. It has, in fact, all the functions and capabilities of 
a special-purpose operating system* Since most computers currently 
used in engineering environments already have operating systems 
built arotand a FORTRAN compiler, TRASYS is designed to augment the 
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existing FORTRAN system. Hence, the TRASYS library serves as an 
extension to the existing FORTR^ library, and the TRASYS program 
serves as a preprocessor to (i.e., it preceeds) the existing 
FORTRAN compiler. This augmentation arrangement is illustrated 
in Figure 1-3. 



Figure 1-3 Internal Flow of TRASYS 

When using the f till capability of TRi^YS , the user will be re- 
o[uired to exert a programming effort of sorts, in a language 
consisting of FORTRAN statements and problem-oriented TRASYS 
statements that are FORTRAN-related. This, together with the 
wide variety of options and features offered by the system, sug- 
gests an app’-opriate word of caution; TRASYS is a comprehensive 
system that cannot be mastered overnight. Nevertheless, to help 
the novice usjf, we have attempted to default much of the required 
input to normally used values so that the user need not define 
them. 
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PART 1 “PREPROCESSOR 


SEGMENT DEFINITIONS 
SEGMENT NAME: TRASYS 


PURPOSE : Driver segment of the TRASYS preprocessor (see Fig. Il-l) . 


CAI.LIN6 SEGMENT : 

None 



SEGMENTS CALLED: 

START 

LOGIGG ■ 

RAPUP 


DATARD 

TPGEN 



Segment NAME-: .start 

PURPOSE: This segment is. the main driver segment for the pfe- 

processor ini talisation, model coilectors source editor, 
and edit output tape genei-ator segments. 


CALLING SEGMENT; 

TRASYS 


SEGMENTS CALLED; 

INITAL 

SEDIT 


SEGMENT NAME : INIT AL 

PURPOSE: This segment initializes the preprocessor-labeled common, 

wtites the TRASiYS banner oh the output file and reads 
in and processes the OPTION. DATA block, 

CALLING SEGMENT; START ' 

SEi^NTS CALLM None 


REV 1 




TJI5ASYS 

CO.O) 


I DATAED 


LOGICCi 


DCMEJRD 

( 2 , 7 ) 


INITAL 

( 1 , 1 ) 


LOGICl 


QUAWRD 

(2,1) 


AWRD 


I 0 GIC 2 


L0GIC3. 


gHDWRD 


S.RFCED 


QRSERD 


Figicra II- 2 Diagram of Dreproce^soi^ySigfrmvti^'^:^^^ 


FLUXRD 

(? 

,1 

1 

























SEGMT . WAtlE : SEDIT 


PURPOSE: Thisr segment’ .performs ’the sotirce edf t functioii of the 

preprocess, or. Input data are read from the INPUT, 
CMERG, EtlERG,’ and RSI files and' output is written 
to, the DA.TAI- and. ESO files , - 

calling .SEGMENT^ S^RT 

SEG^TS, . . Norie- 


SEGMENT NAME: DATARD ■ 

PURPOSE: This segment is the Diain driver segment that reads in 

and processes the user’ s data input block. 

CAIilNG SEGMENT; ' 'TRASYS • ; ' 

SEGMENTS GAtl-ED;; . ^ SDPSSE TO# . 

' Ecs;rd ^ :eeU:E:P : " 

SKIED .DGMNRD , 'CRSPR^^^ 

■ SRFCRD. , ..lERD ' , . .. 


SEGldENT NAtlE : O^^ &NED 

PURPOSE : This segment reads iil and processes the user's QUANTI 

TIES DATA it.put block. 


CALLING SEGMENT: DATARD 

SEGMENTS CALLED: None 


SEGMENT MAMEr - M 

PURPOSE; This segiiieat reads in aad pirpcesses the user’s ARRAY . 

DATA input block. 

CAU.IHG ^GmENT;: dm 

SE6MENTS CALLED^ None . , 

SEGI'IEWT W^4E: SKIRD 

PURPOSE: This segment reads in and processes the user- input "I" 

and *'K” cards of the user’s SURFACE DATA input block. 

CALLING SEGMENT: DATARD 

SEGtiENTS CALLED; None 

SEGMENT NAME ; ' ' SRFCBD 

PURPOSE: This segment reads in the "S”, ”R", "D", ”N’’, and "B" cards 

of the user’s SURFACE DATA input block. These data are 
combined with the output of segment SKIRD and are output 
on file NRIOS for final processing in segment SDPSS2. 

CALLING SEGl-IENT; DATARD 

SEGMENTS CALLED;. None 

SEGMENT NAME; SDPSS2 

PURPOSE ; This segment comp 1 e tes the processing of the user’s 
SURFACE DATA input block. 

CALLING SEG1.1ENT: DATARD 

SEGMENTS CALLED; None 
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SEGMENT NAME; BCSRD 


PURPOSE: This segment reads in and processes the user’s BCS DATA 

input block. This segment also reads the NODE/BCS 
directory generated by segment SDPSS2 and processes 
it with the user's input data to form the communica- 
tion link between the BCS data and the SURFACE data. 

CALLING SEGMENT: DATARD 

SEGMENTS CALLED: None 


SEGMENT NAME : DCMNR D 

PURPOSE: This segment reads in and writes out to the system 

output file the user's DOCUMENTATION DATA input block. 

CALLING SEGMENT: DATARD 

SEGMENTS CALLED: None 

SEGMENT NAME: FFRD 

PURPOSE: This segment reads in and processes the user's FORM 

FACTOR DATA input block, 

CALLING SEGMENT: DATARD 

SEGMENTS CALLED; None 

SEGMENT NAME : SHOWE D 

PURPOSE; This segment reads in and processes the user's SHADOW 
FACTOR DATA input block. 

CALLING SEGMENT: DATARD 

SEGMENTS CALLED: None 
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SEGMENT NAME : FLUKED 


PURPOSE; This segment reads in and processes the user’s FLUX 
DATA input block. 

CALLING SEGMENT: DATAED 

SEGMENTS CALLED: None 



8 


SEGMENT NAME: CRSPRD 

PUFEOSE; This segment reads in and processes the user's CORRE- 
SPONDENCE DATA input block. 

CALLING SEGMENT: DATARD 

SEGMENTS CALLED: Nona 


SEGMENT NAME: LOGICO 

PURPOSE: This segment is the driver segment for the segments 

that read in and process the user's OPERATIONS DATA 



block. 


CALLING 

SEGMENT : 

TRASYS 

SEGMENTS 

CALLED: 

LOGICl 

L0GIC2 

LOGICS 


SEGMENT NAME; LOGICl 

PURPOSE: This routine reads in the user's OPERATIONS DATA block 

and sets up all the variables needed to write the 
ODPROG segment of the processor. 

CALLING SEGMENT: LOGICO 

SEGMENTS CALLED: None 
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SEGMENT NAME; L0GIC2 


PURPOSE; This routine reads the output of segment LOGiGl and 
writes to the processor compile file CMPL, the main 
processor segment, TRASYS, and the ODPROG subsegment. 

CALLING SEGMENT: LOGICO 

SEGMENTS CALLED: None 


SEGMENT NAME; LOGICS 

PURPOSE: This segment reads in and processes the user’s SUBROUTINE 

DATA block and writes to the processor compiler file 
CMPL all subsegments of the processor that require com- 
pilat ion. 

CALLING SEGMENT; LOGICO 
SEGMENTS CALLED; None 


SEGMENT NAME; TP GEM 

PURPOSE; This routine writes the needed driver information of 
the processor to the sequential data file SQNTL. 

CALLING SEGMENT: TRASYS 

SEGMENTS CALLED: None 


SEGMENT NAME ; RAP UP 

PURPOSE; This segment wraps up the preprocessor execution phase. 
CALLING SEGMENT: TRASYS 

SEGMENTS CALLED: None 
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B. SUBROUTINE AND FUNCTION DESCRIPTIONS - PREPROCESSOR 

ROUTINE NAME: AAAAAA 

DESCRIPTION ; This routine initializes the variables containing the 
la,st program modifications number and date. 

CALLING SEQUENCE: CALL AAAAAA (V, D) 

V - Last version modification number 
D - Date of last modification 

REFERENCED BY: SEGMENT ROUTINE 

INITAL INITAL 


ROUTINE NAME : ABNORMl (CDC system routine) 

DESCRIPTION: This routine when called causes the program to termi- 

nate abnormally with error traceback. 

CALLING SEQUENCE: CALL ABNORMl (PI, P2, P3) 

PI - The name of the calling subroutine; left- 
justified, Hollerith input 
P2 - A decimal number, maximum of 88, which is 
used as an error number. Must not be 0 
P3 - The error message; lef t-justified zero- 
filled. The message must be terminated 
with 4 octal zeros in the rightmost 
position of a word 

REFERENCED BY; SEGMENT ROUTINE 

TRASYS ABTl 

INITAL OPTNRD 
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ROUTINE NAME : ABTl 


DESCRIPTION: This is an abnormal exit routine for the TRASYS pre- 

processor, 

CALLING SEQUENCE: CALL ABTl (N, NO) 

N - Name of the routine exiting from. The 

name is left-justified and the rest of the 
field blanked within the word 
NO “ Type of error 

=1 User input error 
= 2 Bad source edit input tape 
= 3 TAPE/DISK/DRUM read error 
= 4 Program limitations exceeded 
= 5 Job field length too short 
= 6 Programmer error 


REFERENCED BY: SEGMENT ROUTINE SEGMENT ROUTINE SEGMENT ROUTINE 


TRASYS 

HCARD 
INRDD 
TPERR 
SORTD 
S ERROR 

SRFCRD 

SEDIT 

SEDIT 

SDPSS2 


INRDPP 

INRDE 

BCSRD 

QUANRD 

QUANRD 

QUANSD 

FFRD 

ARRYRD 

ARRYRD 

SHDWRD 

SKIRD 

IPl 



IP2 



KP2 



SRFCSl 

PLUXRD 

PLUXRD 

CPDRC3 

REFCD 

CRSPRD 

CRSPRD 

DUPSRF 

LOG ICO 

WCMMN 


LOGICl 

LOGICl 

SDPSS2 


LPl 

BCSRD 


ORBGEN 

BCSPl 


ORB IN 

BCSP2 


ORBNPL 

FFPROC 


ORBPL 

SHDmD 


ORBPLC 

BLDGEN 

SFRSI 

RAPUP 

RAPUP 
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ROUTINE NAME; AERYRD 

DESCRIPTION: This routine reads in and processes the user array 

input block. 

CALLING SEQUENCE; GALL ARRYRD 
REFERENCED BY: SEGMENT ROUTINE 

DATARD DATARD 

FILES: NRI02" Writes three records to the random l/O file 

Record 1 ~ Array names /input card sequence number 
Record 2 - Array starting positions with respect to 
the output array common/length of the 
array +1 

Record 3 - Array data 

NSCl “ Temporary file to store the array data 
NSG2 " Temporary file to store the array lengths 
NODT - System output file 


ROUTINE NAME: ASKCRD 

DESCRIPTION: This routine decodes the user input source edit con- 

trol cards for the preprocessor source editor. 

CALLING SEQUENCE; CALL ASKCRD (lER) 

lER “ Error detection flag 
= 0 No error found 
= 1 One or more errors found 

KEY VARIABLES: Card data passed to this routine is common/ CARD/ 

variables NCDTYP and IND. Decoded information is 
retu3med in cononon/CARD/ variables IND^ NDOT, IDOT, 
JDOT. 

REFERENCED BY: SEGMENT ROUTINE 

SEDIT SEDIT 

FILES; NOUT - System output file 


t 
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ROUTIHE NME: ATJTOCM 

DESCRIPTION: This routine sets up a correspondence data file on 

NGBIRR for the automatic combination of surfaces 
generated by a user-input polygon. 

CALLING SEQUENCE! CALL AUTOCM 

REFERENCED BY; SEGMENT ROUTINE 

SDPSS2 POLYGN 

FILES: NGBIRR 

ROUTINE NAME; BANHPP 

DESCRIPTION: This routine prints the banner page for the TRASYS 

preprocessor on the output file. 

CALLING SEQUENCE: CALL BANNPP 

REFERENCED BY: SEGMENT ROUTINE 

INITAL INITAL 

FILES: NOUT - System output file 
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ROUTIHE NAME': gCSPl 

DEScMptION; This ibtitiirie pio^ the Sirs t pass on the user- 

. iaipht BGS data bleckh 

CALLING SEQUEUGE: CALL BC?P1 (NIX) 

1^^ used- in dynamic storage, blank 

common variable (IX) 

O 

KEY VARIABLE; IX - Data once converted . from user input are 

■ passed to the BCSP2 routine via this blank 
common variable 

REFERENCED BY; SEGMENT ROUTINE 

BCSRD BCSRD 

FILES; NOUT - System output file 


ROUTINE NAtiE; BCSP2 

DESCRIPTION; This routine processes the block coordinate system 
input into data values and a BCS directory. Data 
are input to this routine in blank common generated 
by routine BGSPl. Data leave the routine via 
blank common. 


GALLING SEQUENCE; CALL BCSPl (NIX, NBSD, NBND, NBSV* NBEV) 

NIX - Length of vector data in blank conpion 
input to this routine 

NBSV - Starting word of BCS data values in 
blank common (always 1) 

NBEV - Length of BCS data values in blank 
common 

NBSD - Starting word of BCS directory in 
b lank common 

■ NBWD - Length of BCS directory in blank common 


REFERENCED BY: SEGMENT 


ROUTINE 


BCSRD BCSRD 


FILES: NGUT - System, output file 
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ROUTINE NME: BGSRD 

DESCRIPTION: This roudine and the routines that it calls read 

in and process the user-input BCS DATA block and 
tie the user- input surface data to the BCS input 
data. Besides the user BCS data inputs this rou- 
tine reads a random I/O record that was output by 
the surface data block's pass 2 processor. 

CALLING SEQUENCE: CALL BCSRD 

KEY VARIABLES: NRXS — Random l/O record number of record written 

by surface data pass 2 processor 
NLIS - Length of the random I/O record 

REEERENCED BY: SEGMENT ROUTINE 

DATARD DATARD 

FILES: NOUT - System output file 

NRIO - Random I/O file 


ROUTINE NAME; BLDGEN 

DESCRIPTION: This routine decodes the "BUILD” cards found in the 

operations cata block and writes the corresponding 
calls to processor subroutines BUILDC and ADD, 

CALLING SEQUENCE: CALL BLDGEN 

REFERENCED BY: SEGMENT ROUTINE 

LOGICl LPl 

FILES; NOUT - System output file 

NSCl - Sequential scratch file 
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ROUTINE NAME: BOX 

DESCRIPTION: This routine sets up the Euler angles, position vector, and 

program-compatible surface description parameters for the top 
of a five- or six-sided box input by the point method. 

calling SEQUENCE: GALL BOX (Pll , P12, P13 , P31, P32, P33) 

Pll, E12, P13 - Xj Y, and Z components of a vector 
CO linear with the surface coordinate system Y-axis 

P31, P32, P33 - Xj Y, and Z components of a vector col- 
inear with the surface coordinate system X-axis 

REFERENCED BY: SEGMENT ROUTINE 

SDPSS2 SDTPS2 

FILES; NOUT - System output file 



ROUTINE NAME: BOXGEN 

DESCRIPTION: This routine generates Euler angles, position vector, and 

surface description parameters for the sides and bottom 
(B0X6) of a box based on similar parameters set up by BOX 
for the top of the box. 

CALLING SEQUENCE: CALL BOXGEN (FU, P12, P13 , P31, P32, P33, PHI, PSI, 

OMG, ISURF) 

Pll, P12, P13 - X, Y, and Z components of a vector colinear 
with the surface coordinate system Y-axis 

P31, P32, P33 - X, Y, and Z components of a vector colinear 
with the surface coordinate system X-axis 

PHI, PSI, OMG Euler angles necessary to rotate the ICS, 
BCS, or CCS into the surface coordinate system 

ISURF “ Counter indicating which side of box is being 
generated (1 ISURF ^ 6) 

REFERENCED BY; SEGMENT ROUTINE 

SDPSS2 SDTPS2 
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ROUTINE NAME: BSRGHD 

DESGRIPTION: This routine perfpxnns a table lookup for exact equal com- 

pares. The table must be a doublet array and the routine 
uses the binary search technique. 


GALhING SEQUENCE: GALL BSRGHD (NAME, ITABLE, NROINT, MIDPT, IAN) 

NAME - Independent variable name searched for 
ITABLE - Doublet table (NAME./ ANSWER) 

NPOINT - Number of words in the table (2 * number of 
entries) 

MIDPT - Binary midpoint of table (calculated in this 
routine of MIDPT - 0) 

IAN - Dependent variable found (ANSWER) 

= 0 if no match is found 


REFERENCED BY: SEGMENT 


ROUTINE 


INITAL OPTNRD 

SRFCRD SRECSl 

CDPRG3 

FFRT FFREAD 

SHDWRD SFTBRD 

SFRSI 

FLUXED DIBLDR 


ROUTINE NAME: CALAC 

DESCRIPTION: This is a function subroutine that calculates surface and 

node areas. 


CALLING SEQUENCE: CALAC (ILK, ALPR, BMIN, BMAX, GMIN, GMAX) 


ILK - Surface type 

1 = Rectangle 

2 = Disc 

3 = Trapezoid 

4 = Cylinder 

5 = Cone 

6 = Sphere 

7 = Paraboloid 

ALPH - \ 

BMIN - I 

BMAX - ) description parameters 

GMIN - V 
GMAX - I 


REFERENCED BY: SEGMENT 


ROUTINE 


SDPSS2 SDTPS2 
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ROUTiWE NAME; gALAI^G 


DESCRIPTION: Given the diresction cbsines . of SGS X and Y-;a;;icis^ 

IGSi BGS-, -on CCS> this; ,E(iu|iiie^ c : a; set Of' .Eulen . , 

angles to rotate the ICS, BCS, or CCS into the SCS . 

CALLING SEQUENCE : CALL GALANG (Pll , PI2 E13, P31, :P32 , P33;,^ ROfZ^ . 


Pll, P12, P13 - X,. y , and Z . eompohents of a yectbr with 
magnitude BMAX that is col inear V7ith the SCS Y-axis 

P31, P32, P33 - X, Y, and Z components of a vector with 
magnitude GMA3C that is co linear with the SCS X-axis 


ROTZ, ROTY, ROTX - Euler angles to rotate ICS, BGS, or 
CCS into the SCS 


REFERENCED BY; SEGMENT 
SDPSS2 


ROUTINES 


BOX 

PARAB 

BO'XGEN 

TRAPZ 

CONE 

RECT 

CYLNDR 

SPHERE 

DISC 


IMAGES 



FILES : NOUT - System output file 


ROUTINE NAME: CALB 

DESCRIPTION: This routine calculates the BMIN and BMAX values for each 

node. 


CALLING SEQUENCE; 


CALL CALB (DB, BETA, IB, BMAXT, NVB, IG, DVMUL) 


DB - Measure of node width in the beta direction 

BETA - Measure from the edge of the surface to the center 

of the current node in the beta direction 
IB “ Sequence number of current node in the beta direction 
BMAXT - Temporary storage of BMAX value for surface 
NVB - Number of nodes in the beta direction 

IG - Sequence number of current row of nodes in the gamma 

direction 

DVMtJL - Length unit multiplier from a D-card . 


REFERENCED BY : SEGMENT 

SDPSS2 


ROUTINE 

SDTPSZ 


FILES: NOUT - System output file 
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ROUTINE NAME: CALG 

DESGRIPTION: This . routine calculates- the, GKtN '.and- SliiS va^ fcr: 

each node. , 

CAEIING SEqUENGE.; CALL CALG (DG/GANMA,'IG, GlliXT, NVGyDVMUt) 

. . _ ^ - Measure of node length in ganuna direction 

. GAMH-JA - Measure from the beginning edge Of the 
- : surface :tp th& cehte^^ current node 

, . in the gamma direction 

IG - Sequence number of the current node in 

GMAXT - Temp or ary : s torage Of GMAX value for . 
surfaCe-- 

WG j: -Numher of nodes in gairana direction 

D\fMUL - Length unit : mu ^ a D-card 

REFERENCED BY; SEGMENT ROUTINE^ 

SDPSS2 SDTPSY 
FILES : NOUT - System Output file 
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ROUTINE FAME: CDRRCl 


DESCRIPTION: This rbatine perfoxnis the first-pass conversion of a card 

data input that was read using an "A” format to Hollerith, 
integer, and floating-point words 


CALLING SEQUENCE: GALL GDPRGl 


KEY VARIABLES; 


REFERENCED BY: 


IND - Array containing the card image data in 
Hollerith format 

NDOT - Array containing the conver-ted values 
I DOT - Array describing what is contained in the 
NDOT array. 

JDOT - Column location in which the values in the 
NDOT array started on with respect to the 
IND arrays 

ID - A single word containing the number of words 


used 

in the NDOT 

array 


SEGMENT 

ROUTINE 

SEGMENT 

ROUTINE 

trasys 

HGARD 

SlffiWRD 

SHOWED 

QUANRD 

QUANRD 


SFTBRD 

ARRYRD 

AKSYRD 

ELUXRD 

FLUIOID 

SKIRD 

SKIRD 


DIBLDR 


KPl 

CRSPRD 

CRSPRD 


IPl 

LOGICl 

LPl 

BLDGEN 


SRFCRD 

SFRCRD 

ORBGEN 


SRFCSl 

LOGIC 2 

L0GIC2 


SRFCBC 

LOGICS 

LOGICS 

BCSRD 

BCSPl 



FFRD 

FFREAD 

FFRD 




FILES; NOUT - System output file 
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ROUTINE NAME : GDPRG3 


DESCRIPTION: This routine compares all Hollerith words in the NDOT array 

to a directory array residing in the dynamic storage area 
for a matching name. If it finds one, then the name in the 
NDOT array is replaced by the corresponding data value or 
values. This routine should be called after CDPRCl and 
before CDPRC2. 


CALLING SEQUENCE: CALL CDPRG3 


KEY VARIABLES: NDOT - Output array from CDPRCl and input for this routine 

Also the output from this routine 
IDOT Ss JDOT - Same as NDOT 

IX - Dynamic storage (blank common), containing the con- 
stant directory and constant data 
IXSD ~ Index pointing to first word of directory 
IXWD - Length of directory 

IXSV - Index pointing to first word of data values 


REFEREI'IGED BY: SEGMENT 


ROUTINE 


SKFCRD SRFCSl 


FILES: NOUT - System output file 


ROUTINE NAME : CHEC 

DESCRIPTION: This routine checks the validity of surface description 

parameter values. 

CALLING SEQUENCE: CALL CHEC 

REFERENCED BY: SEGMENT ROUTINE 

SDPSS2 SDTPS2 

FILES: NOUT - System output file 

ROUTINE NAME : ^EAR 

ENTRY POINT : CHECTP 
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ROUTINE NAME; CDPRC2 


DESCRIPTION: This routine processes the arithmetic calculations (*, /, +) 

within a data field input on a user input card. This routine 
is usually called after subroutine CDPRCl, but in some instances 
it is cp.lled after CDPRC3. 


GALLING SEQUENCE: GAI.L GDPRC2 


KEY VARIABLES; KDOT, IDOT, JDOT and ID have the same function in this 

routine as they do in CDPRCl. This routine may or may not 
condense the arrays, depending on the type of input. 


REFERENCED BY: SEGMENT ROUTINE 


TRASYS HOARD 

QUANRD QUANRD 

ARRYED ARRYRD 

SKIRD SKIRD 

KPl 
IPl 

SRFCRD SRFCRD 

SRFCSl 
SRFGBC 

BCSRD BCSPl 

FFRD FFRD 

PFREAD 

SHDlfflD SHDWRD 

SFTBRD 

FLUXRD FLUXRD 

DIBLDR 
CRSPRD CRSPRD 

LOCI Cl ORBGEN 

BLDGEN 

PILES; NOUT “ System output file 


ROUTINE NAME : CLEAR 

DESCRIPTION: This routine frees the buffers and rewinds the 

unit so that it may be re-used, (Generally 
used on scratch files.) 

CALLING SEQUENCE: CALL CLEAR (lUNIT) 

lUNIT - Unit number to be freed. 

REFERENCED BY: SEGMENT ROUTINE 

FFRD PFRPROC 

SRFCRD SRFCRD 
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ROUTINE NAME 


WFFCOH 


DESCRIPTION: This routine name is an entry point into the routine WCOM, 

When this entry point name is called j the Fortran-labeled 
common for the processor PFPROG segment is written to the 
NCMPL file. 

CALLING SEQUENCE: CALL WFFCOM 

REFERENCED BY; SEGMENT ROUTINE 

LOGICO WPROG 

LOGICS LOGICS 

FILES: NCMPL - The file that the generated processor Fortran cards are 

written to 


ROUTINE NAME; WSFCOM 

DESCRIPTION; This routine name is an entry point into the routine WCOM. 

Wlien this entry point name is called, the Fortran-labeled 
common for the processor SFPROG segment is written to the 
NCMPL file, 

CALLING SEQUENCE; CALL WSFCOM 

REFERENCED BY: SEGMENT ROUTINE 

LOGICO WPROG 

LOGICS LOGICS 

FILES; NCMPL - The file that the generated processor Fortran cards are 
written to 
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ROUTINE NAME; 


WNPGQM 


DESCRIPTION: This routine naiiis is sin entry point into the routine WCOM. 

When this entry point is called,, the Portran-labeled common 
for the processor NPPROG segment is written to the NCMPL 
file. 

CALLING SEQUENCE: CALL WNPCOM 

REFERENCED BY; SEGMENT ROUTINE 

LOGICO WPROG 

LOGICS LOGICS 

FILES: NCMPL - The file that the generated processor Fortran cards are 

written to 


ROUTINE NAME; WOPCOM 

DESCRIPTION; This routine name is an entry point into the routine WCOM. 

When this entry point name is called, the Fortran- labeled 
common for the processor OPPROG segment is written to the 
NCMPL file. 

CALLING SEQUENCE: CALL WOPCOM’ 

REFERENCED Bv- SEGMENT ROUTINE 

LOGICO WPROG 

LOGICS LOGICS 

FILES: NCMPL - The file that the generated processor Fortran cards are 

written to 


ROUTINE NAME: WDICOM 

DESCRIPTION: This routine name is an entry point into the routine WCOM. 

\-Ihen the entry point is called, the Fortran-labeled common 
for the processor DIPROG segment is written to the NCMPL 
file , 

CALLING SEQUENCE: CALL WDICOM 

REFERENCED BY; SEGMENT ROUTINE 

LOGICO WPROG 

LOGICS LOGICS 


FILES : NCMPL 


The file that the generated processor Fortran cards are 
written to REV 1 

11-23 


I 


ROUTINE NAME; WGBGOM 


DESCRIPTION: This routine name is an entry point into the routine 

WGOM, When this entry point name is called, the Fortran- 
labeled common for the processor* GBPROG segment is written 
to the NGMPL file. 

I CALLING 3EQTXENGE: CALL WGBGOM 

REFERENCED BY: SEGMENT ROUTINE 

LOGICO WPROG 

LOGICS LOGICS 

FILES: NGMPL - The file that the generated processor Fortran cards are 

written to 


I ROUTINE NAME: WGMCOM 

DESCRIPTION: This routine name is an entry point into the routine WCOM. 

When this entry point name is called, the Fortran-labeled 
I common for the processor CMPROG segment is written to the 

NGMPL file. 

I GALLING SEQUENCE: CALL WCMCOM 

REFERENCED BY: SEGMENT ROUTINE 

LOGICO WPROG 

LOGICS LOGICS 

FILES; NCMPL - The file that the generated processor Fortran cards are 
written to 


I ROUTINE NAME: WAQGOM 

DESCRIPTION: This routine name is an entry point into the routine WCOM. 

When this entry point name is called, the Fortran- labeled 
common for the processor AQPR0G segment is written to the 
NCMPL file. 

I CALLING SEQUENCE: CALL WAQCOM 

REFERENCED BY: SEGMENT ROUTINE 

LOGICO WPROG 

LOGICS LOGICS 

FILES: NCMPL - The file that the generated processor Fortran cards are 

written to 
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ROUTINE NAME: 


WQOGOM 


DESCRIPTION: This routine name is an entry point into the routine WGOM. 

When this entry point name Is called, the Port ran- labeled 
common for the processor QOPROG segment is written to the 
NGMPL file. 

CALLING SEQUENCE: CALL WQOCOM 

REFERENCED BY: SEGMENT ROUTINE 

LOGICO ^WROG 

L0GIC3 LOGICS 

PILES: NGMPL - The file that the generated processor Fortran cards are 

written to 


ROUTINE NAME: WRBCOM 

DESCRIPTION: This routine name is an entry point into the routin e WCO M 

When the entry point name is called, the Fortran-labeled 
common for the processor RBPROG segment is written to 
NGMPL file. 

CALLING SEQUENCE: CALL WRBCOM 

REFERENCED BY: SEGMENT ROUTINE 

LOGICO WPR06 

LOGICS LOGICS 

FILES: NGMPL - The file that the generated processor Fortran cards are 

written to 


ROUTINE NAME: WRCGOM 

DESCRIPTION: The routine name is an entry point into the routine WCOM. 

When this entry point name is called, the Fortran-labeled 
common for the processor RCPROG segment is written to the 
NCMPL file. 

CALLING SEQUENCE: CALL WRCCOM 

REFERENCED ,BY: SEGMENT ROUTINE 

LOGICO WPROG 

LOGICS LOGICS 

FILES: NCMPL - The file that the generated processor Fortran cards are 

written to 
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ROUTINE NAME: WDRCOM 

DESCRIPTION: This routing riante is an entry point into the routine WCOM. 

When this entry point name is called, the Fortran” labeled 
coiMnon for the processor DEPROG segment is written to the 
NCMPL file. 

GALLING SEQUENCE: CALL WDRCOM 

REFERENCED BY: SEGMENT • ROUTINE 

LOGIGO WPROG 

LOGICS LOGICS 

FILES: NCMPL - The file that the generated processor Fortran cards are 

written to 


ROUTINE NAME: COME 

DESCRIPTION: This routine converts point input for cones to program- 

compatible surface description parameters and sets up 
Euler angles and a position vector to transform the ICS, 
BCS, or CCS into the PCS of the cone. 

CALLING SEQUENCE: CALL CONE 


REFERENCED BY: 

SEGMENT 

ROUTINE 


SDPSS2 

SDTPS2 

FILES: NOUT - 

System output 

file 
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ROUTINE NAME: GRSPRD 

DESCRIPTION: This routine and the routines that it calls read in and 

'process the user^.s CORRESPONDENCE DATA bio cfe. 

CALLING SEQUENCE: CALL GRSPRD 

REFERENCE: SEGMENT ROUTINE 

MTARD DATARD 

FILES: NGBIRR - Contains the correspondence data output from this routine 

for the processor phase. 

NOUT - System output file 

ROUTINE NAME; CELHDR 

DESCRIPTION: This routine converts point input for cylinders to pro- 

gram-compatible surface description parameters and sets 
Euler angles and a position vector to transform the ICS, 

BCS, or CCS into the SCS of the cylinder. 

CALLING SEQUENCE; CALL CYLNDR 

REFERENCED BY: SEGMENT ROUTINE 

SDPSS2 SDTPS2 

FILES: NOUT - System output file 


ROUTINE NAME: PATARP 

DESCRIPTION: This routine is the driving segment that calls the other 

segments that read in and process the data blocks of the 
user’s input data deck. It does not read in the logic blocks. 

calling SEQUENCE: GALL DATARD 

REFERENCED BY: SEGMENT ROUTINE 

TRASYS TRASYS 

FILES: NOUT - System output file 
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ROUTINE NAME: DATE (ODC system routine) 

DESCRIPTION: This routine returns the current date that the job was run 

on the computer. The data is returned in display code in 
the following format.... 

bm/DD/YY. 

CALLING SEQUENCE: CALL DATE (DTE) 

DTE - Returned data variable 
REFERENCED BY: SEGMENT ROUTINE 

INITAL INITAL 

ROUTINE NAME: DCMNRD 

DESCRIPTION: This routine processes the user-input documentation data 

block, 

CALLING SEQUENCE: CALL DCMNRD 

REFERENCED BY: SEGMENT ROUTINE 

DATARD DATARD 

FILES; NOUT - System output file 
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ROUTINE NAME: DDIMP 


DESCRIPTION: This routine prints consecutive core memory words to 

the system output file in octal format. This routine is 
used to dump the error trace information of the pre- 
processor. 

CALLING SEQUENCE; CALL DDUMP (ISADD, lEADD, NN) 

ISADD - Address of word to start dumping 
lEADD - Address of word to end dumping 
NN - Integer number to print to identify the 
printed dump 

REFERENCED BY: SEGMENT ROUTINE SEGMENT ROUTINE SEGMENT ROUTINE 


SEDIT 

SEDIT 

SKIRD 

SKIRD 

CRSPED 

CESPRD 


EPTAPE 


IP 2 




OTWTE 


KP2 

LOGICl 

LOGICl 


MERGE 

SRFCRD 

SRFCRD 



qUANRD 

QUANRD 


SRFCSl 


LPl 



BCSRD 

BCSRD 

TPGEN 

TPGEN 


QUANSD 





arryrd 

ARRYRD 

SHDWRD 

SHDWRD 

FFRD 

FFRD 




SFTBRD 


FFPROC 




SFRb ■ 


FFREAD 




SFWPL 

TRASYS 

TRASYS 


SFCIO 

FILES: NOUT - System output file 
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ROUTINE NAME; 


DIBLDR 


DESCRIPTION; This routine process one set of user's flux data input 
in the FLUX DATA input block. The output of this 
routine is the flux data restart file, NDIR. 

CALLING SEQUENCE: CALL DIBLDR (NODIM, IQDS, IQDR, IQDP, NN, NN2, INIT, 

ISTEP) 

NODIN - Starting location for the unsorted node array 
IQDS - Starting location for the solar data 
IQDR “ Starting location for the albedo data 
IQDP “ Starting location for the planetary data 
NN - Number of nodes in the node array 
NN2 - Number of xjords in the node directory (2 * NN) 
INIT - Value to initialize the IQDS, IQDR, IQDP arrays 
to 

ISTEP - Logic block step number to which this set of 
flux data pertains 

REFERENCES : SEGMENT ROUTINE 

FLUXRD FLUXRD 

FILES: NOUT - System output file 

NDIR - Flux data restart file 


ROUTINE NAME: DIRCS 


DESCRIPTION; Given the Euler angles necessary to rotate one coordinate 
system into another, this routine calculates the corre- 
sponding direction cosines. 

CALLING SEQUENCE; CALL DIRCS (II, J J , KK, PHI, PSI, OMI, TRAN) 


II, JJ , KK - Integers (1, 2, or 3) defining the order 
in which the rotations' PHI , PSI, and OMI 
are to be performed 

PHI, PSI, OMI - Euler angles defining the rotations 

about the Z, Y, and X-axes, respectively 
TRAN - A three-by- three matrix of direction cosines 


REFERENCES: SEGi^lENTS ROUTINES 

SKIRD IPl 

SDPSS2 SDTPS2 

BOX 
IMAGES 
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ROUTINE NAME; DISC 

DESCRIPTION: This routine converts point input for discs * to program- 

compatible surface description parameters and sets up 
Euler angles and a position vector tp transfo;rm the IGS, 
BCS» or CCS into the SCS of the disc. 

CALLING SEQUENCE: CALL DISC 

REFERENCED BY; SEGMENT ROUTINE 

SDPSS2 SDTPS2 

FILES; NODT - System output file 


ROUTINE NAME: DUPSRF 

DESCRIPTION; This routine duplicates previously input surfaces and 
modifies surface description parameters as specified 
by the user. 

CALLING SEQUENCE; CALL DUPSRF (ISSX» lESX, lER) 

ISSX - The index of the starting location where 
previously input surface descriptions are 
loaded into blank common 
lESX - The index of the location of the end of 
the surface description data in blank 
common for a previously input surface 
lER - An indicator of errors encountered 

REFERENCED BY: SEGMENT ROUTINE 

SRFCRD SRFCSl 

FILES; NOUT - System output file 

NRIOS - Scratch random access file 
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I DESCRIPTION: This routine positions the NRSI file at the beginning of the 

request so that source editing of the user's data 

^ cah.^Begin. 

CALL EPTAPE (NU, NiBiE, lER) . 

NES'I - ■ Restart 'injpnt'^’^ 

NAME " Name ofmbdel requested for editing 
lER - Error return flag . ' 

= 0 Ready for editing , 

^ 0 Not ready for editing 

REFERENCED BY: SEGMENT - RODTINE 

SEDIT SEDIT 

^ FILES: NODT ~ System output file 

I HRS I - Restart input file 


GALLING SEQUENCE: 



ROUTINE NAME: FDPRC 

DESCRIPTION: This routine does the **, - calculations between 

constants input within a single data field on a user- input 
data card. 

calling SEQUENCE: CALL FDPRC (TEL, IPR, NT-JP) 

Starting location within NDOT array to start 
performihg calculations 

Ending location to stop processing calculations 
Number of words left after calculation processing- 
that contains data within the data field processed 

REFERENCED BY: SEGMENT ROUTINE 

TEASYS GDPRG2 

FILES: NOUT - System output file 


IPL - 

IPR - 
NWP - 
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ROUTI'NE NAME : FFPROC 


DESCRIPTION: This routine reads in and processes cards in the form 

factor data block and writes the information to the 
NFFR file. 

CALLING SEQUENCE: CALL FFPROC (NN2,ISKIP) 

NN2 - Twice the number of nodes in the surface 
data block 

ISKIP - Flag to indicate form factor data proces- 
sing was skipped due to insufficient 
’ field length. 

REFERENCED BY; SEGMENT ROUTINE 

FFRD FFRD 

FILES: NOUT - System output file 

NRIOS - Random access data file 
NPPR - Form f-(-tor data file 
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ROUTINE NAME : FFRD 


DESCRIPTION: This routine decodes all form factor data block input 

cards. 

CALLING SEQUENCES: CALL FFRD 

REFERENCED BY: SEGMENT ROUTINE 

DATARD DATARD 

FLIES: NOUT - System output file 

NFFR “ Form factor data file 


ROUTINE NAME ; FFREAD 

DESCRIPTION: This routine decodes form factor data cards that 

contain node-to-node form factor data. 

CALLING SEQUENCE: CALL FFREAD <NN2 , NODATA, lERCOR) 

NN2 - Twice the numbet of nodes in the surface 
data block 

NODATA - Counter for number of new data entries 
written to NSCl 

lERCOR - Counter for number of card input errors 
found 

REFERENCED BY: SEGMENT ROUTINE 

FFRD FFPROC 

FILES: NOUT - System output file 

NSCl - Sequential scratch file 
NFFR - Form factor data file 
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ROUTIHE NAME; 


FFSORT 


DESCRIPTION; This routine sorts the form factor input data into 
the order expected by the form factor calculation 
segment. 

CALLING SEQUENCE: CALL FFSORT (lA, IB, NN) 

lA “ Unsorted array 
IB - Sorted array 

NN - Number of values in lA and IB 
REFERENCED BY: SEGI-IENT ROUTINE 

FFB,D FFPROC 

FILES : None 

ROUTINE NAME ; FLUXRD 

DESCRIPTION: This routine and the routines that it calls read in 

and process the user’s input FLUX DATA block. 

CALLING SEQUENCE: CALL FLUXRD 

REFERENCED BY: SEGMENT ROUTINE 

DATARD DATARD 

FILES: NOUT - System output file 
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ROUTINE NAME : GNRSO 

DESCRIPTION: This routine generates the first file of the RSO (restart) 

tape. 

CALLING SEQUENCE: CALL GNRSO 

REFERENCED BY: SEGMENT ROUTINE 

SEDIT SEDIT 

FILES: NRSO Restart output tape 


ROUTINE NAME; HOARD 

DESCRIPTION: This routine decodes the header cards of the user's input 

data deck when the card is read from the NDI file. 

CALLING SEQUENCE: CALL HOARD (N) 

N - 1 Card to be processed has been read prior to 
calling this routine 

N = 0 Card is to be read before processing 


REFERENCED BY: 

SEGMENT 

ROUTINE 


START 

HEADCD 


INITAL 

OPTNRD 


SEDIT 

SEDIT 


datard 

DATARD 


QUANRD 

QUANRD 


ARRYRD 

ARRYRD 


SKIRD 

SKIRD 


SRFCRD 

SRFCRD 

SRFCSl 

SRFCBC 


BCSRD 

BCSRD 


DCMNRD 

DCMNRD 


FFRD 

FFRD 

FFREAD 


SHDWRD' 

SHDWRD 

SFTBRD 


FLUXRD 

FLUXRD 

DIBLDR 


CRSPRD 

CRSPRD 


LOGICl 

LPl 

BLDGEN 


LOGICS 

LOGICS 

FILES ; NOUT - 
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System output 

file 



ROUTINE NAME: 


HEADCD 






DESCRIPTION: This routine decodes the HEADER cards o£ the user input 

data deck when the card is to be read directly from the 
system input file. 

.CALLING SEQUENCE: CALL HEADCD (N) 


N = 1 Card to be processed has been read prior to 
calling this routine 

N = 0 Card is to be read before processing 


REFERENCED BY: SEGMENT ROUTINE 


INITAL INIXAL 


FILES: NIN System input file 

NOUT “ System output file 


ROUTINE NAME: HTI 


DESCRIPTION: This routine converts a Hollerith word containing a positive 

integer number to an integer number. 

CALLING SEQUENCE: CALL HTI (IN, lOUT, IFLAG) 

IN - Word to be converted 

lOUT - Converted word 

IFLAG ^ = 0 Conversion was complete 

= 1 Conversion was not complete 

REFERENCED BY: SEGMENT ROUTINE 

INITAL OPTNRD 
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ROUTINE NAME; IMAGES 


DESCRIPTION: This routine images previously input surfaces in reference 

planes that have been set up by subroutine REFGD and writes 
the surface description data for the images on the random 
access file. 


GALLING SEQUENCE: CALL IMAGES (ISAVE, NOSF, NOND, IMAG, IREF, NLEN, 

INDEXS, INDEXN, lER) 


ISAVE - 

NOSF 

NOND 

IMAG 

IREF 

NLEN 
INDEXS - 

INDEXN - 

lER 


An array of previously input surface numbers 
and their corresponding random access file 
record numbers 

Surface counter (sequence number) 

Node counter (sequence number) 

ID number of previously input surface that 
is to be imaged 

ID number of reference plane in which IMAG 
is to be imaged 

Length of a random access record 

An array of random access record numbers for 

surfaces 

An array of random access record numbers for 
surfaces 

An indicator of errors encountered 


REFERENCED BY: SEGMENT ROUTINE 

SDPSS2 SDTPS2 

FILES: NOUT - System output file 

NRIO - Random access file. 


ROUTINE NAME : INFORM 

DESCRIPTION: This routine writes out the user information contained 

in the INFO file. 

CALLING SEQUENCE: CALL INFORM 

REFERENCED BY: SEGMENT ROUTINE 

INITAL INITAL 

FILES; NINFO - Information file 

NOUT - System output file 
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ROUTINE NAME: 


INITAL 


DESCRIPTION: This routine initializes the preprocess.or'-labeled Common, 

and calls the routines to write the TEASTS banner on the 
system output file and read in the OPTION DATA block. 

CALLING SEQUENCE; GALL INITAL 

REFERENCED BY: SEGMENT ROUTINE 

start start 
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ROUTINE NAI4E : INRDB 

DESCRIPTION: This subroutine name is an entry point into routine INRDD, 

CALLING SEQUENCE: CALL INRDB 

REFERENCED BY: SEGMENT ROUTINE 

BCSRD BCSRD 

BCSPl 


ROUTINE NAME: INRDD 

DESCRIPTION: This routine reads in the NDI file that contains the user- 

input data cards. 


CALLING SEQUENCE: CALL INRDD 

CALL INRDB \ 

CALL INRDO / 

CALL INRDDC ? Entry Points 
CALL INRDSF \ 

CALL INRDOD * 

REFERENCED BY: SEGMENT ROUTINE SEGMENT ROUTINE 


I 


TRASYS 

HCARD 

SRFCRD 

SRFCRD 

QUANRD 

QUANRD 


SRFCSl 

ARRYRD 

ARRYRD 


SRFCBC 

SKIRD 

SKIRD 

CBSPRD 

CRSPRD 


KPl 




IPl 




FILES : NOLT 

NDI 


System output file 

Current model source data f"f le 


ROUTINE NAME: INRDE 

I DESCRIPTION: This routine reads in the NRSI file for the TRASYS source 

editor. 

CALLING SEQUENCE: CALL INRDE 

REFERENCED BY: SEGMENT ROUTINE 

SEDIT SEDIT 

I FILES: NRSI - Input restart data file 

NOUT - System output file 
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IRBilO ■ ■ 

RiESGl^PlliONl; . n^ie -an .entiy ,jp into rontane Il^RDD. 

REFERENCED BY: SEGMENT ROUTINE 

FFRD FFRD 

: 

:aUXEi) FLUXED 

dibedr 

LOGICS LOGICS 

ROUTINE NAME: INRDOD 

DESCRIPTION! This subroutine name is an entry point into routine INRDD. 

CALLING SEQUENCE: CALL INRDOD 

REFERENCED BY: SEGMENT ROUTINE 

LOGICl LPl 

BLDGEN 

ROUTINE NAME: IHRDPP 

DESCRIPTION: This routine reads in the MDI file for the source editor 

print /punch routine. 

CALLING SEQUENCE; CALL INRDPP 

REFERENCED BY: SEGMENT ROUTINE 

SEDIT PTPHSF 

FILES: NDI - Input source edit data file 


ROUTINE NAME: IHRDSF 

DESCRIPTION: This subroutine name is an entry point into the subroutine 

INRDD. 

CALLING SEQUENCE; CALL INRDSF 
referenced BY: SEGMENT ROUTINE 


SHDWRD 


SHDIORD 

SFTBRD 
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ROUTINE NME: IPX 

DESCRIPTION: This routine performs the first-pass processing of the 

intermediate coordinate system input cards that the user 
input in the SURFACE DATA block, 

GALLING SEQUENCE: GALL IPl (NIX) 

NIX - Index of the last word used in the dynamic 
storage array (IX) on completion of this 
routine. 

REFERENCED BY: SEGMENT ROUTINE 

SKIED SKIRD 

FILES; NOUT - System output file 


ROUTINE NAME: IP2 

DESCRIPTION; This routine performs the second-pass processing of the inter” 
mediate coordinate system input data that the user input 
in the SURFACE DATA block, 

CALLING SEQUENCE: CALL IP2 (NIX) 

NIX “ Index of the last word in dynamic storage con- 
taining the data generated and passed by IPl. 


KEY VARIABLES: NISV 

NIEV 
NISD 
NIWD 


Index of the starting word in dynamic storage that 
contains the output intermediate coordinate data 
Index of the last word in dynamic storage that 
contains the output intermediate coordinate data 
Index of the starting word in dynamic storage that 
contains the output intermediate coordinate directory 
Length of the intermediate coordinate directory 


REFERENCES: SEGMENT ROUTINE 


SKIRD SKIRD 


FILES; NOUT - System output file 
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ROUTINE NA^ : JQBHQ (GDC system routine.) 

DESGRIPTIDH; This routine returns the current job number of the job as 
the computer hnows it, , The job number is returned in dis- 
play code in the following fomat: 

UUUUUCCbb 

where UDUUU are user-assigned and cc are computer-assigned. 
CALLING SEqUMCE: CALL JOBNO (N) 

N = Returned variable for job number 
REFERENCED BY: SEGMENT ROUTINE 

INITAL INITAL 


ROUTINE NAME: RPl 


DESCRIPTION: This routine performs the first-pass processing of the 

constant input cards the user inputs in the SURFACE DATA 
blocks. 


GALLING SEQUENCE: CALL KPl (NKREC, NICNAM, NKVAL, IBUF) 


HKREC - Number of records written to the NSCl file 
NKVAL - Total number of constant values, including 
the integer count 

NKNAM - The number of constant names 
IBUF - Scratch array that may or may not be used as 

a communication link between routine KPl and 
KP2 


REFERENCED BY: SEGMENT ROUTINE 

SKIRD SKIRD 

FILES: NOUT - System output file. 

NSCl - Scratch file 1 used to communicate between routines KPl 
and KP2. 
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ROUTINE NAME: 


KPIE 


DESCRIPTION: This name is an entry point into routine KPl and is called 

to flush the IBUP array when the last constant data card 
has been processed. 

CALLING SEQUENCE: CALL KPIE (tKREC, IBUP) 

NKREC - Number of records written to the NSCl file 
IBUF “ Scratch array that may or may not be used 
as a communication link between routines 
KPl and KP2 

REFERENCED BY: SEGMENT ROUTINE 

SKIED SICIED 

FILES: NSGl ~ Ccratch file 1 used to communicate between routines 

B.P1 and KP2 


ROUTINE NAME: KP2 

DESCRIPTION: This routine performs the second-pass processing of the 

constant input cards the user inputs in the SURFACE DATA 
blocks. 

I CALLING SEQUENCE: CALL KP2 (NKREC, NKNAM, NKVAL, IBUF) 

Argument description is same as for routine KPl, 


KEY VARIABLES: 

NIiSD - 

Index of the starting word in dynamic storage 
that contains the output constant data 


NKWD . “ 

Length of the constants directory 


KKSV - 

Index of the starting word in dynamic storage 
that contains the output constant values 


NKEV - 

Index of the last word in dynamic storage that 
contains the output constant values 

REFERENCED BY: 

SEGMENT 

ROUTINE 


SKIRD 

SKIED 


FILES: NOUT - System output file 

NSCl - Scratch file 1 used to communicate between routines KPl 
and KP2 
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ROUTINE NAME: LFILL 

DESCRIPTION: This function routine left- justifies a word that is filled 

with zeros to the left of the data. 

CALLING SEQUENCE: I = LFILL (IN, IFILL) 

IN - Word to be left-justified 
IPILL - The character to fill the word with right 
of the data 


REFERENCED BY; 

SEGMENT 

ROUTINE 

SEGMENT 

ROUTINE 


TRASYS 


QUANED 

QUANRD 



CDPRCl 

ARRYRD 

ARRYRD 



PTS^ER 

SKIRD 

SKIRD 


INITAL 

OPTNRD 


IPl 




SRFCRD 

SRFCSl 





GDPRC3 




BCSRD 

BCSRD 



H10T6 


BCSPl 




SHDWRD 

SHDWRD 




LOGICl 

LOGICl 


SEDIT 

SEDIT 


BLDGEN 



MERGE 

LOGICS 

LOGICS 


CRSPRD 

CRSPRD 




FFRD 

FFRD 



ROUTINE NAME: 

LLSF (CDC 

system routine) 




DESCRIPTION: This routine is a function that left-shifts a word with 

wrap-a-round on the pushoff bits, It is usually a system 
routine. 

calling SEQUENCE: I = LLSF (IWORD, IB ITS) 

IWORD - Word to be left-shifted 
IBITS - Number of bits to shift 

REFERENCED BY: 


SEGMENT 

ROUTINE 

SEGMENT 

ROUTINE 

SEGMENT 

ROUTINE 

TRASYS 

HCARD 

SEDIT 

SEDIT 

BCSRD 

BCSRD 


INRDD 




BCSPl 


CDPRCl 




BCSP2 


LFILL 


ASKCRD 

DCMNRD 

DCMNRD 


PTSVER 


PTPHSF 

FEMFRD 

FRMFED 


BSRGHD 

QUANRD 

QUANRD 


READVF 

START 

HEADCD 

ARRYRD 

ARRYRD 

SHDWRD 

SHDWRD 

INITAL 

OPTNRD 

SKIRD 

SKIRD 


SFTBRD 


HTI 


KPl 

FLUXRD 

FLUXRD 


H10T6 

SRFCRD 

SRFCRD 


DIBLDR 


OPTNCV 


SRFCSl 

CRSPRD 

CRSPRD 

MCOLL 

MCOLL 


CDPRCS 

LOGICl 

LPl 


MCOLCV 

SDPSS2 

SDTPS2 

LOGIC2 

L0GIC2 




IMAGES 

LOGICS 

LOGICS 
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ROUTINE NAME: LNGCOR ' 

DESCRIPTION: This rdutine prints the core lengths of the processor- 

segments as calculated by the preprocessor. 

CALLING SEOTIENGE; CAI;L ; 

REFERENCED BY: SEGMENT RQUTII^ 

LOGICS LOGICS 


ROUTINE NAME: LQCF (GDC system routine) 

DESCRIPTION: This routine is a function that returns the absolute 

memory word address of the argument passed to it. 

The routine is usually a Fortran-callable system routine. 

CALLING SEQUENCE; I = LOCF (j) 

J = The word whose address is wanted 


REFERENCED BY: 


SEGMENT 

ROUTINE 

TRASYS 

. NS GO PE 


DDUMP 


BtlFCOR 

start 

ALTFILE 

inital 

INITAL 

SKIED 

IP2 


RP2 

SDPSS2 

SDPSS2 

BCSRD 

BGSP2 

EAPUP 

RAPUP 


ROUTINE NAME: LOGICO 

DESCRIPTION: This routine is the controlling driver for the processing 

of the user's logic block, input 


CALLING SEQUENCE: CALL LOGICO 


REFERENCED BY : SEGMENT ROUTINE 


TEASYS . TRASYS 

FILES: NOUT - System, output file 

REV 1 
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EDUTiNE KAtffi: LQGICl 

DESGRIPTIONi: This routine and the routines that it calls read in the 

array, quantities, and BCS directories. It then checks 
for duplicate names between the directories ' and reserve 
variable list and creates the conmvon array consisting of 
names that will be used to write the processor common 
arrays. This routine also reads in the user's OPERATION 
DAtA block and counts the step cards and the substep cards 
within a step. This information will be used to create 
the computed GO TO statements in routine L06IC2, 

GALLING SEQUENCE: GALL LOGICl 

REFERENCED BY: SEGMENT ROUTINE 

lOGICO LOGlCO 

FILES: NOUT -System output file 

NRIOZ - Preprocessor and processor random I/O communication file 


ROUTINE NAME: L0GIC2 

DESCRIPTION: This routine generates the main program segment of the pro- 

cessor and also reads in the OPERATION DATA output from 
rou ne LPl and processes it into valid Fortran routines. 
Thes routines are written to file NCMPL. 


CALLING SEQUENCE: CALL L0GIC2 

REFERENCED BY: SEGMENT ROUTINE 

LOGICO LOGICO 


FILES : NOUT - 

NSGl “ 

NCMPL - 


System output file 

Scratch file containing user’s operation data output by 
routine LPl 

The file that the generated processor Fortran cards are 
written to 


ROUTINE NAME: 


LOGICS 


DESCRIPTION: This routine and the routines that it calls read in and 

process the user's SUBROUTINE DATA block. This routine 
also writes all of the primary routines of all the seg- 
ments, along with the user-input subroutines to the NCMPL 
f i].e. 

CALLING SEQUENCE: CALL LOGICS 

REFERENCED BY: SEGMENT ROUTINE 

LOGICO LOGICO 

FILES: NOUT - System output file 

NCMPL - The file that the generated processor Fortran cards are 
written to 


ROUTINE NAME : LPl 

DESCRIPTION: This routine reads the user's OPERATION DATA block in pass 1 

for the number of steps and number of substeps within a step. 
Th_e ORJBGEN user-input card is also expanded in this routine 
to standard type of input. The step directory is written to 
file KSQNTL. 

CALLING SEQUENCE: CALL LPl 

REFERENCED BY: SEGMENT ROUTINE 

LOGICl LOGICl 

FILES: NOUT -System output file 

NSCl -Scratch file which the operation data block cards are 
written to for pass 2 processing 

NSQNTL-Preprocessor and processor sequential l/O communication file 
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ROUTINE NMIE: 


LRSF 


i 

\ 

DESCRIPTION: This routine is a function that right-shifts a word 

with sign extension. This routine is usually a system 
function. ‘ 

CALLING SEQUENCE; I = LRSF (IWORD, TSiTS) 

IWORD = Word to be right-shifted 
IBITS = Number of bits to shift 



REFERENCED BY: 


SEGMENT 

ROUTINE 

TRASYS 

INITAL 

BSRCHD 

OPTNRD 

QUAKRD 

QUANRD 

ARRYRD 

ARRYRD 

SKIRD 

KPl 

SRFCRD 

SRFCSl 

BCSRD 

BCSRD 

SHDNRD 

SFUNPK 

FLUXRD 

SFTBRD 

DIBLDR 

SEDIT 

MERGCD 

FFRD 

MERGE 

FFPROC 

L0GIC2 

FFEEAD 
FFSORT 
LOGIC 2 


e 


\ 
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ROUTIHE name: MERGCD 


DESCRIPTION: This routine processes the CMERG source editor directives. 

Once the CMERG card is decoded, the NCMG file is searched 
to find the requested file and cards , and the requested 
cards are then written to the NDI file. 

I CALLING SEQUENCE: CALL MERGCD (NCARD) 

NCARD = Last new edit sequence number used 

REFERENCED BY: SEGMENT ROUTINE 

SEDIT SEDIT 

FILES: NOUT - System output file 

NCMG - User-input card image supplementary input file 


ROUTINE NAME; MERGE 

DESCRIPTION; This routine processes the EMERG source editor directives. 
Once the EMERG card is decoded, the NEMG file is searched 
for the requested model and cards, and the requested cards 
are then written to the NDI file. 

CALLING SEQUENCE; CALL MERGE (NCARD) 

NCARD = Last new edit sequence number used 

REFERENCED BY: SEGMENT ROUTINE 

SEDIT SEDIT 

FILES; NOUT - System output file 

NEMG - User-input source edit-formatted supplementary input file 
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ROTJTIHE NAME: 


NONDP 


DESCRIPTION: This function routine, when given the internal 

sequence number of the node, returns the node 
number . 

CALLING SEQUENCE: I = NONDF (NPf) 

REFERENCED BY: SEGMENT ROUTINE 

SDFSS2 SDTPS2 
IMAGES 
AUTOCM 
STUFDT 


ROUTINE NAME: OPENMS (CDC system routine) 

DESCRIPTION: This routine opens a mass storage random I/O file. 

CALLING SEQUENCE: CALL OPENMS (U, IX, LWGTH, T) 


U - Unit designator 

IX - First word address in central memory of 

the array that will contain the file index 
LNGTH - Length of the index 

LNGTH “ (no, of records +1) for number index. 
LNGTH - 2 * (no, of records -fl) for name 
index 

X _ X = 0 file is referenced through a number 

master index 

T = 1 file is referenced through a name 
master index (TRASYS uses only a numbered 
index) 


referenced BY: SEGMENT 


ROUTINE 


INITAL INITAL 

SRFCKD SRFCRD 

shdwrd shdwrd 
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ROUTINE NAME: 


OPTNCV 


DESCRIPTION: This routine decodes the user's input cards that are input 

in the OPTION DATA block. 

CALLING SEQUENCE: CALL OPTNCV 

REFERENCED BY: SEGMEN T ROUTINE 

INITAL OPTNRD 

ROUTINE NAME: OPTNRD 

DESCRIPTION; This routine reads in and processes the user-input OPTION 
DATA block. 

CALLING SEQUENCE: CALL OPTNRD (I 1ST) 

REFERENCED BY: SEGMENT ROUTINE 

INITAL INITAL 

FILES: NOUT - System output file 

NIN - System input file 

ROUTINE NAME: ORBGEN 

DESCRIPTION: This routine and the routines that it calls process the 

user's logic input ORBGEN card, 

CALLING SEQUENCE: CALL ORBGEN 

REFERENCED BY; SEGMENT ROUTINE 

LOGICl LPl 

FILES; NOUT - System output file 

NSCl - Scratch file that the output of logic pass 1 is written to 
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ROUTINE NAME: ORBIH 


) 


DESCRIPTION: This routine writes the operations data block code for 

automatic generation of inertial (sun or star) oriented 
orbits. 

CALLING SEQUENCE: CALL QRBIN (IREP, TRU, DVC, NPOINT, ILORC) 


IREP 

TRU 

DV6 

NPOINT - 
ILORC - 


Repeat flag for complete orbits 

Initial true anomaly 

True anomaly increment 

Number of orbit points generated 

Flag to eliminate AQCAL calls from operations 

data block code 


REFERENCED BY: SEGMENT 


ROUTINE 


LOGICl ORB GEN 


FILES : 


NOUT - System output file 

NSCl - Scratch file that the output of logic pass 1 is written to 


ROUTINE NAME: ORENPL 

DESCRIPTION: This routine writes the operations data block code for auto- 

matic generation of heliocentric orbits. 

galling SEQUENCE; CALL ORBNPL (IREP, TRU, DV&, NPOINT, ILORC) 

IREP - Repeat flag for complete orbits 
TRU - Initial true anomaly 

DV6 - True anomaly increment 

NPOINT - Number of orbit points generated 
ILORC “ Flag to eliminate AQCAL calls from operations 
data block code 

REFERENCED BY: SEGMENT ROUTINE 

LOGICl ORBGEN 

FILES: NOUT - System output file 

NSCl - Scratch file that the output of logic pass 1 is written to 
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ROUTINE NAME: 


ORBPL 


DESCRIPTION: This routine writes the operations data block code for 

automatic generation of noncircular, planet-oriented orbits, 

•CALLING SEQUENCE: GALL ORBPL (IREP, TRU, DV6, WPOINT, ILORC) 



IREP \ 
TRU j 

DV6 \ 

NPOINT 1 
ILORC ! 

Reference ORBIN 

REFERENCED BY: 

SEGMENT 

ROUTINE 


LOGICI 

ORBGEN 

FILES: NOUT - 

System output file 


NSCl - Scratch file that the output of logic pass 1 is 
written to 


ROUTINE NAME: ORBPLC 

DESCRIPTION: This routine writes the operations data block code for 

r automatic generation of circular, planet-oriented orbits. 

CALLING SEQUENCE: CALL ORPLC (IREP, TRU, DV6, NPOINT, ILORC) 

IREP 
TRU 
DV6 

NPOINT 
ILORC 

REFERENCED BY: SEGMENT ROUTINE 

LOGICI ORBGEN 

FILES: NOUT - System output file 

NSCl - Scratch file that the output of logic pass 1 is written to 


1 

> Reference ORBIN 

) 
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ROUTINE NAME 


: OTWTE 

DESCRIPTION; This routine writes the NDI file called by the TRASYS 
source editor. 

galling SEQUENCE; CALL OTWTE (IILUSH) 

lELDSH =0 Not the last call to this routine 
IFLDSH = 1 The last call to this routine. Flush the 
lOB array to the NDI file 

KEY VARIABLES; lOB — An array used for collecting card images until 

lOT words have been filled 

lOT - Last word to fill in the lOB array before flush- 
ing the data to the NDI file 

REFERENCED BY; SEGMENT ROUTINE 

SEDIT SEDIT 

MERGE 
MERGCD 

FILES: NDI - Source edit-formatted output file used as input file for 

data and logic preprocessor segments 


ROUTINE NAME; PAGEH 

DESCRIPTION: This routine ejects the page and writes the page heading 

on all pages generated by the preprocessor, 

galling SEQUENCE: CALL PAGEH (N) 


N - 0 Do not write card column designator line 
N = 1 Write card column designator line 


REFERENCED BY; SEGMENT 

ROUTINE 

SEGMENT 

ROUTINE 

TRASYS 

INRDD 

SRFGRD 

SRFCSl 


SERROR 


SRFCBC 

START 

HEADCD 

BCSRD 

BCSRD 

INITAL 

OPTNRD 

DCMNRD 

DCMNRD 


INFORM 

FFRD 

FFRD 

SEDIT 

SEDIT 

SHD^ffiD 

SHDWRD 



FLUXRD 

FLUXRD 


MERGE 

CRSPRD 

CRSPRD 


PTPHSF 

LOGICl 

LOGICl 


MERGCD 

L0GIC2 

LOGIC2 


PTHSTY 

LOGICS 

LOGICS 

LNGGOR 

DATARD 

DATARD 

RAPUP 

RAPUP 

QUANRD 

QUANRD 



ARRYRD 

ARRYRD 



SKIRD 

SKIRD 

PTCKSl 




NOUT “ System output file 
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FILES ; 


ROUTINE NAME: PARAB 


DESCRIPTION: This routine converts point input for paraboloxj^:: to program- 

compatible surface description parameters and se;s up Euler 
angles and a position vector to transform the luS, BCS, or CCS 
into the SCS of the paraboloid. 

GALLING SEQUENCE; GALL PARAB 

REFERENCED BY: SEGMENT ROUTINE 

SDPSS2 SDTPS2 

FILES: NOUT - System output file 


ROUTINE NAME: FDIMP (GDC system routine) 

DESCRIPTION; This routine dumps the main memory on the system output 
file. 

calling SEQUENCE: CALL PDUMP (A, B, C) 

A - Starting location to start dumping 
B - Last word to be dumped 
C “ Mode in which to dump the words 

0 = Octal dump 

1 = Real dump 

2 “ Integer dump 

3 = Octal dump 

4 = Octal dump 

REFERENCED BY; SEGMENT ROUTINE 


SDPSS2 

SDTPS2 


IMAGES 

FRMFRD 

FRMFRD 

BLDR 

FINAL 

MERGF 

FLUXED 

FLUXED 

DIBLDR 


FILES : NOUT - System output file 
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ROUTINE NAME: POLYGN 


DESCRIPTION: This routine generates N - 2 triangles from point input, where 

N is the numher of points , and sets up Euler angles and a posi- 
tion vector for each triangle generated, 

CALLING SEQUENCE: C^L POLYGN (ISUEP) 

ISUEF - A counter, from 1 to N-'2, indicating which 
triangle is being generated 

REFERENCED BY: SEGMENT ROUTINE 

SDPSS2 SDTPS2 

FILES: NOUT - System output file 


ROUTINE NAME: POSIT 


DESCRIPTION: This routine transposes a vector in the IGS, BCS, or CCS 

into the SCS of a reference plane (for imaging purposes) , 
negates the Z component, and transforms it back into the 
IGS, BCS, or CCS. 

CALLING SEQUENCE: CALL POSIT (X, Y, Z, TRAN) 


X, Y, Z - Vector components in the ICS, BCS, or CCS 
TRAN - A 3x3 matrix of direction cosines to 

transform a vector from the IGS, BCS, or 
CCS to the SCS 


REFERENCED BY: SEGMENT ROUTINE 

SDPSS2 IMAGES 
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RODTINE NAME: 


PPTIM (CDC system routine) 


DESCRIPTION; This routine will return the accumulated peripheral 
processor time, in integer seconds, incurred up to 
the time of the call. 

CALLING SEQUENCE: CALL PPTIM (N) 

N = Return integer preprocessor time 

REFERENCED BY; SEGMENT ROUTINE 


SDPSS2 

S DTPS 2 

SDPSS2 

IMAGES 

TRASYS 

TRASYS 

START 

START 

DATARD 

DATARD 

LOGICO 

LOGICO 

TPGEN 

TPGEN 

RAPUP 

RAPUP 


ROUTINE NAME: PRNTCK 

DESCRIPTION: This routine outputs the traced node/surface data if 

the variable ITRC25 is set in the OPTIONS DATA block. 

CALLING SEQUENCE; CALL PRNTCK (ITP) 

ITP - Flag defining node dump or surface dump 

REFERENCED BY: SEGMENT ROUTINE 

SDPSS2 SDTPS2 

IMAGES 

FILES: NOUT - System output file 
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EDUTIISE NAME: FTCHEK 


DESGSIPTION: This routine, given an array of points, checks for duplicates, 

and returns ICHEK as a flag, 

CALLING SEQUENCE: CALL PTGHEK (N, IT, JT, X, ICHEK) 

N " Number of points to be checked 

IT - Skip this point 

JT “ Skip this point 

X - Array of point data 

ICHEK - 0 = Paints not same 

1 = Points same 

REFERENCED BY: SEGMENT ROUTINE 

SDPSS2 BOX 

POLYGN 

TRAPZ 

RECT 


ROUTINE NAME: PTCKSl 

DESCRIPTION: This routine writes out the surface ICS directory, ICS 

values, constant-directory and the constant values. 

This routine is only executed when the surface data pass 
1 error trace flag is turned on in the OPTIONS DATA block. 

CALLING SEQUENCE: CALL PTCKSl 

REFERENCED BY: SEGMENT ROUTINE 

SKIRD SKIRD 

FILES: NOUT - System output file 
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ROl 


PTDIR 


DESCR, '’TION ; This routine writes the directory record of a source 
edit-"forinatted tape to the system output file. 

CALLING SEQUEWgE; CALL PTDIR (NU, IDIR) 

KfTJ " Pile that contained the directory to be 
printed ■ 

IDIR “ The array containing the directory to be 
printed 

FILES: ROUT ~ System output file 


ROUTINE NAME : PTHSTY 

DESCRIPTION: This routine writes the history record of a source 

edit-formatted tape to the system output file. 

CALLING SEQUENCE: CALL PTHSTY (IHS, ITTPFT) 

IHS - The array that contains the history 
record to be printed 

ITTPFT - Approximate tape footage used by model 
being printed 

REFERENCED BY; SEGMENT ROUTINE 

GNEDO GNEDO 

FILES: NOUT - System output file 
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i 


I 

j 

ROUTINE NAME: PTPHSF j 

DESCRIPTION: This routine prints and/or punches the NDI source edit- 

formatted file. 

j 

CALLING SEQUENCES: CALL PTPHSF (ISEQ) j 

ISEQ =0 No sequencing in Columns 73-80 is wanted 

on the punched output ! 

ISEQ ^ 0 Sequencing in Column 73-80 is wanted on 
the punched output 

REFERENCED BY: SEGMENT ROUTINE 

SEDIT SEDIT 

FILES: NOUT - System output file 


ROUTINE NAME: PTSVER 

DESCRIPTION: This routine writes the character mark under the bad 

character or field when an error is encountered. This 
routine also saves the bad card edit sequence line number 
so that it can be printed at the end of the preprocessor 
execution. 

CALLING SEQUENCE: PTSVER (IC, lER) 

IC “ Column to print error character mark 

IC = 1 is actual character 7 on the input card. 
lER - =0 Fatal error was encountered. 

7 ^ 0 Caution error was encountered. 

REFERENCED BY: SEGMENT ROUTINE 

TRASYS SERROR 

FILES: NOUT - System output file 
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ROUTINE NAME : QUANRD 


DESCRIPTION: This routine processes the user-input QUANTITIES DATA 

block. 

CALLING SEQUENCE; CALL QUANRD 

REFERENCED BY: SEGMENT ROUTINE 

DATARD DATARD 

FILES: NOUT - System output file 

NSGl - Scratch file used to write the user quantity values. 

This file is written in pass 1 and read in pass 2 
NRIO - Preprocessor and processor random I/O communication 
file 


ROUTINE NAME: QUANSD 

DESCRIPTION: This routine writes to the NRIO file the control quantities 

directory and values when the user does not input a 
QUANTITIES DATA block. 

CALLING SEQUENCE; CALL QUANSD 

REFERENCED BY: S EGMENT ROUTINE 

QUANRD QUANRD 

FILES: NOUT - System output file 

NRIO - Preprocessor and processor random I/O communication file 


ROUTINE NAME: RAP UP 

DESCRIPTION: This routine writes the time and core statistics at the 

end of the preprocessor output. This routine also terminates 
the preprocessor either normally or in the case of fatal 
errors, abnormally. 

CALLING SEQUENCE: CALL RAPUP 

REFERENCED BY; SEGMENT ROITIINE 

TRASYS TRASYS 

FILES: NOUT - System output file 
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ROUTINE NAME: RDPSN 


DESCRIPTION: This routine, given the dimension of a node, defines the 

radius of a sphere Large enough to encompass the node and 
the position vector to the center of the sphere. 

calling SEQUENCE: CALL RDPSN (VEC, BETA, GAMMA, DB, DG) 

VEC(3) - Array containing (X, Y, Z) position to center 
BETA, gamma - Length along side of node 
DB, DG - Distance to center of node 

REFERENCED BY: SEGMENT ROUTINE 

SDPSS2 SDTPS2 

ROUTINE NAME: RDPSS 

DESCRIPTION: This routine computes the radius and corresponding position 

vector used to encompass the surface and locate the center. 

calling SEQUENCE: CALL RDPSS (VEC) 

VEC - Array of position vectors 

REFERENCED BY; SEGMENT ROUTINE 

SDPSS2 SDTPS2 


ROUTINE NAME: READMS 

DESCRIPTION: This routine reads a logical record from a random I/O 

mass storage file. 

CALLING SEQUENCE; CALL READMS (U, FWA, N, K) 

U - Unit designator. 

FWA - Address in central memory of first word of record 
N - Number of central memory words in the record to 
be transferred 

K - Number index of record or name index of record to 
be read. (In TRASYS, K = a number index always) 

REFERENCE..' BY: SEGMENT ROUTINE 

SRFCRD DUPSRF 

SDPSS2 IMAGES 

BCSRD BCSRD 

SHDWRD SFWPL 

LOGICl LOGICl 
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ROUTINE NAME: 


READVF 


DESCRIPTION: This routine reads in header form factor data and loads the 

various arrays generating all requests for data. 

CALLING SEQUENCE: CALL READVP (NODINl* JEOF, INDEX, NN2, NN, IFRST, NMAX) 


NODINl 

JEOF 

INDEX 

NN2 

NN 

IFRST 

NMAX 


Array of packed node numbers and values 
= 3HEND End of data, 9^ 3HEND Not end of data 
Number of values stored in NODINl 
Number of nodes * 2 
Number of nodes 

Flag indicating first time called 
Maximum number of words in data block 


REFERENCED BY; SEGMENT ROUTINE 


FRMFRD BLDR 


FILES: NOUT - System output file 


ROUTINE NAME: SECT 

DESCRIPTION: This routine converts point input for rectangles to program- 

compatible surface description parameters and sets up Euler 
angles and a position vector to transform the ICS, BCS, or 
CCS into the SCS of the rectangle. 

CALLING SEQUENCE: CALL EECT 

REFERENCED BY: SEGMENT ROUTINE 

SDPSS2 SDTPS2 

FILES: NOUT - System output file 


ROUTINE NAME: REFCD 

DESCRIPTION: This routine sets up surface description information for 

reference planes and writes it on the random access file. 

CALLING SEQUENCE: CALL REFCD (ISSX) 

ISSX " The index of the starting location where surface 
description data for the reference plane is to 
be loaded in blank common 

REFERENCED BY: SEGMENT ROUTINE 

SRFCRD SRFCSl 

FILES: NRIO - Random access file 
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ROUTINE NAME: RETUBNS (GDC system function) 


DESCRIPTION: This routine releases a local file from a job. 

CALLING SEQUENCE: CALL RETURNS (N) 

N = Name of file to be released 


REFERENCED BY: SEGMENT ROUTINE 

SRFCRD SRFCRD 


ROUTINE NAME : S CHECK 


DESCRIPTION: This routine checks the validity of user- input surface de- 

scription parameters and sets the default values, 

CALLING SEQUENCE: CALL SGHECK (ISSX, lERR) 


ISSX - 
IERR - 

The index of the starting 
description data in blank 
Counter for the number of 
in the user- input surface 

location of surface 
common 

fatal errors found 
descriptions 

REFERENCED BY: SEGMENT 

ROUTINE 


SRFCRD 

SRFCSl 



FILES; NOUT - System output file 


ROUTINE NAME: SDPSS2 

DESCRIPTION: This routine sets up blank common indexes for use in pro- 

cessing the surface description data in pass 2. 

CALLING SEQUENCE: CALL SDPSS2 

REFERENCED BY: SEGMENT ROUTINE 

DATARD DATARD 

FILES: NOUT - System output file 
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ROUTINE NAME; SDTPS2 

DESCRIPTION; This routine reads user- input surface description data from 
a sequential scratch file that was written in surface data 
pass 1, converts the data to program compatible surface 
description parameters, and writes the data on the random 
access file. 

CALLING SEQUENCE; CALL SDTPS2 (INDEXS, INDEXN, XSAVE) 

INDEXS - Array of BCS names and corresponding random 
access record numbers for surfaces 

INDEXN - Array of random access record numbers for 
nodes 

ISAVE - Array of surface numbers and corresponding 
random access record numbers 

REFERENCED BY: SEGMENT ROUTINE 

SDPSS2 SDPSS2 

FILES: NOUT - System output file 
NRIO - Random access file 
NSGEl - Scratch file 


ROUTINE NAME: SECOND 

DESCRIPTION: This routine returns the central-processor time from the 

start of job, in seconds, in floating-point format, accurate 
to one thousandth of a second 

calling SEQUENCE: CALL SECOND (T) 

T = Variable that central processor seconds will be retumed to 

REFERENCED BY; SEGMENT ROUTINE 

TRASYS TRASYS 

START START 

DATARD datard 

LOG I CO LOG I CO 

TPGEN TPGEN 

RAPUP RAPUP 
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ROUTINE NAME: SEDIT 


DESCRIPTION: This routine and the routines that it calls perfprm the 

preprocessor source edit function of the TRASYS program, 

CALLING SEQUENCE: CALL SEDIT 

REFERENCED BY: SEGMENT ROUTINE 

START START 

FILES: NIN - System input file 

NOUT “ System output file 


ROUTINE NAME; S ERROR 


DESCRIPTION: This routine is called when a caution or error condition 

results when preprocessing the user's input. This routine 
and the routines that it calls store the error accounting 
information and account for proper line-page format on 
printed error messages. 

CALLING SEQUENCE; CALL SERROR (NTYP, NARRW, ILINE, NCOL, KER) 


REFERENCED BY: 


NTYP = 0 Caution message is to be printed 
= 1 Error message is to be printed 
NARRW - Column to print error character 
ILINE - Number of error message lines to be printed 
NCOL =0 No output column number on page heading 

= 1 Print output column numbers on page heading 


KER 

- The returned 

caution 

or error number 

SEGMENT 

ROUTINE 

SEGMENT 

ROUTINE 

TRASYS 

HCARD 

SDPSS2 

BOX 


CDPRCl 


CALB 


GDPRC2 


galg 


FDPRC 


CHEG 

START 

HEADCD 


CONE 

inital 

OPTNRD 


cylndr 

SEDIT 

SEDIT 


DISC 

IMAGES 

PARAB 


EPTAPE 


polygn 


ASKCRD 


TRAPZ 


INRDE 


RECT 


MERGG 


SPHERE 

datard 

DATARD 

STORDT 

BCSRD 

BCSRD 

quanrd 

QUANRD 


BCSPl 

ARRYRD 

ARRYRD 

FRMFRD 

BCSP2 

FRMFRD 
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REFERENCED BY; (continued) 


SEGMENT 

ROUTINE 

SEGMENT 

ROUTINE 

SKIRD 

SKIRD 


FINAL 

READVF 


KPl 

SHDWRD 

SHDWRD 


I PI 


SFTBRD 


IP2 




KP2 



SRFCRD 

SRFCRD 

FLUXRD 

FLUXED 


SRFCSl 


DIBLDR 


SRFCBC 

CDPRC3 

CRSPRD 

CRSPRD 


DUPSRF 

LOGIGO 

LOGIGO 


SGHEGK 

LOGIGl 

LOGIGl 

LPl 


SRFCDM 

LOGIC2 

ORBGEN 

LOGIC2 



LOGICS 

LOGICS 


FILES; NOUT - System output file. 


ROUTINE NAME: SFCIQ 

DESCRIPTION: This routine copies the unused or requested SAVE models 

on the NSHADI file to the NSHADO file. 

CALLING SEQUENCE: CALL SFCIO (ISHAD2, NSF) 

ISHAD2 = Shadow factor request directory 

NSF = Number of shadow factor models on NSHADI file. 

REFERENCED BY: SEGMENT ROUTINE 

SHDWRD SHDWRD 

FILES: NSHADI ~ Shadow factor use- input file 

NSHADO - Shadow factor user-output file 
NOUT - System output file 
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ROUTINE NMIE: 


SFRSI 


DESCRIPTION: This routine reads in the requested shadow factor data from 

the NSHADI file and writes the data to the NPLSR file, 

CALLING SEQUENCE; CALL SFRSI (NDEl, NDE2, NDE3, NDE4, lECC, IPCC, NODEA, 

NODET) 

NDEl “ Last word of the IX array that contains the 
node array directory input by the user 
NDE2 - Last w6rd of the IX array that contains the 

node array directory input from the NSHADI file 
NDE3 “ Last word of the IX array that contains the 
shadow factor request directory 
NDE4 - Last word of the IX array that contains the 
scratch random I/O file index array (used on 
CDC systems only) 

lECC - Output array from this routine that will contain 
the unpacked cone-clock values 
IPCC - Array containing the cone-clock values in packed 
form read from the NSHADI file in this routine 
NODEA - Number of nodes in the node array input from 
user input data 

NODET - Nuirfser of nodes in the node array input from 
the NSHADI file 

ICEY VARIABLES: IX - Preprocessor blank common dynamic storage array 

REFERENCED BY: SEGMENT ROUTINE 

SHDWRD SHDWRD 

FILES: NOUT “ System output file 

NSHADI - Shadow factor user input file 
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ROUTINE NAME: 5ETBRD 


DESCRIPTION: This routine reads and processes the TABLE and REC0MP 

shadow factor data block input cards for a. single node and 
writes the processed data to the NRIOS file. 

calling SEQUENCE: CALL SFTBRD (NDEl, NDS3 , lECC) 

NDEl " Last word of the IX array that contains the 
node array directory input by the ujer 

NDS3 - Starting address of the shadow factor request 
directory 

lEGC - Scratch array to contain the cone- clock value 
for writeout to the NRIOS file 

REFERENCED BY: SEGMENT ROUTINE 

SHDWRD SHDWRD 

FILES: NOUT - System output file 

NRIOS - Scratch random I/O file that the user-input shadow 
factor data are written to 
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ROUTINE NAME: 


SFWFL 


DESCRIPTION; This routine writes the NPLSR file from the shadow factor 
data that reside on the WRIOS file. 


CALLING SEQUENCE: CALL SFWPL (NDEl, NDS3, lECC, KECC) 


NDEl - 

Last word of the IX array that contains the 
node array directory input by the user 
This array is also the node output directory 

NDS3 - 

Starting word of the IX array that contains 
the shadow factor request directory 

lECC - 

Scratch array to read in and write out 

of 

KECC " 

Length of the records to be read off the 
NRIOS file 

REFERENCED BY: SEGMENT 

ROUTINE 

SHDWRD 

SHDWRD 


FILES: NOUT - System output file 

NRIOS - Scratch random I/O file 

NPLSR - File used to transmit the shadow factor data from the pre-* 
processor to the processor 


ROUTINE NAME: SHDWRD 

DESCRIPTION: This routine and the routines that it calls read in and 

process the SHADOW DATA user’s input block, combine it 
with the shadow factor data residing an a shadow factor 
input tape, and output the data on a file to be processed 
in the processor phase. 

CALLING SEQUENCE: CALL SHDWRD 

REFERENCED BY: SEGMENT ROUTINE 


DATARD DATARD 


FILES: 


NOUT - System output file 

NPLSR - File used to transmit the shadow factor data from the 
preprocessor to the processor 
NSHADI- The user-input shadow factor input tape 
NRIOS “ Scratch random I/O file 
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ROUTINE NAME: SKRILE 


DESCRIPTION: This routine skips to and positions after an end-of-file 

mark on a specified mass storage unit. 

calling SEQUENCE; GALL SKPILE (U, N) 

U - Unit name or number that is to have the file 
skipped 

N - The number of files to skip 


REFERENCED BY: SEGMENT ROUTINE 


SEDIT EPTAPE 

GNEDO GNEDO 

DATARD STORDT 

SHDWRD SHOWED 

SFCIO 


ROUTINE NAME: SKIRD 

DESCRIPTION: This routine, along with the routines that it calls, reads 

in and processes the I and K cards that the user inputs in 
the SURFACE DATA input block. 

CALLING SEQUENCE: CALL SKIRD 

REFERENCED BY: SEGMENT ROUTINE 

DATARD DATARD 

FILES; NOUT - System output file 


ROUTINE NAME; SQRTD 

DESCRIPTION: This routine numerically sorts a doublet array. The original 

order is not preserved on equal comparisons. 

GALLING SEQUENCE: GALL SORTD (lA, NA) 


lA - 

The array to be sorted. The sorted array is 


also returned 

in this array 

NA - 

The length of 

the array to be sorted 

REFERENCED BY: SEGMENT 

ROUTINE 


INITAL 

OPTNRD 


QUANRD 

QUANRD 


ARRYRD 

ARRYRD 


SKIRD 

SKIRD 


SRFCRD 

SRFCSl 
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REFERENCED BY: 


(continued) 


SEGMENT 

ROUTINE 

BCSRD 

BCSRD 

FRMFRD 

FRMFRD 

SHDWRD 

SHDWRD 

SFRSI 

FLUXRD 

FLUXRD 

LOGICl 

LOGICl 

LPl 




ROUTINE NAME: SORTDB 

DESCRIPTION; This routine numerically sorts a doublet array. The origi- 
nal order is preserved on equal comparisons. 

GALLING SEQUENCE; CALL SORTDB (lA, IB, NN) 


lA - The array to be sorted. The sorted array is 
also returned in the array 

IB - A scratch array equal in length to the lA array 
NN - The length of the array to be sorted 


REFERENCED BY: SEGMENT ROUTINE 


FRMFRD BLDR 


ROUTINE NAME: SORTS 

DESCRIPTION: This routine numerically sorts an array. The original order 

is not preserved on equal comparisons, 

GALLING SEQUENCE: CALL SORTS (A, JJ) 


A - The array to be sorted- The sorted array is 
also returned in this array 
JJ - Length of the array to be sorted 


REFERENCED BY: SEGMENT ROUTINE 


BCSRD BCSRD 
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ROUTINE NAME: 


SPHERE 


DESCRIPTION: Converts point input for spheres to program-compatible 

surface description parameters and sets up Euler angles 
and a position vector to transform the ICS, BCS, or CCS 
into the SCS of the sphere. 

CALLING SEQUENCE: CALL SPHERE 

REFERENCED BY: SEGMENT ROUTINE 

SDPSS2 SDTPS2 

FILES: NOUT - System output file 

ROUTINE NAME: SRFCBC 

DESCRIPTION: This routine processes the BCS user-input cards input in 

the BCS DATA block. 

CALLING SEQUENCE: CALL SRFCBC 

REFERENCED BY: SEGMENT ROUTINE 

SRFCRD SRPCRD 

FILES: NOUT - System output file 


ROUTINE NAME: SRPCDM 


DESCRIPTION; This routine processes the DIMENSION variable and its 
associated data values that the user inputs as part of 
the surface description input. 

CALLING SEQUENCE; CALL SRFCDM (KTYPE, IBM, JDM) 


KTYPE 


IDM 

JDM 


The type of surface being processed: 


1 = Rectangle 6 = 

2 = Disk 7 = 

3 = Trapezoid 8 = 

4 = Cylinder 9 = 

5 = Cone 10 = 


Length of the JDM array 
Array containing the data 


Sphere 

Paraboloid 

5- Sided box 

6- Sided box 
Polygon 

to be processed 


REFERENCED BY: SEGMENT ROUTINE 


SRFCRD SRFCSl 


FILES; NOUT - System output file 
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ROUTINE NAME: 


SRFCRO 


DESCRIPTION: This routine and the routines that it calls read in and 

process the surface description data the user inputs in 
the SURFACE DATA block. 

GALLING SEQUENCE: CALL SEFCRD 

REFERENCED BY: SEGMENT ROUTINE 

DATARD DATARD 

FILES: NOUT - System output file 

NRIOS - Scratch random I/O file 


ROUTINE NAME: SRFCSI 

DESCRIPTION; This routine reads in the S, R, N, and D cards of the user- 
input SURFACE DATA block, decodes these cards, and writes 
the surface information to file NRIOS for input to the sur- 
face data pass 2 processing. 

CALLING SEQUENCE: CALL SRFCSI 

REFERENCED BY; SEGMENT ROUTINE 

SRFCRD SRFCRD 

FILES; NOUT - System output file 

NRIOS - Temporary random I/O scratch file for passing the surface 

description information to the surface data pass 2 processor 


ROUTINE NAME: START 

DESCRIPTION: This routine calls routines that initialize the label 

commons, x<7rite the TRASYS preprocessor banner page on the 
system output file, read in and process the user- input 
OPTION DATA block, and perform the model collecting and/or 
source editing functions. 

CALLING SEQUENCE: CALL START 

REFERENCED BY: SEGMENT ROUTINE 

TRASYS TRASYS 




ROUTINE NAME : STORDT 


DESCRIPTION: This routine writes all header records to internal data 

files and keeps count of the internal files. 

CALLING SEQUENCE: CALL STORDT (NUNIT, ISTEF, LABEL 1 , LABEL2) 

NUNIT - Unit to contain data 
ISTEP - Step number 
LABELI ” Identifies 1 
LABEL2 - Identifies 2 

REFERENCED BY: SEGMENT ROUTINE 

SDPSS2 AUTO CM 

FRMFRD FINAL 

SHDWRD SHDWRD 

FLUXRD DIBLDR 

GRSPRD CRSPRD 


ROUTINE NAME: STUFDT 

DESCRIPTION: This routine stores all data directly applicable to each 

node/surface tliat requires no conversion determined in 
Surface data pass 2. 

GALLING SEQUENCE: CALL STUFDT 

REFERENCED BY; SEGMENT ROUTINE 

SDPSS2 SDTPS2 


ROUTINE NAME: TIME (CPC system routine) 

DESCRIPTION: This routine will return the current clock time in Hollerith 

code of the format. . . 

HH.MM.SS (CDC) HHMMSS (UNTVAC) 

CALLING SEQUENCE; CALL TIME (IME) 

TME - I’he returned time variable 

REFERENCED BY: SEGMENT 

TRASYS PAGEH 
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ROUTIM NAME : TPERR 


DESCRIPTION: This subroutitie writes out a message whenever a tape/ 

file read error is encountered and aborts run, 

CALLING SEQUENCE; CALL TPERR (ISUB,NUNIT) 

ISUB “ Name of subroutine where read error was 
encountered 

NUNIT - Unit being read when read error was 
encountered 

FILE: NOUT - System output file 


REV 
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ROUTINE NAME: TPGEN 


DESCRIPTION: This routine writes the variable record pointers for the 

random I/O file and the TRASYS processor logical file 
designators to the NSQNTL file. 

CALLING SEQUENCE: CALL" TPGEN 

REFERENCED BY; SEGMENT ROUTINE 

TRASYS TRASYS 

FILES: NSQNTL - The preprocessor and processor conununication file 


ROUTINE NAME : TRAPZ 

DESCRIPTION: This routine converts point input for trapezoids to program- 

compatible surface description parameters and sets up Euler 
angles and a position vector to transform the ICS, BCS , or 
CCS into the SCS of the trapezoid. 

CALLING SEQUENCE: CALL TRAPZ (ISCS) 

ISCS - Flag to indicate if points were input in the 
surface coordinate system 

REFERENCED BY: SEGMENT ROUTINE 

SDPSS2 S DTPS 2 

FILES: NOUT - System output file 


ROUTINE NAME: TRASY S 

DESCRIPTION; This routine is the major segment driver for the TRASYS 

preprocessor. It calls five subsegments; and initializa- 
tion segment, data read segment, logic read setment, 
communication file initialization segment, and the wrapup 
segment . 

CALLING SEQUENCE: None 

REFERENCED BY: None (this is tiie main preprocessor segment.) 
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ROUTINE NAME; TRNGLE 

DESCRIPTION; This is an entry point in subroutine TRAPZ that converts 
point input for triangles to program-compatible surface 
description parameters and sets up Euler angles and a 
position vector to transform the ICS, BCS, or CCS into 
the SCS of the triangle. 

CALLING SEQUENCE: CALL TRNGLE (ISCS) 

ISCS - Flag to indicate if points were input in the 
surface coordinate system 

REFERENCED BY: SEGMENT ROUTINES 

SDPSS2 IMAGES 

POLYGN 

FILES: NOUT - System output file 


ROUTINE NAME: TRS3 

DESCRIPTION: This routine transforms a vector in an SCS to a vector in 

the ICS, BCS, or CCS. 

CALLING SEQUENCE: CALL TRS3 (X, Y. Z, A, B, C, RX, RY , RZ, TRAN) 

KEY VARIABLES; X, Y , Z - Components of vector in ICS, BCS, or CCS 

A, B, C “ Components of vector in SCS 

RX, RY , RZ - Position vector of SCS origin in ICS, BCS, 
or CCS 

TRAN - 3x3 transformation matrix 

REFERENCED BY: SEGMENT ROUTINES 

SDPSS2 BOX 

IMAGES 

SDTPS2 


11-79 


ROUTINE NAKE: WCMfrlN 


DESCRIPTION; This routine writes to file NCMPL a Fortran-labeled 

common statement based on the information in the argu- 
ment list. 

CALLING SEQUENCE: CALL WCMMN (NAME, LNM, N, NT) 

NAME - Common name 

LNM - Array containing a list of names to be 
the variables within the common being 
written 

N - Length of the LNM array 
NT - Type of common to write 

= 0 Single-value variable name 
= 1 Single-dimension (vector) names. The 
LNM under this option is in format... 

Word 1 NAME 

2 Dimension 

3 i^IA^^E 
(etc) 

REFERENCED BY: SEGMENT ROUTINE 

LOGICO WCOMO 

FILES; NCMPL - The file that the generated Fortran cards are written to 


ROUTINE NAME; WCOMO 

DESCRIPTION: This routine writes the TRASYS main-program segment commons 

to the NCMPL file. 

CALLING SEQUENCE: CALL WCOMO 

REFERENCED BY: SEGMENT ROUTINE 

LOGICO WPRGO 

WRITNP 
L0GIC2 L0GIC2 

L0GIC3 LOGIC 3 

FILES; NCMPL - The file that the generated Fortran cards are written to 
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ROUTINE NAME : WCOM 


DESCRIPTION: This routine is a collection of entry points that write 

the Fortran common statements of the processor subsegments 


CALLING SEQUENCE: WCOM is not called, only the entry points are called. 


ENTRY POINTS: GALL 

CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 


CM30 - FFPROG segment 
CM40 “ SFPROG segment 
CM50 - NPPROG segment 
CM60 “ OPPROG segment 
CM70 - DIPROG segment 
CM80 - GBPROG segment 
CM90 “ RKPROG segment 
CMIOO " AQPROG segment 
CMllO - QOPROG segment 
CM120 - RBPROG segment 
CM140 - RCPROG segment 
CM150 - DRPROG segment 


REFERENCES: None 


ROUTINE NA1*IE: WGQTO 

DESCRIPTION: This routine writes to the file NCMPL a Fortran-computed 

GO TO statement based on the information in the argu- 
ment list. 

CALLING SEQUENCE: CALL WGOTO (NSS, NIS, NSN, NAME) 

NSS - Starting statement number 
NIS - Number to increment the statement number 
NSN - Number of statement numbers wanted 
NAME “ Name of the variable to place on the 
computed GO TO statement 

REFERENCED BY: STATEMENT ROUTINE 

LOGICO WRITNP 

LOGIC2 LOGIC2 

FILES: NCMPL - The file that the generated processor Fortran cards 

are written to 

REFERENCES: None 
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ROUTINE NAME : WMAP 


DESCRIPTION; This routine writes the load map directives for the 
TRASYS processor. 

CALLING SEQUENCE: CALL WMAP 

REFERENCES : SEGMENT ROUTINE 

LOGIC2 L0GIC2 

LOGIC3 LOGICS 

FILES: NMAP - File containing processor load map directions 


REV 1 
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ROUTINE NAME: VIPRGQ 


DESCRIPTION: This routine writes the processor main-program segment to 

the NCMPL file. 

CALLING SEQUENCE; CALL WPRGO 

REFERENCED BY: SEGMENT ROUTINE 

LOGICO WRITNP 

LOGIC2 LOGIC2 

FILES: NCMPL - The file that the generated processor Fortran cards 

are written to 


ROUTINE NAME : WPRG2 

DESCRIPTION; This routine name is an entry point into the routine WROG. 

tJhen this entry point name is called, the Fortran processor 
subsegment driver RDPROG is written to the file NCMPL. 

GALLING SEQUENCE: CALL WPRG2 

REFERENCED BY: SEGMENT ROUTINE 

LOGICO WRITNP 

LOGICS LOGIC3 

FILES: NCMPL - The file that the generated processor Fortran cards are 

written to 
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ROUTINE NAME: WPEG3 


DESCRIPTION: This routine name is an entry point into the routine 

l-JPROG. When this entry point name is called, the Fortran 
processor subsegment driver FFPROG is written to the file 
NCMPL. 

CALLING SEQUENCE: CALL WPRG3 

REFERENCED BY: SEGMENT ROUTINE 

L0GIC3 L0GIC3 

FILES: NCMPL - The file that the generated processor Fortran cards are 

written to 


ROUTINE NAME: WPRG4 

DESCRIPTION: This routine name is an entry point into the routine WPROG. 

When this entry point name is called, the Fortran processor 
subsegment driver SFPROG is written to the file N(HPL. 

CALLING SEQUENCE: CALL WPRG4 

REFERENCED BY: SEGMENT ROUTINE 

LOGICS LOGICS 

FILES; NCMPL - The file that the generated processor Fortran cards are 
written to 
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ROUTINE NAME: WPRG5 


DESCRIPTION: This routine name is an entry point into the routine WPROG. 

When this entry point name is called j the Fortran processor 
subsegment driver NPPROG is written to the file NCMPL. 

CALLING SEQUENCE: CALL WPRG5 

REFERENCED BY: ‘ SEGMENT ROUTINE 

LOGICO WRITNP 

L0GLC3 L0GIC3 

FILES; NCMPL “ The file that the generated processor Fortran cards are 
written to 


ROUTINE NAPiE: WPRG6 

DESCRIPTION: This routine name is an entry point into the routine WPROG. 

When the entry point name is called, the Fortran processor 
subsegment driver OPPROG is written to the file NOIPL. 

CALLING SEQUENCE: CALL OTRG6 

REFERENCED BY: SEGMENT ROUTINE 

LOGICS L0G1C3 

FILES; NCMPL - The file that the generated processor Fortran cards are 
written to 


ROUTINE NAME : WPRG7 

DESCRIPTION: This routine name is an entry point into the routine 

WPROG. When this entry point name is called, the Fortran 
processor subsegment driver DIPROG is written to the NCMPL 
file. 

CALLING SEQUENCE: CALL WPRG7 

REFERENCED BY: SEGMENT R OUTINE 

LOGICS L0GIC3 

FILES; NCMPL - The file that the generated processor Fortran cards are 
written to 
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EDUTINE mME: WFEG8 


DESCRIPTION: This routine name is an entry point into the routine WPROG. 

\ihen this entry point name is called, the Fortran processor 
subsegment driver GBPROG is xjritten to the KCMPL file. 

CALLING SEQUENCE: CALL WPRG8 

REFERENCED BY: SEG^^ENT ROUTINE 

LOGIC! LOGIC! 

FILES: NCMPL - Tlie file that the generated processor Fortran cards are 

written to 


ROUTINE NAME; HPRG9 

DESCRIPTION: This routine name is an i^utry point into the routine liTPROG. 

I'Jhen this entry point name is called, the Fortran processor 
subsegment driver RKPROG is written to the NCMPL file. 

CALLING SEQUENCE: CALL WPRG9 

REFERENCED BY: SEGMENT ROUTINE 

LOGIC! LOGIC! 

FILES: NCMPL - The file that the generated processor Fortran cards are 

written to 

ROUTINE NAME : MPRGlO 

DESCRIPTION; This routine name is an entry point into the routine WPROG. 

Wlien this entry point name is called, the Fortran processor 
subsegment driver AQPROG is written to the NCMPL file. 

CALLING SEQUENCE: CALL WPRGIO 

REFERENCED BY: SEGMENT ROUTINE 

LOGIC! LOGIC! 

FILES: NCMPL - The file that the generated processor Fortran cards are 

written to 
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ROUTINE NAME: WPRGll 


DESCRIPTION: This routine name is an entry point into the routine WPROG. 

When this entiy point name is called, the Fortran processor 
subsegment driver QOPROG is written to the NCMPL file. 

CALLING SEQUENCE: CALL WPRGll 

REFERENCED BY: SEGMENT ROUTINE 

LOGICS LOGICS 

FILES: NCMPL - The file that the generated processor Fortran cards are 

written to 


ROUTINE NAME: WPRG12 

DESCRIPTION: This routine name is an entry point into the routine WPROG. 

When this entry point name is called, the Fortran processor 
subsegment driver RBPROG is written to the NCMPL file. 

CALLING SEQUENCE: CALL WPRG12 

REFERENCED BY: SEGMENT R OUTINE 

LOGICS LOGICS 

FILES: NCMPL - The file that the generated processor Fortran cards are 

written to 


ROUTINE NAME : WPBG13 

DESCRIPTION: The routine name is an entry point into the routine WFROG. 

IJhen this entry point name is called, the Fortran processor 
subsegment driver PLPROG is written to the NCMPL file. 

GALLING SEQUENCE: GALL WRG13 

REFERENCED BY: SEGMENT ROUTINE 

LOGICS LOGICS 

FILES; NCMPL *• Tbe tile that the generated processor Fortran cards are 
written to 
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ROUTINE NAME 


NPRG14 


DESCRIPTION: This routine name is an entry point into the routine t^PROG. 

When this entry point name is called, the Fortran processor 
subsegment driver RCPROG is written to the NCMPL file, 

CALLING SEQUENCE: GALL WPRG14 

REFERENCED BY: SEGMENT ROUTINE 

LOGICS LOGICS 

FILES: NCMPL - The file that the generated processor Fortran cards are 

written to 


ROUTINE NAME : WPRG15 

DESCRIPTION: This routine name is an entry point into the routine WPROG. 

When this entry point name is called, the Fortran processor 
subsegment driver DRPROG is written to the NCMPL file. 

CALLING SEQUENCE: CALL WPRG15 

REFERENCED BY: SEGMENT ROUTINE 

LOGICS LOGICS 

FILES: NCMPL - The file that the generated processor Fortran cards are 

written to 


ROUTINE NAME : WPROG 

DESCRIPTION: This routine is a collection of entry points that writes 

the Fortran main routines of the processor subsegments. 

CALLING SEQUENCES: WPROG is not called, only the entry points are called 

ENTRY POINTS: CALL WPRG2 - RDPROG segment 

CALL WPRG3 ~ FFPROG segment 
CALL WPRG4 - SFPROG segment 
CALL WPRG5 - NPPROG segment 
CALL WPRG6 - OPPROG segment 
CALL WPRG7 - DIPROG segment 
CALL WPRG8 - GBPROG segment 
CALL WPRG9 “ RKPROG segment 
CALL WPRGIO " AQPROG segment 
call WPRGll “ QOPROG segment 
CALL WPRG12 - R8PROG segment 
CALL WPRG13 - PLPROG segment 
CALL WPRG14 - RCPROG segment 
CALL WPRG15 - DRPROG segment 
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ROUnUE NAME: WRXTNP 


DESCRIPTION: This routine writes the processor ODPROG subsegment, which 

overrides the user-input LOGIC DATA block whan an error is 
encountered in the user input. This generated Fortran 
operation data block causes the node plotter to plot each 
o£ the BCS-defined data blocks. 

CALLING SEQUENCE: CALL WRITNP (NBCS, NBGSA) 

NBCS - Number of BCS-defined surface systems to be 
plotted 

NBGSA - The array that contains the names of the 
BCS surface systems 

REFERENCED BY: SEGMENT ROUTINE 

LOGICO LOGICO 

FILES: NOUT - System output file 

NCMPL - The file that the generated processor Fortran cards are 
written to 


ROUTINE NAME; WSCOPE (CPC only) 

DESCRIPTION: This routine displays a fo-word message, with the last word 

zeroed, at the computer operator’s console. The only infor- 
mation passed to this routine is program status information. 

CALLING SEQUENCE: CALL WSCOPE (MSG) 


MSG - Message to be written 


REFERENCED BY: 


SEGMENT 

ROUTINE 

TRASYS 

TRASYS 

START 

START 

DATARD 

DATARD 

LOGICO 

LOGICO 

TPGEN 

TPGEN 
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FILE DEFINITIONS - PREPROCESSOR 

FILE NAME : CMERG 

PROGRAM VARIABLE NAME; NCMG 

UNIT REFERENCE (UNIVAC/JSC) ; 12 

PURPOSE; This file is the user's card image merge file* 
SEGMENT REFERENCE; SEDIT (READ) 


FILE NAME; CMPL 

PROGRAM VARIABLE NAME; NCMPL 

UNIT REFERENCE (UNiVAC/JSC) : 20 

PURPOSE; This is a compile file that contains the processor FORTRAN 
routines generated by the preprocessor. 

SEGMENT REFERENCES: LOGICO (WRITE) 

LOGIC2 (WRITE) 

LOGICS (WRITE) 


FILE NAME; DATAI 
PROGRAM VARIABLE NAME: NDI 


UNIT REFERENCE (UNIVAC/JSC) : 4 


PURPOSE: This file is generated in the SEDIT segment from the INPUT, 

EDITI, CMERG, and EMERG files. It contains the users input 
source data, which are read as the input data file in all 
other segments of the preprocessor. 


SEGMENT REFERENCES: SEDIT (WRITE) 

GNEDO (READ) 
QUANRD (READ) 
ARRYRD (READ) 
SKIRD (READ) 
SRFCRD (READ) 
BCSRD (READ) 


DCMNRD (READ) 
FEMFRD (READ) 
SHDWRD (READ) 
FLUXRD (READ) 
CRSPRD (READ) 
LOGICl (READ) 
LOGICS (READ) 
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FILE NAME; DIR 


PROGRAM VARIABLE NAME: NDIR 

UNIT REFERENCE (UNIVAC/JSC) : 22 

PURPOSE; This Is the flux data restart file, which contains the flux 
data request matrix that is generated by the preprocessor 
and communicated to the processor. 

SEGMENT REFERENCES; FLUXRD (WRITE) 

FILE NAME : EMERG 

PROGRAM VARIABLE NAME; NEMG 

UNIT REFERENCE (UNIVAC/JSC) ; 13 

PURPOSE; This file is the user's :dit input merge file. It con- 
sists of a RSI file, previously generated as an RSO file, 
mounted and referred to as EMERG. 

SEGMENT REFERENCES: MCOLL (WRITE) 

SEDIT (READ) 

FILE NAME: FFR 

PROGRAM VARIABLE NAME: NFFR 

UNIT REFERENCE (UNIVAC/JSC) : 21 

PURPOSE: This is the form factor restart file that contains the form 

factor request matrix generated by the preprocessor and 
communicated to the processor. 

SEGMENT REFERENCES; FRMFRD (WRITE) 


FILE NAME; GBIRR 

PROGRAM VARIABLE NAME; NGBIRR 

UNIT REFERENCE (UNIVAC/JSC) : 23 

PURPOSE; This file contains the correspondence data that were processed 
in the correspondence data read segment in the preprocessor. 
These data are communicated to the processor through this file, 

SEGMENT REFERENCES; CRSPRD (WRITE) 
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FILE NAME: INPUT 


PROGRAM VARIABLE NAME: NIN 

UNIT REFERENCE (UNIVAC/JSC) : 5 

PURPOSE: This file is the system input file, usually the card reader. 

SEGMENT REFERENCES: INITAL (READ) 

MCOLL (READ) 

SEDIT (READ) 


FILE NAME; OUTPUT 

PROGRAM VARIABLE NAME: NOUT 

UNIT REFERENCE (UNIVAC/JSC) ; 6 

PURPOSE: This file is the system output file, usually the line printers. 

SEGMENT REFERENCES: All (WRITE) 

FILE NAME ; PLSR 

PROGRAM VARIABLE NAME; NPLSR 

UNIT REFERENCE (UNIVAC/JSC) ^ 25 

PURPOSE: This file contains the restart shadow factors that are gener- 

ated by the preprocessor and passed to the processor. 

SEGMENT REFERENCES: SHDWRD (WRITE) 

FILE NAME: PUNCH 

PROGRAM VARIABLE NAME: NPNCH 

UNIT REFERENCE (UKIVAC/JSC) : 7 

PURPOSE; This file is the system punch file. 

SEGMENT REFERENCES: SEDITI (WRITE) 
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FILE NAME: TAPEl 


PROGRAM VARIABLE NAME: NS Cl 

UNIT reference (UNIVAC/JSC) : 1 

PURPOSE: This file is scratch file 1, which is used as temporary 

storage for intermediate data generated within the preprocessor 
segments. It is also used to communicate data between segments. 

SEGMENT REFERENCES: QUANRD (WRITE & READ) FRMFRD (WRITE & READ) 

A-RRYRD (WRITE & READ) LOGICl (WRITE) 

SKIRD (WRITE & READ) L0GIC2 (READ) 


FILE NAME: TAPE2 

PROGRAM VARIABLE NAME: NSC2 

UNIT REFERENCE (UNIVAC/JSC) : 2 

PURPOSE: T’lis is scratch file 2, which is used as temporary storage 

for intermediate data generated within the preprocessor 
segments. Is is also used to communicate data betii^een segment. 

SEGMENT REFERENCES: QUANRD (WRITE & READ) 

ARRYRD (WRITE & READ) 

FRMFRD (WRITE & READ) 


FILE NAME: TAPE 3 

PROGRAM VARIABLE NAME: NSC3 

UNIT REFERENCE (UNIVAC/ JSC) : 3 

PURPOSE: This is scratch file 3, which is used as temporary storage 

for intermediate data generated within the preprocessor 
segments. It is also used to communicate data between segments. 

SEGMENT REFERENCES: FRMFRD (WRITE & READ) 
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FILE NAME ; RIO 

PROGEIAM VARIABLE NAME: NRIO 

UNIT REFERENCE (UNIVAC/JSC) : 8 

PURPOSE: This is the main working random I/O file used in communi- 

cating data from one segment to another within the prepro- 
cessor, and passing array data, quantities data, and surface 
data to the processor. 

SEGMENT REFERENCES: INITAL (WRITE) SDPSS2 (WRITE) 

QUANRD (WRITE) BCSRD (WRITE) 

ARRYRD (WRITE) LOGIC 2 (WRITE & READ) 


FILE NAME : RI02 

PROGRAM VARIABLE NAME: NRIO 2 

UNIT REFERENCE (UNIVAC/ JSC) : 10 

PURPOSE; This file is used to store data generated by SRFCRD for 
use by SDPSS2. 

SEGMENT REFERENCES: SRFCRD (WRITE) SDTPSS2 (READ) 


FILE NAME; RIO S 

PROGRAM VARIABLE NAME: NRIOS 

UNIT REFERENCE (UNIVAC/ JSC) : 9 

PURPOSE; This file is the scratch random I/O file. It is’used as 
temporary storage for intermediate data generated within 
the preprocessor and is also used to communicate data 
between segments. 

SEGMENT REFERENCES: SRFCRD (WRITE) 

SDPSS2 (READ) 


FILE NAME : RSI 

PROGRAM VARIABLE NAME : NRSI 

UNIT REFERENCE (UNIVAC/JSC) : 14 

PURPOSE: This file is the master restart input file, 

SEGMENT REFERENCES: EDIT 

FILE NAME : R50 

PROGRAM VARIABLE NAME : NRSO 

UNIT REFERENCE (UNIVAC/jSC) : 15 

PURPOSE: This file is the master restart output file. 

SEGMENT REFERENCES; START 

FILE NAME : SQNTL 

PROGRAM VARIABLE NAME: NSQNTL 

UNIT REFERENCE (UNIVAC/ JSC) : 16 

PURPOSE: This is the sequential conmiunication file between the 

preprocessor and processor. 

SEGMENT REFERENCES: TPGEN (WRITE) 

FILE NAME: MAP 

PROGRAM VARIABLE NAME: NMAP 

UNIT REFERENCE (UNIVAC/JSC) ; 27 

PURPOSE: This file contains the processor map directives. 

SEGMENT REFERENCES: L0GIC2, LOGICS 
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D. VARIABLE DEFINITIONS - PREPROCESSOR 


LABELED COMMON /CARD/ 

NCDTYP - Variable that contains the first 6 characters of the 
last input card read 

IND - Array variable that contains the 66 characters from Col. 

7 thru Col. 72 of the last input card read 
NNED - New edit number of the last card read 

NOED " Old edit number of the last card read 

LMOD - Modification label of the last card read 
ID - Number of decoded words in the NDOT and IDOT arrays 

NDOT - The array containing the decoded words of the last input 

card that was decoded 

IDOT - The array containing the type of word decoded in the 
NDOT array 

= -(IRI) Integer word 

= -ClRF) Floating-point word 

= -(IRS) Special character 

= + NUMBER Hollerith word character count 

JDOT - The array containing the start Col- numbers of the input 
card data fields 

IBE - Last word in the INB buffer array that has been processed 

IBT “ Last word in the INB buffer array filled with input cards 

INB - Buffer array that the 46 cards are read into 
lEOF - End of file flag 

= 0 No end of file encountered 
= 1 End of file has been encountered 
NER “ Subroutine independent error flag 

NRECN - The record number currently being processed on file DTI 
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LABELED COMMON /CHRCTR/ 


NIRA 

NlRB 

N2RG 

NlRD 

NlRE 

NlRF 

NlRG 

NlRH 

NlRI 

MlRK. 

NlRL 

NlRM 

NlRN 

NIRO 

RlRP 

NlRQ 

NIRR 

WIRS 

NlRT 

NlRU 

NlRV 

NIRW 

NlRX 

NlRY 

NlRZ 

NIRO 

N1R9 

NlRPR 

NIRDO 

NIRAS 

NlRSL 

NlRPL 

NIRMI 

NIRCM 

WlREQ 

NIRLP 

NlRRP 

NIRBK 
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IRA 
IRB 
IRC 
IRD 
IRE 
IRF 
IRG 
IRH 
IRI 
IRK 
• IRL 
IRM 
IRN 
IRO 
IRP 
IRQ 
IRR 
IRS 
IRT 
IRU 
IRV 
IRW 
IRX 
IRY 
IRZ 
IRO 
1R9 
IR. 
1R$ 
IR* 

w 

1R+ 

IR- 

IR, 

1R= 

IRC 

IR) 

IR 



LABELED COMMON /CONST/ 


NPAGE 

WLINE 

MLINE 

MR77 

MR67 

MR127 

MRBLK. 

NBLANK 

MAXBC 

MAXF 

NA 

ISPCL 

INTEG 

IFLT 

MLABEL 

NONE 

INPTSF 

lEQFF 

FFZERO 

FFNAC 

MR47 


- Page number of the current page being written to the 
output file 

- Line number of the last line written to the output file 

- Maximum number of lines allowed on an output page less 
3 lines 

- Variable used in masking operations 

- Variable used in masking operations 

- Variable used in masking operations 

- Variable used in checking for blank characters 

- A word containing an all-blank Hollerith code 

- Maximum length of blank common 

- Maximum core field length possible for the run 

- Contains the Hollerith code (N/A) 

- Denotes a special character decoded (-1RS) 

- Denotes an integer word decoded (-1RI) 

•' Denotes a floating-point number decoded (-1RF) 

- Current edit modifier label 

- Contains the Hollerith code (none) 

- Multi-shadow factor input block flag 

- Equivalent form factor flag 

- Zero form factor flag 

- Flag to indicate if node array encountered in form factor 
data 

- Variable used in masking operations 
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LABELED COMMON /CPCORE/ 


TCPST 

ICEMC 

TCPE 

ICRE 

TCPQ 

ICRQ 

TCPA 

ICRA 

TCPSl 

ICRSl 

TOPS 2 

ICRS 2 

TCPB 

ICEB 

TCPPF 

ICRFF 

TCPSW 

ICRSW 

TCPF 

ICRF 

TCPC 

ICRC 

TCPO 

ICRO 

TCPSR 

ICRSR 

TCPD 

I CRD 

TCPS 

ICRS 

TCPR 

ICRR 

IPPST 

IPPIN 

IPPE 

IPPQ 

IPPA 

IPPSl 

IPPS2 

IPPB 

IPPFF 

IPPSW 

IPPF 

IPPC 

ippn 

IPPSR 

IPPD 


- Central processor (GP) time of start of preprocessor 
“ Minimum core for model collect editing 

- CP time for source editing 

- Minimum core for source editing 

- CP time for compiling quantities data 

- Minimum core for compiling quantities data 
“ CP time for compiling array data 

- Minimum core for compiling array data 

- CP time for compiling surface data in pass 1 

“ Minimum core for compiling surface data in pass 1 

- CP time for compiling surface data in pass 2 

- Minimum core for compiling surface data in pass 2 

- CP time for compiling BCS data 

- Minimum core for compiling BCS data 

- CP time for compiling form factor data 

- Minimum core for compiling form factor data 

- GP time for compiling shadow data 

- Minimum core for compiling shadow data 

- CP time for compiling flux data 

" Minimum core for compiling flux data 

- CP time for compiling correspondence data 

“ Minimum core for compiling correspondence daua 

- CP time for compiling operation data 

- Minimum core for compiling operation data 

- CP time for initialization - RSTRTO and SQNTL 
" Minimum core initialization - RSTRTO and SQNTL 

- CP .time for processing documentation data block 

- Minimum core for processing documentation data 

- CP time for preprocessing subroutine data 

- Minimum -core for compiling subroutine data 

- CP time for preprocessing restart data 

- Minimum core for processing restart data 

- Preprocessor (PP) time of start of preprocessor 

- PP time for initialization 

“ PP time for source editing 

- PP time for compiling quantities data 

- PP time for compiling array data 

- P? Time for compiling surface data in pass 1 

” PP time fcr compiling surface data in pass 2 

“ PP time for compiling BCS data 

“ PP time for compiling form factor data 

- PP time for compiling shadow data 
“ PP time for compiling flux data 

- PP time for compiling correspondence data 
“ PP time for compiling operations data 

" PP time for initialization of SQNTL and RSTRTO 

- PP time for documentation data block 
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LABELED COMMON /GPCORE/ (continued) 


IPPO - PP time for processing documentation, data block 

IPPS “ PP time for compiling subroutine data 

IPPR - PP time for dbmpillng restart data 

ICPIN - CP time for initializatio.i 

TCRIN - Minimum Core for initialization 
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LABELED COMMON /CVCONS/ 


NBITS - Maximuin number of bits in a computer word 
NBITC - Maximum number of bits in a Ho llerith character 
M^1T72 - Minimum number of cbmputer words tq read in Col. 1 
Col. 72 of an input card 

NW7T72 - Minimum number of computer words to read in Col. 7 
Col. 72 of an input card 

NCHAR - Number of Hollerith characters in a computer word 


LABELED, COMMON /CVFAC/ 

DTR - Degrees to radians 
B.TD - Radians to degrees 
PI - PI 


LABELED COMMON /DIMES/ 


NSRFC " Number of surfaces read in by the preprocessor 
M'lDSMX - lfe.ximum number of wc’ is on longest surface record 
NOFBCS - Number of BCS read in by the preprocessor 
NTSPSF - Total number of specular surfaces 


thru 

thru 
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LABELED COMMON /ERROR/ 

NERR - Total number o£ errors detected by the preprocessor 

NGTW - Total number of cautions detected by the preprocessor 

KERR - Array of line numbers where errors exist 

KGTN - Array of line numbers where cautions exist 

lERR - Number of error words filed in the KERR array 
ICTN - Number of caution words filed in the KGTN array 


LABELED COMMON /FILE/ 


NSCl 

NSC2 

NSC3 

NDI 

NIN 

NOUT 

NRIO 

NRIOS 

NCMG 

NEMG 

NRSI 

NRSO 

NSQNTL 

NRTI 

NCMPL 

NFFR 

NDIR 

NGBIRR 

NPLSR 

NTQR 

NDSERl 

NTJSER2 

NFE 

NDILB 

NGBIR 

NGBSO 

NPLS 

NTQ 

NRTO 

NBGDOU 

NPNCH 

NRIO 2 

NINFO 


- Scratch sequential file 1 

- Scrtach sequential file 2 
Scratch sequential file 3 

- Output of source editor and input to preprocessor 

- Card reader input file 
” Printer output file 

- Main working random I/O file 

- Scratch random I/O file 

- Card input merge file 
~ Edit input merge file 

- Restart file input (permanent) 

- Restart file output (permanent) 

- Sequential communication file 
“ Restart file input (temporary) 

- Generated FORTRAN Program output file 

- Form factor restart file 

- Flux data restart file 

“ Correspondence data file 

- Shadow factor restart file 

- Not used 

- User file 1 

- User file 2 

- Not used 

- Not used 

- Not used 

- Not used 

“ Shadow factor file for processor 
“ Not used 

- Restart file output (temporary) 

- MITAS-SINDA BCD interface tape 

- System punch output file 

- Working random I/O file 

- User information file 
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LABELED COMON /FILET/ 

KRSI - 14 
KRSO = 15 
KRTI = 17 
KRTO = 18 
KCKG = 12 
KEMG = 13 
KBCDOU =29 
KTRAJ = 28 


LABELED COMMON /FNAME/ 

These variables are used to transfer file information from the options 
card read-in routine to the edit and history routines 

JCMG = IH 
JEMG = IH 
JRSI = IH 
JRSO = IH 
JRTI = IH 
JBCDOU = IH 
JTRAJ = IH 
JUSERl = IH 
JUSER2 = IH 


LABELED COMMON /HEAD/ 

LMDTE - Date of last modification 

I.vRSN “ Current TRASYS version and modification number 

NDTE - Job run data 

NTME “Job run start time 

NJOB - Job number 

NTITLE - Array containing the TRASYS internal title that is printed 
at the top of each output page 

ITIILE - Array containing the user-input title for the primary 

model. This title is printed at the top of each printed 
page 


LABELED COMMON /IDSTR/ 

IDSTR - Array containing guiding directives for all information 
on the files created by the input data blocks. 
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LABELED COMMON /OPTION/ 


IRSTRT “ Flag for designating a restart or start run 

IPRNTS - Flag for type of source edit printing 

IPKCHS “ Flag for type of source edit punching 

IGO - Flag for editing and preprocessing the data, or just 

editing the data and terminating 
IPRNTE - Flag for type of edit directive output 
IRLBLS - Flag for relabeling the source 
MDLNM - Primary model name that is to be processed 
MDLNHM - Input model name to be changed to MDLNNM after editing 
IDMPD ~ Flag for printing documentation data block 
lERCNT - Preprocessor error - normal continuation flag 
MAXFL - Maximum field length for processor execution 
NOINP - Not used 

IPINFO - Flag for printing information data file 
ERPLOT - Generate node plot in case of errors 


LABEI.ED COMMON /RNDMIO/ 

MRREC - Maximum number of records that can be written to the 
random I/O file, RIO 

NRREG - Number of records that have been written to the random 
I/O file, RIO 

IRINDX " Not used 

MRREC2 - Maximum number of records that can be written to random 
I/O file, RI02 

NRREC2 - Number of records that have been written to random I/O 
file, RI02 

JRINDX - Not used 
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LABELED COMO^l /MDMRC/ 


H 


NLGQV 

NRUQD 

NLUQD 

HLUQV 

NRAND 

NRAPD 

NLAD 

NLAV 

NRXS 

NLIS 

NLIN 

MBCSD 

NLBCSD 

NLBCSR 

NLSD 

NLBCSN 


“ Number of words In the control constant value array 
record 

“ Random I/O record that the user constant directory 
array record Is written to 

“ Number of words in the user constant directory array 
record 

- Number of words in the user constant value array record 

- Random I/O record that the user array name directory 
array is written to 

“ Random I/O record that the user array position directory 
array is written to 

- Number of words in the user array directory array record 

- Number of words in the user array value array record 

- Random I/O record that the surface index is written to 

- Number of words in the surface index record 

“ Number of words in the node index record 

- Random I/O record that the BCS directory is written to 

- Number of words in the BCS directory record 

" Number of words in the BCS index record 

“ Number of words in the step directory 

- Number of BCS names written to the NRBCSN record 


LABELED COMMON /TRACE/ 


ITRC12 

- 

Trace 

printout 

flag 

for 

segment 

MCOLL 

ITRCED 

“ 

Trace 

printout 

flag 

for 

segments 

i SEDIT & GNEDO 

ITRCDM 

- 

Not used 





ITRC21 

“ 

Trace 

printout 

flag 

for 

segment QUANRD 

ITRC22 

- 

Trace 

printout 

flag 

for 

segpient 

ARRYRD 

ITRC23 

- 

Trace 

printout 

flag 

for 

segment 

SKIRD 

ITRC24 

- 

Trace 

printout 

flag 

for 

segment 

SRFCRD 

ITRC25 

- 

Trace 

printout 

flag 

for 

segment 

SDPSS2 

ITRC26 

- 

Trace 

printout 

flag 

for 

segment 

BCSRD 

1111027 

- 

Not used 





ITRCFF 

- 

Trace 

printout 

flag 

for 

segment 

FRMFRD 

ITRC29 

- 

Trace 

printout 

flag 

for 

segments 

FLUXED and CRSPRD 

ITRC31 

- 

Trace 

printout 

flag 

for 

segment 

LOGICl 

ITRC32 

- 

Trace 

printout 

flag 

for 

segment 

L0GIC2 

ITRC33 

- 

Trace 

printout 

flag 

for 

segment 

LOGICS 

ITRCTG 

- 

Trace 

printout 

flag 

for 

segment 

TP GEN 

ITRCRS 

- 

Trace 

printout 

flag 

for 

segment 

RESTRT 

ITRCON 


Trace 

flag to be checked for 


ITRCSF 

- 

Trace 

printout flag 

for 

segment 

SHDWRD 
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EMTRY POINTS - PREPROCESSOR 


EMTRY FOIMT 

GHECTP 

IHRDB 

INRDDC 

INRDO 

INRDOD 

KPIE 

WAQCOM 

WCMGOM 

WDICOM 

WDRCOM 

WFFGOM 

WGBCOM 

WMFGOM 

WNPCOM 

WOPCOM 

WPLCOM 

WQOGOM 

WRBCOM 

WRCCOM 

WSFCOM 

WAQSEG 

WCMSEG 

WDISEG 

WDRSEG 

WFFSEG 

WGBSEG 

^#1FSEG 

MPSEG 

WOCSEG 

WOPSEG 

WPLSEG 

WQOSEG 

l^RBSEG 

l-TRCSEG 

WRDSEG 

WSFSEG 


SITBRODTIHE 

GHEC 

INRDD 

IHRDD 

INRDD 

INRDD 

KPl 

WCOM 

WGOM 

WCOM 

WCOM 

WCOM 

WGOM 

WCOM 

WGOM 

WCOM 

WCOM 

WCOM 

WCOM 

WCOM 

WGOM 

WPROG 

WPROG 

WPROG 

WPROG 

OTROG 

WPROG 

^^r?ROG 

WPROG 

WPROG 

WPROG 

WPROG 

WPROG 

WPROG 

WPROG 

OTROG 

WPROG 


III. PART 2 - PROCESSOR LIBRARY 


A. SEGMENT DEFINITIONS 
SEGMENT NAME : AQCAL 

PURPOSE: This segment computes absorbed heat rates in two wavebands, 

accounting for specular and diffuse reflection (Fig. III-l), 

RESTRICTIONS: Appropriate direct -irradiation and gray-body factors 

must be located on the DI, GBIR, and GBSO files 
before execution- 

CALLING SEQUENCE; L AQCAL 

OUTPUT; Absorbed total heat rates written to file TQ. 


SEGMENT NAME: CMCAL 

PURPOSE: This segment combines form factors according to correspon 

pondence data, resulting in a form factor matrix of 
reduced size. 

RESTRICTIONS: Requires form factors on file FF. 

CALLING SEQUENCE: L CMCAL 

OUTPUT: Writes combined form factors to file FF. 

SEGMENT NAME: DICAL 

PURPOSE: This segment computes solar, planetary, and albedo 

irradiation incident on spacecraft nodes (Fig. III-2) . 

RESTRICTIONS: Execution is possible after previous calls have been 

made to define spacecraft geometry, location in space, 
characteristics and distances of heat source bodies, 
and computation accuracy parameters. 

CALLING SEQUENCE; L DICAL 

OUTPUT: Printed, punched, written to files DI, RTO, and PLS. 


III-l 



Figure IIJ-2 Segment AQCAL Ftou Diagram 
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DM^AIN 


CALL DREDRQ 


CALL DRIMAG 


KR = 1 
IN = 1 




Figure III~3 Segment DECAL Flow Diagram 
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SEG^E1^T NAME* FFCAL 


PURPOSE: This segment calculates all form factors for the active 

configuration (Fig. III-4) . 

EESTRIGTIONS None 

CALLING SEQUENCE; L FFCAL 

OUTPUT; printed, punched, written to file FF. 

SEGMENT NAME; DRCAL 

PURPOSE; This segment computes the components of direct irradiation 
that reach the nodes via specular reflections in other 
nodes. 

RESTRICTIONS: Requires direct irradiation data on file DI. 

CALLING SEQUENCE: L DRCAL 

OUTPUT: Writes direct irradiation that includes specular reflection 

components to file DI. Print file output. 


SEGMENT NAME : RBCAL 

PURPOSE; This segment computes the components of internode form 
factors that result from images of nodes in specular 
surfaces. 

RESTRICTIONS: Requires direct form factor data on file FF. 

CALLING SEQUENCE: L RBCAL 

OUTPUT: Writes modified form factors, containing specular compon- 

ents, to the FF file. 



Figure III-4 Segment FFCAL Flow 
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Figure IU~S Segment RBCAL Flow Diagram 
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SEGMENT NAME: GBGAL 


PURPOSE: This seginent computes and stores the gray-body factor 

matrix (Eig. TlI-6) . 

RESTRICTIONS: Requires form factor data on file FE. 

CALLING SEQUENCE: L GBGAL 

OUTPUT: Writes gray-body factor matrices to GBIR and/or GBSO. 


SEGMENT NAME: NPLOT 

PURPOSE: This segment generates pictorial plots of nodal sur- 

faces (Fig. III-7) - 

RESTRI CHONS : None 

CALLING SEQUENCE: L NPLOT 

OUTPUT: Printed, plot file. 


SEGMENT NAME: OPLOT 

PURPOSE: This segment generates pictorial plots of the space- 

craft in orbit (Fig. III-8) . 

RESTRICTIONS : None 

CALLING SEQUENCE; L OPLOT 

OUTPUT: Printed, plot file. 


SEGMENT NAME: PLOT 

PURPOSE; This segment generates function vs time plots of ab- 
sorbed and incident heat rates and fluxes. When used 
in conjunction with operations block FORTRAN that writes 
data to a plot data unit, this segment provides a gen- 
eral X vs y plot capability (Fig. HI-9) . 

RESTRICTIONS: Reference Subroutine PLUATA 

CALLING SEQUENCE: L PLOT 


OUTPUT: Printed, plot file. 
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Figure III-3 Segment OFLOT Flow Diagram 
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Figure III-9 Segment PLOT FtoiD Diagram 
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SEGMENT NAME : QOCAL 


PURPOSE: This segment accesses absorbed flux data and generates 

orbital average and absorbed flux vs time arrays (Fig. III-IO) . 

RESTRICTIONS: Requires data on the TQ file. 

CALLING SEQUENCE: L QOCAL 

OUTPUT; Printed, punched, BCD tape (File BCDOU) , 


SEGMENT NAME; RCCAL 

PURPOSE: This segment computes radiation conductors and simpli- 

fies and condenses these conductors using the ERN and 
MESS techniques (Fig, III-ll) , 

RESTRICTIONS: Requires data on the GBIR file. 

CALLING SEQUENCE: L RCCAL 

OUTPUT; Printed, punched, BCD tape (file' BCDOU) . 


SEGMENT NAME : RKCAL 

PURPOSE; This segment computes radiation conductor values (Fig. III-12) . 
RESTRICTIONS: Requires data on the GBIR file. 

CALLING SEQUENCE: L RKCAL 

OUTPUT: Printed, punched, BCD tape (file BCDOU). 


SEGMENT NAME : SFCAL 

PURPOSE: This segment computes tables of internode blockage 

(shadow) factors for use in direct irradiation calcu- 
lations. When a complete shadow factor tape is supplied 
SFCAL is executed to pass the shadow tables to file PLSR 
and initialize DICAL to compute irradiations using the 
shadow tables (Fig. III-13) . 

RESTRICTIONS : None 

CALLING SEQUENCE; L SFCAL 

OUTPUT: Printed, binary tape (files NPLS and RSO). 
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Figure III-IO Segment QOCAL Flob) Diagram 
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Figure III-12 Segment RKCAL Flow Diagram 
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B. SUBROUTINE AND FUNCTION DESCRIPTIONS - PROCESSOR LIBRARY 
ROUTINE NAME : ABSBEA 

DESCRIPTION: This routine locates plotting beam at the absolute 

raster coordinate (K, Y) based on 0.0 X, Y 1.0. 

CALLING SEQUENCE; CALL ABSBEA (X, Y) 

REFERENCES; NPMAIN, NPINFO, OPMAIN, OPINFO, PLGRID 


ROUTINE NAME: ADSURE 

DESCRIPTION: This routine provides the user with the capability 

to add an adiabatic "closure” surface to the problem 
geometry. 

CALLING SEQUENCE; CALL ADSURF (NBCSN»NGONFG, ARAS) 

NBGSN - Name of block coordinate system under 
which adiabatic surface is defined 
NCONFG - Configuration name to locate appro- 
priate form factor matrix 
ARAS - Area of adiabatic surface 

REFERENCE: ODPROG (Users operations data) 


ROUTINE NAME : ADD 

DESCRIPTION; This is an entry point in subroutine BUILDC that adds 
additional nodes/surfaces to the existing configura- 
tion (see BUILDC). 

CALLING SEQUENCE; CALL ADD (BCSN) 

BCSN - Left- justified, blank-filled block 
coordinate system name 

REFERENCES: User call in the Operations Data Block 

FILES : NRAN 
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R0UTI1«E NAME: 


AQGMPR 


DESCRIPTION: This routine combines absorbed solar, albedo and 

planetary heat rate data according to correspondence 
data, 

CALLING SEQUENCE: CALL AQGMPR 

REFERENCE: AQMAIN 

FILES : NS CRl , NGBIRR , NOUT 


ROUTINE NAME; AQDATA 

DESCRIPTION: This is a user-called routine to define control constants. 

CALLING SEQUENCE: CALL AQDATA (IGBI, IGBS, RS , RA, RP) 

IGBI “ Step number containing gray-body IR 

IGBS - Step number containing gray-body SOMR 

RS - Solar multiplier 

RA - Albedo multiplier 

RP - Planetary multiplier 

REFERENCES; ODPROG (user’s operations data) 


ROUTINE NAME ; AQEND 

DESCRIPTION: Routine for user logic before termination of AQCAL link. 

REFERENCES: AQMAIN 


ROUTINE NAME; AQHEAD 

DESCRIPTION; This routine writes to output file the default param- 
eters and flags set for AQCAL execution. 

CALLING SEQUENCE; CALL AQHEAD 

REFERENCE : AQMAIN 
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ROUTINE NAME: AQMAIN 


DESCRIPTION: This routine computes absorbed Q data from direct incident 

data and solar/IR gray bodies and stores the data on file 
NTQ. 


CALLING SEQUENCE: CALL AQMAIN 


KEY VARIABLES; 


QAS “ Absorbed solar array 
QAR - Absorbed albedo array 
QAP - Absorbed planetary array 
GBIR ^ Gray-body IR array 
GBSO - Gray-body solar array 
QDS - Direct solar 
QDR - Direct albedo 
QDP - Direct planetary 


REFERENCES; AQPROG 


FILES: Reads NGBIR and NGBSO, NDI writes NTQ 

NOUT, NRSO 


ROUTINE NAME; AQNCHK j 

DESCRIPTION; This routine checks the node array on the data file with 
the active node array to determine if they contain the 
same configuration. If the configurations are different, 
ACNCHK outputs both arrays and aborts, if both configura- 
tions are the same, it returns to the main program, 

CALLING SEQUENCE: CALL AQNCHK (ITEMP, Nl) 

ITEMP - Array of node numbers from data file 
Nl “ Number of nodes 


KEY VARIABLES : NODE - Node number array of actual configuration 

NNOD - Number of active nodes in configuration 

REFERENCES : AQMAIN 


ROUTINE NAME: AQPRE 

DESCRIPTION; This Is a user-called routine in the AQCAL link prior 
to any computations. 

GALLING SEQUENCE; CALL AQPRE 

REFERENCES ; AQMAIN 
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ROUTINE NAME: 


AQERT 


DESCRIPTION: This routine prints out the solar, albedo and planetary 

and total absorbed heat rates. 

CALLING SEQUENCE: CALL AQPRT (IN) 

IN - Node sequence number 

REFERENCE : AQCMPR 

FILE : NOUT 


ROUTINE NAME; AQPROG 

DESCRIPTION; This routine is generated by the preprocessor to call 
into AQMAIN and the remainder o£ the link. 


CALLING SEQUENCE; CALL AQPROG 


KEY VARIABLES : AQQDS 

AQQDR 
AQQDP 
QAS 
QAR 
QAR 
GBSO 
GBIR 
AQTEMP 


- Long direct solar in array NNOD 

- Long direct albedo in array NNOD 

- Long direct planetary in array NNOD 

- Long absorbed solar in array NNOD 

- Long absorbed albedo in array NNOD 

- Long absorbed planetary in array NNOD 

- Long gray -body solar in array NNOD 

- Long gray-body IR in array NNOD 

“ Working storage area in array NNOD 


REFERENCES; TRASYS (root segment) 


ROUTINE NAME: BANNLB 

DESCRIPTION; This routine writes the TRASYS processor banner on 
the output file. 


CALLING SEQUENCE: CALL BANNLB 

REFERENCES ; RDMAIN 
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ROUTINE NAME: BtDPLT 


DESCRIPTION: This routine reads plot data from scratch disk 1 (WSGRl) 

and writes plot data in the proper format on output and 
plot disk IPLUNT. It also scales the X and Y values. 


CALLING SEQUENCE: GALL BLDPLT (DATA, NV, NTIME, TIME) 


DATA 

NV 

NTIME 

TIME. (1, x) 
TIME (2, x) 


- Working storage array 

- Number of lines per grid 

- Number of data points per line 

- Independent variable (time) 

“ Associated step number 


REFERENCES : PLLOAD 


FILES: NSCRl is defined in PLLOAD, IPLUNT is written in plot format. 

NOUT, NSCR2 


ROUTINE NAME : BUILDC 


DESCRIPTION: This routine activates and initializes the configuration 

1 block coordinate system name per call. It also defines 
arrays NODE, AREA, EMISS, TRIR, TESO, SRIR, and SRSO, and 
Initializes configuration counters NN, NS, NSPEC, and NONLY. 

CALLING SEQUENCE: CALL BUILDC (BCSN)/CALL ADD (BCSN) 

BCSN - Left-justified, blank-filled block coordinate 
system name 


KEY VARIABLES* INDXN - Array of .andom access record pointers to active 

nodes defined by BDILDC/ADD 

INDXS - Array of random access record pointers to active 
surfaces defined by BUILDC/ ADD 
NN - Number of nodes plus shadower-only nodes 
NS - Number of surfaces 

NSPEC - Number of specular surfaces 

NNOD - Number of nodes 

NSURF - Number of surfaces 

NONLY - Number of shadower only surfaces 


REFERENCES: User-called in the operations data block only 

PILES ; NRAN, NOUT 
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ROUTINE NAME: 


CHGBLK 


DESCRIPTION: 


This routine reads block coordinate system dat^^transla- 
tions and rotations) from random access file NRAN, redefines 
them, rewrites them to the NRAN file, and flags the BCS 
directory that the BCS has been changed. 


CALLING SEQUENCE; 


CALL CHGBLK (NAME, X, 
ROTY, ROTZ) 


NAME 

-BCS name 

NROTX - 

X 

“ Translation X 

NROTY - 

Y 

- Translation Y 

NROTZ - 

Z 

- Translation Z 

ROTX - 



ROTY - 



ROTZ - 


Y, Z, NROTX, NROTY, NROTZ, ROTX, 


Rotation order of ROTX 
Rotation order of ROTY 
Rotation order of ROTZ 
Rotation angle about X-axis 
Rotation angle about Y-axis 
Rotation angle about Z-axis 


KEY VARIABLES: NBLKDR (1, I) 

(2, I) 
(3, I) 
(4, I) 


BCS name 

Length of block 

Random record number 

Flag transform has been applied 


REFERENCES: User-called in the Operations Data block 


FILES: NRAN - Rewritten 

NOUT 


ROUTINE NAME ; CLEAR 

DESCRIPTION: This routine frees, clears and returns the data written 
on a unit so it may be reused. (Generally used with 
scratch files.) 

CALLING SEQUENCE; CALL CLEAR (lUNlT) 

lUNIT = Unit number to be freed 
REFERENCED BY; All processor segments 


/ 

ROUTI^JE NAME: CMGMBN 

DESCRIPTION: This routine combines form factors, area and property- 

arrays . 

CALLING SEQUENCE: CALL GMCMBN (NODEG ,AREAC ,EMISSC ,ALPHG,TRIRC , 

TRSDC ,SRIRG ,SRSOC ,NLEN) 

NOD EC - Combined node array 

AREAG - Combined area array 

EMISSC - Combined emissivity array 

ALPHC - Combined solar absorptivity array 

TRIRG - Combined infrared transmissivity array 
TRSOC - Combined solar transmissivity array 
SRIRC - Combined infrared specular reflectivity 
array 

, SRSOC - Combined solar specular reflectivity 
array 

NLEN - Maximum length of form factor record 
on KFF 

REFERENCE: CMMAIN 

FILES: NSCRl, CR2, NSCR3,NFF, NOUT, NRSI, NRSO 


ROUTINE NAME: GMDATA 

DESCRIPTION; This routine is called by the user to set up control 
constants for the CMCAL program segment, 

CALLING SEQUENCE: CALL CMDATA (NFIG,NFCO,NTYP, lAUTO, FFCM,FFPRT) 

KEY VARIABLES: Reference Users Manual, Appendix D 

REFERENCE: ODPROG (User's Operatinns Data) 


ROUTIN?i NAME: CMESUM 

DESCRIPTION: This routine writes combined form factor sums to the 

output file, 

CALLING SEQUENCE: CALL CMESUM (NNOD, NODE) 

NNOD - Number of nodes after combining 
NODE - Combined node array 

REFERENCE: CMCMBN, CMRSI 

FILES : NOUT 
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ROUTINE NAME; 


CMINIT 


DESCRIPTION: This routine sets up temporary node, area and property 

arrays and checks for duplicate node numbers, 

CALLING SEQUENCE: CALL CMINIT (NODEC ,EMMISSC,ALPHC ,TRIRC,TRSOC , 

SRIRC,SE.SOC) 

KEY VARIABLES; Reference subroutine CMGMBN 

REFERENCE: CMMAIN 

FILES; NGBIRR, NOUT, NSGR3, NRSI 


ROUTINE NAME: CMMAIN 

DESCRIPTION: Driver routine for form factor combining segment, 

CALLING SEQUENCE: CALL CMMAIN 

REFERENCE: CMPROG 

FILES: NOUT, NFF, NRSI 

ROUTINE NAME: CMOUT 

DESCRIPTION; This routine writes individual combined form factors 
to the output file. 

CALLING SEQUENCE: CALL CMOUT (IN, JN, NODE, AREA, NNCM) 

IN - Row number 
JN - Column number 
NODE “ Node array 
AREA - Area array 

NNCM - Number of nodes after combining 
REFERENCE: CMCMBN, CMRSI 

FILE: NOUT 
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ROUTINE NAME: 


CMPROG 


DESCRIPTION: This is the preprocess or --generated routine that calls 

in the form factor combining segment. 

REFERENCE: TRASYS ^ ot segment) 


ROUTINE NAME; CMRSI 

DESCRIPTION; This routine reads combined form factors from the 
master restart tape. 

CALLING SEQUENCE: CALL CMRSI (ICOMP, NOD EC, AREAC, EfiXSSC, ALPHC, 

TRIRC, TRSOC, SR IRC , SRSOC) 

KEY VARIABLES: Combined node, area and property arrays (reference 

subroutine CMGMBN) 

REFERENCE; GMMAIN 

FILES: NRSI, NOUT, NFF, NRSO, NSCRl 


ROUTINE NAME; CMRSUM 

DESCRIPTION: This routine ^■nrites combined form factor sums to 

the output file* 

CALLING SEQUENCE: (NODE, AREA, SUM, IN) 

NODE - Combined node array 

AREA - Combined area array 

SUM - Array of combined form factor sums 

IN - Node sequence numbers 

REFERENC E : CMGMBN, CMRS I 


FILE: NOUT 


ROUTINE NAME; 


CMSAVE 


DESCRIPTION; This routine writes the combined node, area and property 
arrays to a random access file. 

GALLING SEQUENCE: CALL CMSAVE (ISAVE) 

ISAVE - Random file record length flag 

REFERENCE: CMCMBN, CMRSI 

FILES : NRARR , NOUT 


ROUTINE NAME: CMTRAN 

DESCRIPTION: This routine transposes the combined form factor 

matr?lx. 

CALLING SEQUENCE: CALL CMTRAN (DATA(l), DATA(Nl) , DATA(l), DATA(Nl) , 

:.. TA(N2), DATA(N3), DATA(N2) , DATA(N3), NB, NND, 
NR, mECT) 


DATA 

“ Array 

containing form factors 

NB 

" Number 

of blocks 

NND 

- Number 

of nodes 

NR 

- Number 

of records per block 

NRECT 

- Number 

of records, total 

REFERENCE: CMMAIN 



FILES; NFF, NSCRl, NSCR2 




ROUTINE NAME: CROSS 


DESCRIPTION: This routine forms a vector cross-product of the form. 


CROS = VECl X VEC2 

CALLING SEQUENCE: CALL CROSS (VECl, VEC2, CROS) 

VECl (3) - X, Y, Z in array 

VEC2 (3) - X, Y, Z in array 

CROS (3) “ X, Y, Z in array 

REFERENCES ; DIL0C2 


ROUTINE NAME: PI CALF 

DESCRIPTION: This routine computes the direct incident albedo radia- 

tion QDR (IN) and direct planetary radiation QDP (IN), forms 
an optimum elemental grid based on altitude and node posi- 
tions/shadowing, and saves elemental planetary factors on 
disk and applies shadowing. It is the main control routine 
for computing albedo and planetary radiation. 


calling SEQUENCE; CALL DICALP (IRSPH) 


KEY VARIABLES: SHADR 

SHADP 
SUMR 
SUMP 

QDR (IN) 
QDP (IN) 
IS FT 
PLTYPE 


Shadow factor for albedo 
Shadow factor for planetary 
Elemental albedo factors 
Elemental planetary factors 
Complete incident albedo for Node IN 
Complete incident planetary for Node IN 
Flag to indicate shadow factor tape 
Flag to indicate to save or read planet 
elemental factors 


REFERENCES : DITYPE 

PILES: NPLS, NOUT, NRSO 
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ROUTINE NAME; 


DIGALS 


DESCRIPTION: This is the main controlling routine to compute direct solar 

incident fluxes. It forms an optimum elemental grid based 
on position and shadowing. 

GALLING SEQUENCE: CALL DICALS 

KEY VARIABLES: ISFT - Flag to indicate shadoxif factor tape 

SHADS “ Shadow factor 

QDS (IN) - Incident solar flux for Node IN 
REFERENCES : DIMAIN 

FILES : None 


ROUTINE NAME; DICCV 

DESCRIPTION; This routine defines the position vector, given clock and 
cone angles and the vector magnitude. 

CALLING SEQUENCE; CALL DICCV (VMAG, CL, CO, PVEC) 

VMAG - Magnitude 

CL - Clock angle in degrees 

CO - Cone angle in degrees 

PVEC(3) “ Position vector 

REFERENCES; DIL0C2 


ROUTINE NAME; DIGLCO 

DESCRIPTION: This routine computes clock and cone angles given a 

position vector, 

CALLING SEQUENCE: CALL DIGLCO (X, Y, Z, CL, GO) 

X, Y, Z “ Components of position vector 
CL, CO - Clock and cone angles 

REFERENCE: DICALP, DICALS, DILOC, DIPRNT 
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ROUTINE NAME : DICOMB 


DESCRIPTION; This routine performs a matrix multiplication of 
the form 

A'B = C 

where 

Aj B, and C are 3x3 matrices. 

CALLING SEQUENCE: GALL DICOMB (A, B, G) 

A - 3 X 3 matrix 
B - 3 X 3 matrix 
C - Resultant 3x3 matrix 

REFERENCES: DILOC, DIL0C2 


ROUTINE NAME; DICOMP 

DESCRIPTION; This is a user routine to define compute flag or step 
numbers , 

CALLING SEQUENCE: CALL DIGOMP (IS, lA, IP) 

IS - Solar flag to compute, zero, and retrieve 
data 

lA - Albedo flag to compute, zero, and retrieve 
data 

IP “ Planetary flag to compute, zero, and 
retrieve data 

REFERENCES: User-called routines in the Operations Data block. 
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ROUTINE NAME : DIDCS 


I 


DESCRIPTION: This routine computes direction cosines, given the Euler 

angles and the order to perform the rotations. 

CALLING SEQUENCE: CALL DIDCS (II, JJ, KK, PHI, PSI,- OMI, TRAN) 


II ) - Order of rotations 

Jj / for angles PHI, PSI and OMI, respectively 

KK ) 


PHI - Rotation angle about Z-axis 

PSI - Rotation angle about Y-axis 

OMI - Rotation angle about X-axis 

TRAN (3,3) - Resultant direction-cosine matrix 


REFERENCES ; None 



ROUTINE NAME : PI DENT 

DESCRIPTION: This routine generates a 3 x 3 identity matrix in A. 

CALLING SEQUENCE: CALL DIDENT (A) 

A - 3 X 3 array withi = j = l, 

REFERENCES: DILOC , DIDT2, SURFP 
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ROUTINE NAME; DIDTl 


DESCRIPTION: This is a user-called routine to define the DI link 

parameters . 

CALLING SEQtfENCE; See users manual, Appendix D 
REFERENCES: User-called routine in the Operations Data block 

ROUTINE NAME: DIDTIS 

DESCRIPTION; This is a user-called routine to define the DI link 
parameters . 

CALLING SEQUENCE; See users manual. Appendix D 
REFERENCES; User called routine in the Operations Data block 

ROUTINE NAME; DIDT2 

DESCRIPTION; This is a user-called routine to define the DI link 
parameters . 

CALLING SEQUENCE: See users manual, Appendix D 

REFERENCES: Operations Data, user-called 

ROUTINE NA^E: DIDT2S 

DESCRIPTION; Short form call to define DI link parameters. 

CALLING SEQUENCE; See users manual. Appendix D 
REFERENCES: User-called routine in the Operations Data block 

ROUTINE NAMES: DIDT3 , D IDT3S 

DESCRIPTION. User called routines to define DICAL segment parameters. 
CALLING SEQUENCES: See Users Manual, Appnedix D. 

REFERENCES: ODPROG (User’s Operations Data) 
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DIELEM 


ROUTINE NAME; 


DESCRIPTION; Computes position and area vectors for all elements on a 
surface given the number of elements in each direction on 
a node, 

CALLING SEQUENCE: CALL DIELEM (ILP, DATA, TRAN, RX, RY, RZ, NTOT) 

ILP 

DATAC5) 

TRAN(3,3) 

RX, RY, RZ 
NTOT 


KEY VARIABLES: NEST 

SFPV (400, 3) 
SFAV (400, 3) 

REFERENCES: DICALS, DICALP 


ROUTINE NAME : DIELSL 

DESCRIPTION: This routine computes the number of elements in tr<io direc- 

tions on a node in an attempt to provide square elements, 
given the total number required. 

CALLING SEQUENCE; CALL DIELSL (NB, NG, ILP, DATA, NTOT) 

NB - Computed number of elements in beta direction 

NG - Computed number of elements in gamma direction 

ILP - Surface type 

DATA(5) - Nodal dimensions 

NTOT - Total number of elements required 

REFERENCES: DIELEM 


ROUTINE NAME: DIEND 

DESCRIPTION: This routine outputs final arrays of data (QDS - incident 

solar, the QDR - incident albedo, and QDP - incident 
planetary) to the NDI disk file for later reference. It 
end-files NRTO if a restart tape is desired and end-files 
NPLS if a PLTYPE (save) flag has been set. 

CALLING SEQUENCE: CALL DIEND 

REFERENCES : DIMAIN 

FILES: Writes the NDI file and closes NPLS 
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- Type of geometric node 

- Dimensions of node 

- Direction cosines of node 

- Translation vector 

- Number of elements node is to be divided 
into 

- Number of elements (maximum = 400) 

- Node elemental position vectox* 

- Node elemental normal vector 



ROUTINE NAME: DIENDP 


DESCRIPTION; This is a user-callable routine to intervene after completing 
alfaedo/planetary computations. 

CALLING SEQUENCE; CALL DIENDP 

REFERENCES : DIMAIN 


ROUTINE NAME: DIENDS 

DESCRIPTION; This is a user-callable routine to intervene after completing 
solar computations. 

CALLING SEQUENCE; CALL DIENDS 

REFERENCES ; DIMAIN 


ROUTINE NAME: DIGTST 

DESCRIPTION; This integer function routine determines if GN< GT< 
GX. If true, the function is set equal to 0; if false 
the value is set equal to 1. 

CALLING SEQUENCE: I = DIGTST (GN, GX, GT) 

GN - Lower bound 
GX - Upper bound 
GT - Intermediate value 
I - 0, 1 flag 

REFERENCES: DISHAD 
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ROUTINE NAME: 


BILOC 


DESCRIPTION: This routine computes necessary orbital parameters from 

input values for PERIOD, TRUEAN, TIMEPR, SUNPV, and various 
orbital transformations, 

CALLING SEQUENCE: CALL DILOC 

KEY VARIABLES: PERIOD - Orbit period 

TRUEAN - True anomaly 
TIMEPR - Present time 
SUNPV *- Stm posit5on vector 

PLDC - Matrix of direction cosines to transform vectors 

in the planet-oriented VCS to the user-defined VCS 

RJiFERENCES : DIMAIN 

FILES : NOUT - Sys tem output file 


ROUTINE NAME; DIL0C2 


DESCRIPTION: This routine computes necessary orbital parameters, given 

the clock and cone angles to the sun and planet. 


CALLING SEQUENCE; CALL DIL0C2 


KEY VARIABLES: 


SUNCI - Sun clock angle 
SUNCO - Sun cone angle 
PL CL - Planet clock angle 
PLCO - Planet cone angle 
lORBIT - Flag for type of orbit 


REFERENCES; DILOC 


FILES : NOUT 

NSCR3 


System output file 
Scratch file 
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ROUTINE NA>IE: 


•DIMAIN 


DESCRIPTION: This routine is the main driving logic for computing direct 

incident fluxes and in^ ’ udes the main computation loops* 

CALLING SEQUENCE; CALL DIMAIN 

KEY VARIABLES; IN — Sequence to node number currently being computed 
REFERENCES: DIPROG (preprocessor-generated) 


ROUTINE NAME; DIOUTP 

DESCRIPTION: This routine allows the user to change the type or form of 

printed or piinched data. It normally calls: 

DIPRTP “ To print the alhedo/planetary fluxes 

DIPNHP - Tp punch the albedo/planetary fluxes 

DITPP - To write the albedo /planet ary fluxes to the RTO 

tape in restart format 

CALLING SEQUENCE: CALL DIOUTP 

REFERENCES ; DIMAIN 


ROUTINE NAME : PI OUTS 

DESCRIPTION; This routine allows the user to change the type or form of 
printed/punched data. It normally calls; 

DIPRTS - To print the solar fluxes 

DIPNHS - To punch the solar fluxes 

DITPS - To write the solar fluxes to the RTO tape in the 

restart format 

CALLING SEQUENCE; CALL DIOUTS 
REFERENCES : DIMAIN 
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ROUTINE NAME: 


DIPLNS 


DESCRIPTION: This routine computes planet position and area vectors 

based on the orbit and accuracy parameters, and determines 
the emissive power of the planet element and view factor 
from the element to the sun. 

CALLING SEQUENCE: CALL DIPLNS (NPEL) 

NPEL - Total number of desired elements on planet 

KEY VARIABLES; PLPVT (400, 3) - Array of planet position vectors 

PLAVT (400, 3) - Array of planet-normal vectors (magnitude 
= area) 

ALBP (400) - Array of planet-to-sun view factors 

PLNF (400) ~ Array of planet- element emissive powers 

REFERENCES: DICALP 


ROUTINE NAME: DIPNHP 

DESCRIPTION; This routine punches albedo and planetary fluxes in a format 
acceptable for restart if the DIPNCS flag has been set to 
3HPUN. 

CALLING SEQUENCE: CALL DIPNHP 

KEY VARIABLES: DIPNHP - Flag to determine if punched cards are requested 

NODE - Array of node numbers 
QDR - Array of albedo values 
QDP - Array of planetary values 

REFERENCES: DIOUTP, DIPTPR 

FILES ; NPUN 


ROUTINE NAME; DIPNHS 

DESCRIPTION: This routine punches solar fluxes in a restart format com- 

plete with proper header cards generated on the initial 
call to the routine. 

CALLING SEQUENCE; CALL DIPNHS 

KEY VARIABLES: 1ST - Flag to determine if this is the initial call 

REFERENCES : DIOUTS, DIPTSR 

FILES : NPUN 
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I 
\ 
\ 

DIPEEP ! 

j 

This routine can be replaced by the user prior to computing 
albedo/planetary fluxes, . ? 

' ; 

CALLING SEQUENCE: CALL DIPREP | 

I 

REESRENCES : . DIMAIN ! 


i 

ROUTINE NAME; DIPRES i 

DESCRIPTION: This routine can be replaced by the user prior to computing ' 

solar fluxes. 

CALLING SEQUENCE: CALL DIPRES 

f 

REEEEENCES : DIMAIN 


ROUTINE NAME : DIPRNT 

DESGRi ION: This routine writes 'rri output default and actual values 

for all DICAL s(?igment control constants. 

CALLING SEQUENCE: GALL DIPRNT 

REFERENCE : DIMAIN 

FILE: NOUT 



u 


ROUTINE NAME; 
DESCRIPTION: 


ROUTINE NAME: DIPRQG 

DESCRIPTION: This is a prep cessor-generated routine that calls in 

main driving logic to perform the direct-irradiation 
computations , 

CALLING SEQUENCE; CALL DIPROG 
REFERENCES: TRASYS (Root Segment) 


1 
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ROUTINE NAME: 


DIPRTP 


DESCRIPTION: This routine prints albedo/planetary fluxes after each 

nodal computation. It may be overridden by the user if 
he desired to change the format. 

CALLING SEQUENCE: CALL DIPRTP 

KEY VARIABLES: INSHAD - Flag to determine if node is in the planet's 

shadow 

ICRD - Restart card number 

IN “ Current node sequence number 

REFERENCES : DIOUTP 


ROUTINE NAME : DIPRTS 

DESCRIPTION: This routine prints solar fluxes after each nodal 

computation. It may be overridden by the user if he 
desires to change the output format. 

CALLING SEQUENCE: CALL DIPRTS 

KEY VARIABLES: INSHAD - Flag to determine if node is in the planet's 

shadow 

ICRD - Restart card number 

IN “ Current node sequence number 

REFERENCES : DIOUTS 

ROUTINE NAME: DIPSHP 

DESCRIPTION: This routine determines possible shadowing surfaces 

between Node IN and the planet element. 

CALLING SEQUENCE; CALL DIPSHP (RADJ, RADI, POSJ, POSI, NST, JST, 

IN, JUMP) 

RADJ - Radius of sphere enclosing surface node 

RADI - Radius of sphere enclosing planet element 

POSJ - Position vector of surface sphere 

POSI - Position vector of planet element 

NST ~ Number of shadowing surfaces 

JST - Number of possible shadowers 

IN - Node sequence number being computed 

JUMP - Flag for using cylinder or cone technique 

KEY VARIABLES: ISEIAD ~ Array of possible shadowers 

JST “ Number of possible shadowers 

REFERENCES: DICALP, SFPL 
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ROUTINE NAME : DlPbliS 


DESCRIPTION: This routine ^Jatermines possible shadowing surfaces 

between Node IN and the sun. 

CALLING SEQUENCE: GALL DIPSUS (RADS, POS, SUKP, NSURF, NSHAD, IN) 

RADS - Radius of sphere encompassing Node IN 
POS - Position vector of center of encompass- 
ing sphere 

SUNP - Sun position vector 

NSURF - Number of shadowing surfaces 

NSHAD - Number of possible shadowing surfaces 

IN - Sequence number of node being computed 

KEY VARIABLES: ISHAD - Array of shadowing surfaces 

REFERENCES ; DICALS 


ROUTINE NAME : DIPTPR 

DESCRIPTION: This routine writes planetary fluxes obtained from 

the irestart tape to the output file. 

CALLING SEQUENCE: CALL DIPTPR 

REFERENCE : DIRDRQ , DIOUTP 

FILES : NOUT 

ROUTINE NAME : DIPTSR 

DESCRIPTION: This routine writes solar fluxes obtained from the 

restart tape to the output file. 

CALLING SEQUENCE: CALL DIPTSR 

REFERENCE : DIRDRQ 

FILES : NOUT 
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ROUTINE NAME: DIRGOS 


DESCRIPTION: This routine computes the direction cosines, given the 

rotation order and angles. 


CALLING SEQUENCE: CALL DIRCOS (II, JJ, KK, PHI, PSI, OM, TRAN) 


II 

JJ 

KK 

PHI 

PSI 

OMI 

TRAN(3,3) 


order of rotation for PHI, PSI, 
and OMI, respectively. 

Rotation about Z (Y to X = positive) 
Rotation about Y (X to Z = positive) 
Rotation about X (Y to Z = positive) 
Resultant direction cosine matrix 


REFERENCES: ORIENT, BUILDG 


ROUTINE NAME: DIRDRQ 

DESCRIPTION; This routine defines the restart request arrays and deter- 
mines the proper initial data based on the compute flags. . 
It also determines if shadowing is computed directly or by 
a table lookup. If computed, DIRDRQ defines arrays with 
all-shadowing surface data. 

CALLING SEQUENCE: CALL DIRDRQ 

KEY VARIABLES: ISOLFL - Flag for computing, storing, or zeroing the 

solar fluxes 

lALBFL - Flag for computing, storing, or zeroing the 
albedo fluxes 

IPLAFL - Flag for computing, storing, or zeroing the 
planetary fluxes 

REFERENCES : DIMAIN 

FILES: NDIR, NDI , NPLS , NOUT, NRAN, MSI, NRSO, NRTI, NRTO 


ROUTINE NAME : DIRPSP 

DESCRIPTION: This routine computes the radius, area, and position 

vector for the planet. 

CALLING SEQUENCE: CALL DIRPSP 

KEY VARIABLES : RADP - Radius of the planet 

AREAPL - Area of visible portion of planet 
P0SP(3) - Position vector of planet 

REFERENCES: DIMAIN, DIPRNT 
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ROUTINE NAME: 


DIRTI 


DESCRIPTION: This routine reads flux data from the temporary 

restart tape. 

CALLING SEQUENCE: GALL DIRTI 

REFERENCE : DIRDRQ 

FILES: NOUT, NDI, NRTI 

ROUTINE NAME: DIRTOP 

DESCRIPTION: This routine writes planetary fluxes to the temporary 

restart file. 

CALLING SEQUENCE; CALL DIRTOP 
REFERENCE: DIOUTP 

FILE : NRTO 

ROUTINE NAME : DIRTOS 

DESCRIPTION: This routine writes solar flux data to the temporary 

restart tape. 

CALLING SEQUENCE: CALL DIRTOS 

REFERENCE : DIOUTS 

FILE: NRTO 

ROUTINE NAME : DIRT? 

DESCRIPTION: This routine reads planet-element factors from a 

previously stored DI call and defines the QDR and QDP 
arrays based on known planet factors. 

CALLING SEQUENCE: CALL DIRTP 

REFERENCES : DITYPE 


FILES : NPLS 


ROUTINE NAME: 


DISFTP 


DESCRIPTION: This routine reads shadow factor tables from the NPLS 

file for use in computing direct fluxes. 

CALLING SEQUENCE: CALL DISFTP (1ST) 

REFERENCES: DICALP, DICALS, DIMAIN 

FILES : NPLS 


ROUTINE NAME : PI TIME 

DESCRIPTION: This routine computes the true anomaly from time- 

dependent orbital characteristics. 

CALLING SEQUENCE: CALL DITIME (TIME, PER, ECC, TRU) 

TIME - Input orbital time 
PER - Orbital period 
ECC - Eccentricity of orbit 
TRU - Output true anomaly 

REFERENCES : DILOC , STFAQ 
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ROUTINE NAME: DITRS3 


DESCRIPTION: This routine converts points from the surface coordinate 

system to the central coordinate system, 

CALLING SEQUENCE: CALL DITRS3 (X, Y, Z, A, B, C, RX, RY, RZ, TRAN) 

Y, Z - Output components of vector in CCS 

Aj B, C - Components of vector to be transformed 

RX, RY, RZ - Position vector of SCS origin in CCS 
TRAN - Direction cosine matrix 


REFERENCES: DIELEM 


ROUTINE NAME: DITRS4 

DESCRIPTION: This routine converts points in the central coordinate 

system to the surface coordinate system. 

CALLING SEQUENCE; CALL DITRS4 (X, Y, Z, A, B, C, RX, RY, RZ, TRAN) 

X, Y, Z - Output components of vector in SCS 

A, B, C - Components of vector to be transformed 

RX, RY, RZ - Translation vector, in SCS 
TRAN - Direction cosine matrix 

REFERENCES; DILOC, DIL0C2 


ROUTINE NAME : DITTP 

DESCRIPTION: See program listing for definition. 

ROUTINE NAME : PI TYPE 

DESCRIPTION: This routine determines the type of computation and calls 

the proper routine. 

CALLING SEQUENCE; CALL DITYPE (IRSPH) 

KEY VARIABLES: PLTYPE - Flag indicating tnat planetary data were pre- 

viously defined. 

REFERENCES : DIMAIN 
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routine NAME; DIVI^PL 


DESCRIPTION; This routine computes elemental planetary factors without 
shadowingj and forms a sum over all the elements, 

CALLING SEQUENCE: CALL DIVWPL (NTOTP,AREASF,NOT,FRACT,FRACTR,FRACTP) 

REFERENCES; DICALP 


ROUTINE NAME : DIWSN 

DESCRIPTION: This routine computes unshadowed elemental view factors to 

e sun and forms a sum over all the elements. 

CALLING SEQUENCE; CALL DIWSN ( AREAS F, NOT, FRAGT) 

REFERENCES ; DICALS 


ROUTINE NAME; DUPNCK 

DESCRIPTION: This is a user-called routine in the Operations Data block 

that determines if a configuration has duplicate node 
numbers. It aborts if duplicates are found. 

CALLING SEQUENCE: CALL DUPNCK 

REFERENCES; Operations Data block 


ROUTINE NAME; DRCALS 

DESCRIPTION; This is the main routine used in calculating the specularly 
reflected components of incident solar flux. 

CALLING SEQUENCE: CALL DRCALS (MIRROR, KN, NSURFS) 

MIRROR - Surface sequence number of current specular 
surface 

KN ~ Index ranging over the number of specular 

surfaces and indicating the current specular 
surface 

NSURFS - Total number of surfaces, plus images of 
surfaces , in MIRROR 


REFERENCES; DRMAIN 
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ROUTINE NAME: DRDATA 


DESCRIPTION: This routine is called hy the user to set up necessary 

control constants to execute program segment, 

CALLING SEQUENCE: CALL DRDATA (NSTP, DIACCLS) 

KEY VARIABLES; Reference Users Manual, Appendix D. 

Reference; ODPROG (User's Operations Data) 


ROUTINE NAME; DRELEM 

DESCRIPTION: Given the number of elements in each direction on a 

node, this routine calculates the area and position 
vectors for each element. 

CALLING SEQUENCE: CALL DRELEM (ILF, DATA, TRAN, RX, RY, RZ, NTOT) 

ILP - Surface type 

DATA - Nodal dimension parameters (ALPHA, 

BMIN, BMAX, GMIN, and GMAX) 

TRAN “3x3 matrix of direction cosines 

RX, RY, RZ - Translation vector 

NTOT - Total number of elements required on 

node 

KEY VARIABLES: NEST - Counter for number of elements on a node 

(maximum allowable = 100) 

SFAV - Array of elemental surface area vectors 
SFPV - Array of elemental position vectors 

REFERENCES ; DRCALS 


ROUTINE NAME; DRELSL 

DESCRIPTION: Given the total number of elements on a node, this 

routine calculates the number of elements in each 
direction so as to make them as nearly square as 
possible, 

CALLING SEQUENCE: CALL DRELSL (NB, NG, ILP, DATA, NTOT) 

NB ~ Number of elements in the beta direction 

NG - Number of elements in the gamma direction 

ILP - Surface type 

DATA “ Nodal dimension parameters (ALPHA, BMIN, 
BMAX, GMIN, and GMAX) 

NTOT - Minimum number of elements to be distributed 
over node 

REFERENCES : DRELEM 
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ROUTINE NAME: 


DREND 


DESCRIPTION: This routine writes direct incident fluxes QDS, QPR, and 

QDP to the NDI disk file (labeled 6HIMAGEQ) to indicate 
that specular components are included. 

CALLING SEQUENCE; CALL DREND 

KEY VARIABLES: QDS “ Incident solar flux 

QDR - Incident albedo flux 
QDP - Incident planetary flux 

REFERENCES : DSMAIN 

FILES: NDI - Disk file for storing direct incident fluxes 

NRSO 


ROUTINE NAME: DRENDS 

DESCRIPTION; This is a user routine that enables the user to inter- 
vene after calculating the incident solar fluxes, 

CALLING SEQUENCE: CALL DRENDS 

REFERENCES : DRMAIN 


ROUTINE NAME: DRGTST 

DESCRIPTION: This is an integer function routine that tests a 

given value, GT, to determine if it falls in the 
range GN<GT<GX. If true, the function value is 0 
(zero); if false, the function value is 1 (one). 

CALLING SEQUENCE: DRGTST (GN, GX, GT) 


GN - Lower bound 
GX - Upper bound 
GT - Value to be tested 


REFERENCES; DRSHAD 


ROUTINE NAME: 


DRIMAG 


DESCRIPTION: This routine images the solar vectors as well as all 

shadowing surfaces, and writes the results on scratch 
file NSCR2 for each specular surface* 

CALLING SEQUENCE: CALL DRIMAG (NST) 


NST 

NSURFS 

IFS, IKS, 

SUNPVT 


- Number of active surfaces 

- Number of active surfaces plus 
number of images in any given 
specular surface 

PR, DSTR, DIMS, PSH, TSTR - 

Shadowing surface description 
parameters 

- Solar vector 


REFERENCES ; DRMAIN 
FILES; NSCR2 - Scratch file 


ROUTINE NAME ; DRMAIN 

DESCRIPTION: This routine contains the main driving logic for computing 

direct incident fluxes, including specular components. 

CALLING SEQUENCE: CALL DRMAIN 

KEY VARIABLES : ISPEC - Array of surface sequence numbers for specular 

surfaces 

SUNPVT - Solar vector 

PLDC - Matrix of direction cosines to transform vectors 
in the planet-oriented VCS to the user-defined 
VCS 

REFERENCES : DRPROG 

FILES: NEAN - Random access file 

NSCR2 - Scratch file - 

NDI, NOUT, NRSI, NRTO I 
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ROUTINE NAME: 


DROUTP 


DESCRIPTION: This is a user routine that can be replaced, to change 

the form of output for albedo and plauet-ary incident 
fluxes. The normal call is to DRPRTP, which prints 
the albedo and planetary fluxes, 

CALLING SEQUENCE: CALL DROUTP 

REFERENCES: DRMAIN 


ROUTINE NAME: DRQUTS 

DESCRIPTION; This Is a user routine that can be replaced to change 
the form of output for incident solar fluxes. The 
normal call is to DRPRTS, which prints the solar fluxes. 

CALLING SEQUENCE: CALL DROUTS 

REFERENCES : DRMAIN 


ROUTINE NAME: DRPOSI 

DESCRIPTION; This routine transforms a vector in the ICS, BCS, or 
CCS to a vector in the SCS of a specular surface, 
negates the Z component, and transforms it back, into 
the ISC, BCS, or CCS. 

CALLING SEQUENCE: CALL DRPOSI (X, Y, Z, TRAN) 

X, Y, Z - Vector components in ICS, BCS, or CCS 
TRAN - Matrix of direction cosines 

KEY VARIABLES: A, B, C - Vector components in the SCS 

REFERENCES : DRIMAG 
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ROUTINE NAME: 


DRPRTP 


DESCRIPTION: This routine prints albedo and planetary incident 

fluxes after each nodal computation. The user can 
override this routine to change the output format if 
he desires. 

calling SEQUENCE: CALL DRPRTP 

KEY VARIABLES; IN - Current node sequence number 

INSHAD - Flag to indicate if vehicle is in planet 
shadow 


REFERENCES : DROUTP 

FILES: NOUT - System output file 


ROUTINE NAME : DRPRTS 

DESCRIPTION: This routine prints the solar incident flux after 

each nodal computation. The user can override this 
routine to change the output format if he desires. 

CALLING SEQUENCE: CALL DRPRTS 

KEY VARIABLES; IN - Current node sequence number 

INSHAD - Flag to indicate if vehicle is in planet 
shadow 

REFERENCES : D ROUTS 

FILES: NOUT - System output file 
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ROUTIHE ME: ORPSHS 


DESGRIETIOH: This routine determines all possible shadowing surfaces 

between Node IN and the image of the sun in specular 
surface MIRROR, 


cabling SEQUENCE: 


CALL DRPSHS (RADS, POS, SUNP, NSURP, NSS, NSHAD, IN, 
MIRROR) 


RADS 
POS , . 
SUNP 

NSURF 

NSS 

NSHAD 

IN 

MIRROR 


- Radius of sphere enclosing Node IN 
^ Position vector of enclosing sphere 

- Image of solar vector in specular surface 
MIRROR 

- Total number of shadowing surfaces 

- Total number of shadowiiig surfaces plus images 
of shadowing surfaces in MIRROR 

- Number of possible shadowing surfaces 

- Sequence number of current node 

- Surface sequence number of current specular surface 


REPERENCES: DRCALS 


ROUTp^ NAME: 


DRRDRO 


DESCRIPTION: This routine initializes the direct incident flux 

arrays QDS, QDR, and QDP to the flux values calculated 
by segment DICAL, 


CATLING SEQUENCE: CALL DRRDRQ (ICST) 


KEY VARIABLES; Ot»S - Incident solar flux array 
QDR - Incident albedo flux array 
QDP - Indicant planetary flux array 

REFERENCES ; DRMAIN 


FILES: NOUT - System output file 

NDI - Disk file for storing direct incident fluxes 


ROUTINE NAME; DRRSI 

DESCRIPTION: This routine reads real-body effect direct fluxes from 

the restart tapes. 

CALLING SEQUENCE: CALL DIRS I (ICOMP, KN) 

ICOMP - Compute/no compute flag 
KN - Node number on RTl tape record 

REFERENCE: DRMAIN 

FILES: NRSI, NRTI 
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ROUTINE NAME; DRSHAfl 

DESCRIPTION; This routine calculates the elemental shadowing, 

where an element is either completely shadowed or not 
shadowed at ail. Shadowing is reduced by the transmis- 
sivity for semitransparent shadowing surfaces. 

CALLING SEQUENCE; CALL DRSHAD (ILKI, RX, RY, RZ, NE, WA, IN, I, NSS, RS, 

MIRROR) 

RX, RY 
NE, WA 

IN 
I 

NSS 
RS 

MIRROR 


RE FERENCES ; DRCAL S 


ROUTINE NAME: DRTRS3 

DESCRIPTION: This routine transforms points in an SCS to points in 

the ICS, BCS, or CCS. 

CALLING SEQUENCE; CALL DRTRS3 (X, Y, Z, A, B, C, RX, RY, RZ, TRAN) 

X, Y, Z - Coordinates of a point in the ICS, BCS 
or CCS 

A, B, C - Coordinates of the point in the SCS 
RX, RY, RZ - Components of the SCS origin position 
vector in the ICS, BCS, or CCS 
TRAN - Matrix of direction cosines 

REFERENCES: DRELEM, DRIMAG, DRSHAD 


, RZ - Components of solar position vector 

- Elemental shadowing factors for IR 
and solar fluxes 

- Surface sequence number 

-/ Sequence number of element 

- Number of possible shadowing surfaces 
“ Square of the magnitude of the solar 

I vector 

- Sequence number of current specular 
/ surface 
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ROUTINE NAME: DRIES 4 


DESCRIPTION; This routine transforms points in the ICS, BCS, or CCS 
to points in an SCS. 

CALLING SEQUENCE: CALL DRTRS4 (X, Y, Z, A, B, C, RX, RY, RZ, TRAN) 

X, Y, Z “ Coordinates of a point in the SCS 
A, B, C - Coordinates of the point in the ICS, 
BCS, or CCS 

RX, RY , RZ - Components of the ICS, BCS, or CCS 
origin position vector in the SCS 
TRAN - Matrix of direction cosines 

REFERENCES; DRIMAG 


ROUTINE NAME : DRVWSN 

DESCRIPTION: This routine calculates the unshadowed form factor from 

a node to the sun. 

calling SEQUENCE: CALL DRVWSN (AREASF, NOT, FRACT) 

AREASF - Surface area of current node 
NOT - Flag to indicate whether node can "see" 
the sun 
= 0 Can see 
= 1 Cannot see 

FRACT - Unshadowed node-to-sun form factor 

references ; DRCALS 


ROUTINE NA:IE: EARTHD 

DESCRIPTION: This routine defines all Earth-oriented constants 

needed by the DI link. 

CALLING SEQUENCE: CALL EARTHD 


KEY VARIABLES; PRAD 

- 

Planet Radius 


SOL 

- 

Solar Constant 


PALB 

- 

Planet Albedo Factor 


WDS 

- 

Planet Darks ide Emissive 

Power 

WSS 

- 

Planet Subsolar Emissive 

Powe r 

GIlAV 

- 

Gravitational Constant 
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REFERENCES: ORBITl, 0RBIT2, DITTP 


ROUTINE NAME: ENDTP 

DESCRIPTION: This routine writes a pseudo end file on the BCD unit 

and iiiforms the user as to type of data the unit contains. 

CALLING SEQUENCE; CALL ENDTP (NTYPE) 

NTYPE “ Flag to indicate type of end message to write 

REFERENCES; QOSAVE, QOCOMB, RCMAIN 


ROUTINE NAME : FFAREA 

DESCRIPTION: This function routine computes the area of a surface, 

given the surface properties. 

CALLING SEQUENCE: AG = FFAREA (ILP, ALPH, BMIN, BMAX, GMIN, GMAX) 

Ac - Surface area 
ILP - Surface type 

ALPH, BMIN, BMAX, GMIN, GMAX - Surface dimensions 
REFERENCES; FFEXPN, FFXPNI 


ROUTINE NAME : FFCAL 

DESCRIPTION: This is the main computation routine for computing the 

firm factor. It determines the proper number of elements 
and computes position area vectors and the final form 
factor, including shadowing. 

call: G SEQUENCE: CALL FFCAL 

REFERENCES: FFPRE 
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ROUTINE NAME: 


FFCALI 


DESCRIPTION: This routine is analagous to subroutine FFCAL. It is 

used for form calculation when reversing the order of 
the i - j node pair due to area difference. 

CALLING SEQUENCE: CALL FFCALI 

REFERENCE: FFPRE 


ROUTINE NAtlE: FFDATA 

DESCRIPTION: This is a user^-callable routine to define the control 

variables in the FF link. 

CALLING SEQUENCE: CALL FFDATA (ACC, ACCS, IRNOSH, RATL, RMIN, IPRNT, 

IPUN, IFlfAC; 

See Users Manual, Appendix D, for definition 
REFERENCES: Operations Data block 


ROUTINE NAME: FFELEM 

DESCRIPTION: This routine computes elemental position and area vectors, 

given the number required in each nodal direction. 

CALLING SEQUENCE: CALL FFELEM (NB, NG, ILP , IC , DATA, POS , ARA, 

TRAN, RX, RY, RZ) 


NB 

NG 

ILP 

rc 

DATA(5) 
POS(500, 3) - 
ARA(500, 3) - 

TRAN(3, 3) - 

RX, RY, RZ - 

REFERENCES: FFCAL, FFEXPN . FFCA 


Number of elements in the beta 
direction 

Number of elements in the gama 

direction 

Surface type 

Position to start storing element in 

the POS and ARA arrays 

Array of node dimensions 

Array of element position vectors 

Array of normal vectors (magnitude 

& area) 

Direction cosine matrix of node 
Position vector of SGS origin in 
the CCS 
J, FFXPNI 
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ROUTINE NAME : FFELSL 


DESCRIPTION: This routine uses the total number of elements on a node 

to compute the number of elements in various directions 
so that the elements are as nearly square as possible. 


CALLING SEQUENCE: CALL FFELSL (NB, NG, ILP , DATA, NTOT) 


NB 

NG 

ILP 

DATA(S) 

NTOT 


- Computed number of elements in the beta 
direction 

- Computed number of elements in the gamma 
direction 

- Surface type 

- Array of surface dimensions 

- Total number of elements to subdivide the 
node into 


REFERENCES: FFCAL, FFEXPN, FFCALI, FFELSL, FFXPNI 


ROUTINE NAME: FFEND 

DESCRIPTION: This routine is called prior to the termination of the 

FF link to write the end-of-£iles to the internal disk 
units and restart files. 

CALLING SEQUENCE: CALL FFEND 

REFERENCES : FFMAIN 

FILES: NFF, 


ROUTINE NA^E : FFEQIV 

DESCRIPTION: This routine retrieves equivalent form factors from 

data files when equivalent form factor flags are found 
in the form factor request matrix. 

CALLING SEQUENCE; CALL FFEQIV (FTEMP) 

FTEMP - Array for use in reading random-access file 


REFERENCES : FFOUT 

FILES: NOUT, NSCR3 
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ROUTINE NAME : FFESUM 


DESCRIPTION: This routine outputs form-iactor sums for all nodes. 

CALLING SEQUENCE: GALL FFESUM 

KEY VARIABLES: SUM - Array of form-factor suras 

REFERENCES: FFEWD 


ROUTINE NAME : FFEXPN 

DESCRIPTION; This routine computs a node-pair form factor with sh 
shadowing, using the subnode technique. Logic is 
entered when the separation distance along the node 
varies by more than FFRATL. 

CALLING SEQUENCE: CALL FFEXPN (POSI, ARAI, FE, FA, RATI, RATJ, 

JFLAG, JST, NEI, NEJ, NSUR, lADDR) 

POSI - Position vector, storing elements on 
node i and j 

ARAI - Area vector, storing elements on node 
i and j 

FE - Resultant IR form factor 
FA - Resultant solar form factor 
RATI - Ratio of rate of change along node i 

RATJ - Ratio of rate of change along node j 

JFLAG - Flag indicating that maximum number of 
elements has been exceeded 
JST - Number of possible shadowers 
NEI - Average number of elements on node i 

NEJ - Average number of elements on node j 

NSUR - Average number of shadowing surfaces 
lADDR - Maximum time return sequence (not used) 

REFERENCES; FFCAL 
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ROUTINE NAME : FFGTST 


DESCRIPTION: This integer function routine determines if gamma 

falls in allowable ranges and sets the functional 
value to 1 if it does and to 0 if it doesn't, 

CALLING SEQUENCE: FFGTST (GN, GX, GT) 

GN - Minimum gamma 
GX - Maximum gamma 
GT - Gamma in question 

REFERENCES; FFSHD 


ROUTINE NAME: FFHEAD 

DESCRIPTION: This routine prints FFcontrol values in summary form on 

an output file. 

CALLING SEQUENCE: CALL FFHEAD 

REFERENCES : FFMAIN 


ROUTINE NAME; FFMAIN 

DESCRIPTION: This is the main controlling logic in the FP link. 

It checks for restart da:a and directs the logic flow. 

CALLING SEQUENCE; CALL FFMAIN 

REFERENCES; FFPROG 
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ROUTINE NAME; 


FFMINR 


DESCRIPTION; This routine determines if the computed form factor is 
less than the control value FFMIN. If less, it rede- 
fines the FFVALS and FFVALI array elements to zero. 

CALLING SEQUENCE: CALL FFMINR 

REFERENCES : FFOUT 


ROUTINE NAME : FFNDP 

DESCRIPTION; This routine provides the user with a node array 

punched in form factor and flux data block format. 

CALLING SEQUENCE: CALL FFNDP 

REFERENCE; ODPROG (Users Operation Data) 

FILES : NOUT 

ROUTINE NAME ; FFOUT 

DESCRIPTION; This is a user-replaceable routine calling for print 
(as a minimum) and punch options. 

CALLING SEQUENCE; CALL FFOUT 

RE FERE NCE S ; FFMAIN 
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ROUTINE NAME 


FFPCH 


DESCRIPTION: This routine punches restart form factors if the FFPNCH 

option is used. The restart form factors are punched 
in a format acceptable to the Header Form Factor Data 
block. 

CALLING SEQUENCE: CALL FFPCH 

REFERENCES: FFOUT 


ROUTINE NAME: FFPRE 

DESCRIPTION: This is a user-definable routine to determine the type 

of computation. The program normally calls FFCAL to 
use a finite-element technique. 

CALLING SEQUENCE: CALL FFPRE 

REFERENCES : FFMAIN 


ROUTINE NAME: FFPROG 

DESCRIPTION: This is a preprocessor-generated routine that defines 

necessary common blocks and calls into the main logic. 

CALLING SEQUENCE: CALL FFPROG 

REFERENCES: TIU^YS (root segment) 

ROUTINE NAME: FFPRT 

DESCRIPTION; This routine prints nonzero computed and predefined 
form factors on an output file. 

CALLING SEQUENCE: CALL FFPRT 

REFERENCES; FFOUT 
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ROUTINE NAME: FFPSHD 


DESCRIPTION: This routine computes the number of possible shadowing 

surfaces between node pairs. 

CALLING SEQUENCE: CALL FFPSHD (RADJ, RADI, POSJ, POSI, NST, JST, 

IN, JN) 

RADJ - Radius of sphere enclosing node j 

RADI - Radius of sphere enclosing node i 

POSI ~ Position vector to center of sphere 
enclosing node i 

POSJ - Position vector to center of sphere 
enclosing node j 

NST - Number of shadowing surfaces 
JST “ Number of possible shadowers 
IN - Sequence number of node i 

JN - Sequence number of node j 

REFERENCES: FFCAL, FFEXPN, FFCALI, FFXPNI 


ROUTINE NAME : FFRDIN 

DESCRIPTION: This routine determines if there is a form factor re- 

start, and defines the shadowing data. 

CALLING SEQUENCE: CALL FFRDIN (ISHDQ, IRSC) 

ISHDO - Flag indicating shadower-only surfaces 
1 - Shadower-only surfaces 
0 - Not shadower-only surfaces 


REFERENCES : FFMAIN 


ROUTINE NAME: FFRDRQ 

DESCRIPTION; This routine reads the restart file of form factors 
and defines know values, or sets a flag to compute 
values for the current row. 

CALLING SEQUENCE; CALL FFRDRQ 

lEOFFR - End-of-file flag for the NFFR unit 

FFVALI - Array of values/ flags to compute IR form factors 

FFVALS - Array of values/flags to compute solar form factors 


REFERENCES ; FFMAIN 

FILES; NOUT, NRSI, NRTI 
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ROUTINE NAHE: EERIN 


DESCRIPTION; This routine reads form factor information from the 
form factor data block via file NFFR. 

CALLING SEQUENCE: CALL FFRIN 

REFERENCE ; FFRDIN, FFRDRQ 

FILE ; NODT 


ROUTINE NAME ; FFRQM 

DESCRIPTION; This routine outputs to internal file NFF a row of 

form factors to be used in matrix form by other links. 

CALLING SEQUENCE: CALL PFROW 

REFERENCES ; EFMAIN, FFRDIN 

FILES: Writes NFF 

NRSO, NRTO, NSCR3 


ROUTINE NAME: FE RPSN 

DESCRIPTION: Given the nodal dimensions, this routine computes the 

radius of an encompassing sphere and the position 
vector to the center of the sphere. 

CALLING SEQUENCE: CALL FFRPSN (RADN, POSN, ILK, DATA, BETA, GAMMA, 

DB, DG) 

RADN Computed radius of sphere enclosing node 

P0SNC3) - Vector to center of sphere 
ILK - Surface type 
DATA(5) - Dimensions of node 

BETA - Length to center along beta direction 
GAMMA - Length to center along gamma direction 
DB - Absolute value of deviation from gamma 

DG - Absolute value of deviation from beta 

REFERENCES: FFEXPN, FFXPKI 


ROUTINE NAME: FFRSMR 


DESCRIPTION; This routine writes the form factor sums to the output 
file when form factors are real from the master restart 
tape. 

CALLING SEQUENCE: CALL FFRSMR 

REFERENCE: FFRDIN 

FILE : NOUT 


ROUTINE NAME: FFRSUM 

DESCRIPTION; This routine is called after completing each row, to 
output the FF sum and time. 

CALLING SEQUENCE: CALL FFRSUM 

REFERENCES: FFMAIN 

ROUTINE NAME: FFRTR 

DESCRIPTION: This routine reads form factor data from the temporary 

restart tape. 

CALLING SEQUENCE: CALL FFRTR 

KEFERENCE: FFRDIN 

FILE: NRTI, NOUT 
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ROUTINE NAME: 


FFSHD 


DESCRIPTION: This routine computes elemental shadowing between node 

IN and node JN. An element is either completely shadowed 
or not shadowed at all. Shadowing is reduced by the 
transmissivity for semitransparent shadowing surfaces. 

CALLING SEQUENCE: CALL FFSHD (ILKI, ILKJ, EX, RY, RZ, POSI, WE, WA, JST, 

IN, JN, I, NSS, RS) 


ILKI 

RX, RY, RZ 

POSI 
WE, WA 

IN, JN 
I 

NSS 

RS 


Type of surface for IN 

Components of vector from element on 

node i to element on node j 

Array of elemental position vectors 

Elemental shadowing factors for IR and 

solar form factors 

Surface sequence numbers 

Element sequence number 

Number of possible shadowers 

Square of the magnitude of the element- 

to-element connecting vector (RX^ + RY^ 

+ RZ^) 


REFERENCES : 


FFCAL, FFEXPN, FFCALI, FFXPNI , FFPCL 


ROUTINE NAME: FFTMCK 

DESCRIPTION; This routine calculates the time remaining in the run 
and compares it with the estimated time for the next 
calculation. It there is insufficient time remaining, 
the run is aborted. 

CALLING SEQUENCE: CALL FFTMCK (lADDR, MAXLFT) 

lADDR - Not used 

MAXLET - Estimated time required for next calculation 
M - Time remaining in run 

REFERENCES: FFCAL, FFEXPH FFCALI, FFXPNI 
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ROUTINE NAME : FFTRS3 


DESCRIPTION: This routine transforms a point in an SCS to points 

in the ICS , BCS , and CCS . 

CALLING SEQUENCE: CALL FFTRS3 (X, Y, Z, A, B, C, RX, RY, R2, TTAN) 

Xj Y, Z “ Coordinates of a point in the ICS, 

BCS , or CCS 

A, B, C - Coordinates of the point in the SCS 
RX, El, RZ - Components of the SCS origin position 
vector in the ICS, BCS, or CCS 
TRAN - Matrix of direction cosines 

REFERENCES: FFELEM, FFEXPN, FFUNIT, FFUNTI, FFXPNI 


ROUTINE NAME: FFUNIT . FFUNTI 

DESCRIPTION: These routines calculate unshadowed form factors 

between quadrilaterals (and triangles) using the 
Nusselt sphere method with a single integration 
over the I-node, 

CALLING SEQUENCE: CALL FFUNIT (NTI, POSI, ARAI, FRACT, NTOTI, 

NTOTJ, NOT) 

KEY VARIABLES: See FPVIEW 

REFERENCES: FFCAL, FFCALI 
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ROUTINE NAME : FFVIEW 


DESCRIPTION: This routine calculates the unshadowed form factor between 

node IN and node JN. It als.o calculates RATI and RATJ, 
which will later be compared with the user-input variable, 
FFRATL, to see if nodes should be expanded into subnodes 
for better accuracy. 


CALLING SEQUENCE: CALL FFVIEW (NTI, NTJ, POSI, ARAI, NTOTI, NTOTJ, ARI, 

ARJ, NOT, FRACT, RATI, RATJ) 


NTI 

NTJ 

POSI 

ARAI 

NTOTI, NTOTJ 

ARI, ARJ 
NOT 


FRACT 

RATI, RATJ 


Initial number of elements on node IN 
Total initial number of elements on node 
IN and node JN 

Array of elemental position vector 

Array of elemental area vectors 

Number of elements required on node 

IN and node JN 

Areas of node IN .■’nd node TN 

Flag to indicate whether nc IN can 

see node JN 

= 0 Gan see 

= 1 Cannot see 

Unshadowed node-to-node form factor 
Ratios of indicated maximum elements 
required in NTOTI and NTOTJ 


REFERENCES: FFCAL, FFCAT.T 


I 
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ROUTINE NAME; FFVWT 


DESCRIPTION; This routine calculates unshadowed form factors between 
subnodes and determines the number of elements required 
based on a weighted-average criterion. 

CALL FFVWT (NTI, NTJ, POST, ARAI , NTOTI, NTOTJ, ARI, 

ARJ, NOT, FRACT) 

NTI - Number of elements on subnode of node IN 

NTJ - Total number of elements on subnode IN 

and subnode JN 

POSI “ Array of elemental position vectors 

ARAI ~ Array of elemental area vectors 

NTOTI, NTOTJ - Number of elements required on subnodes 
of node IN and node JN 
ARI, ARJ - Areas of subnodes 

NOT - Flag to indicate whether subnode IN can 

"see" subnode JN 
= 0 Can see 
= I Cannot see 

FRACT - Unshadowed subnode-to-subnode form factor 

REFERENCES ; FFEXPN , FFXPNI 


ROUTINE NAME; FFXPNI 

DESCRIPTION: This routine is analogous to subroutine FFEXPN and is 

used when t!ie order of the i - j pair is reversed due 
tc area difference. 

CALLING SEQUENCE: SEE FFEXPN 

REFERENCE: FFCALI 
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ROUTI>'E NAME 


FINDS! 


DESCRIPTION: Given the unit, step number, label 1, and label 2, 

this routine positions the file at the proper point 
or returns and end-of-file flag indicating it was 
unable to find the data. 

CALLING SEQUENCE: CALL FINDS! (NUNIT, ISTEP, LABl, LAB 2, lEOF) 

NUNIT - Unit to be searched 

ISTEP - Program step number needed 

LABl - Subidentifier 1 

LAB2 - Subidentifier 2 

lEOF - Flag indicating if data were found 

= 2HN0, No end-of-file encountered 
(found data) 

= 3HYES End-of-file found (no data) 

REFERENCES: FFRDIN,SFMAIN,SFRDIN,DIRDRQ,GBMAIN^ AQMAIN, QOMAIN, 

QOCMBN, PLLOAD, FNDPLP^ RCMAIN.ADSURF, AQCMPR, CMINIT, 
CMMAIN, DIMAIN, DIRTI, DRRDRQ, GBAPRX, PLCMBN, RBRDIN, 
RCINIT,STPAQ 

FILES: All internal data files 

ROUTINE NAME : FLIP 

DESCRIPTION; Given a working storage area, this routine flips a matrix 
in row order on unit NSl to a matrix in column order and 
writes it to unit NS2. 

CALLING SEQUENCE; CALL FLIP (DATA, NV, NTIME, ICOL, NSl, NS 2) 

DATA “ Working storage area 

NV - Number of rows stored at one time 

NTIME - Number of rows 

ICOL - Number of columns 

NSl - Unit containing the input matrix 

NS 2 - Unit containing the output matrix 

REFERENCES : FNDFLP 


ROUTINE NAME: FNDEXP 

DESCRIPTION: Given a floating-point number, this routine converts 

it to scientific notation. 

CALLING SEQUENCE: CALL FNDEXP (X, BASE, lEXP) 

X 

BASE - 0.0 ^BASE <_ 1.0 

lEXP - Exponent such that X = BASE**IEXP 

REFERENCES: QOSBCD, QOAVGS 
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ROUTINE NAME 


FNDFLP 


DESCRIPTION: This routine reads the overall incident solar, albedo, 

planetary, and total fluxes from a disk generated by 
AQCAL or DICAL and stores them on a scratch file by 
sequential orbit points. It also builds titles as a 
function of user input data. 

CALLING SEQUENCE: CALL FNDFLP (DATA, SOL, ALB, PLA, TOT, NV, LCl, LC2, 

LC3, NTIME, ITIME, TIME, NUNIT) 

DATA - Working storage area where DATA (1) = SOL (1) 

SOL Working array DATA (NTIME+1) = ALB(l) 

ALB - Working array DATA (2NTIME+1) = PLA(l) 

PLA - Working array DATA (3NTIME+1) = TOT(l) 

TOT - Working array 

NV - Number of data items to be plotted 

LCl - First control character - I (incident) or A (absorbed) 
LC2 “ Second control character “ F (flux) or R (rate) 

LC3 - Third control character - S (solar), A (albedo), 

P (planetary), T (total) and ALL (all) 

NTIME - Number of orbit positions 

ITIME, TIME - Doublet array containing time and step number 
NUNIT “ Plot output unit 


P^FERENCES ; PLLOAD 
FILES: NSCRl, NSCR2, NOUT 


ROUTINE NAME: FRAME 

DESCRIPTION: This is a plotting routine to advance to a new frame. 

CALLING SEQUENCE: CALL FRAME 

REFERENCES: NPMAIN, OPHAIN, PLDRIV 
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ROUTINE NAME: FRAMEC 


DESCRIPTION: This routine determines the number of frames of plot 

data generated to this point. 

CALLING SEQUENCE; CALL FRAMEC (N, I) 

N “ Number of frames 
I - Dayfile message flag 

= 0 Do not print dayfile message 
= 1 Print dayfile message 

REFERENCES: NPMAIK, OPMAIN, PLDRIV 


ROUTINE NAME: GBAPRX 

DESCRIPTION; This is a user-called routine that computes approxi- 
mately gray-body factors. 

CALLING SEQUENCE; See User's Manual, Appendix D. 

REFERENCE: ODPROG (Users Operations Data Block) 

FILES: NFF, NGBIR, NGBSO, NOUT, NSCRl, NSCR2, NSCR3 


ROUTINE NAME: GBDATA 

DESCRIPTION: This is a user-called routine to define the step 

number and type of gray-bodies to compute. 

CALLING SEQUENCE: CALL GBDATA (WBAND, MODFF, NTYP) 

NSFF - Step number containing form factors 
WBAND - Waveband to compute gray-bodies for 
2HIR, 3HSOL, and 4HB0TH 

REFERENCES; Operations Data Block 
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ROUTINE NAME: GBEND 


DESCRIPTION: This routine can be replaced by the user to intervene 

after completion of the gray-body link. 

CALLING SEQUENCE: CALL GBEND 

REFERENCES: GBMAIN 

ROUTINE NAME ; GBHEAD 

DESCRIPTION: This routine prints the input parameters defined by 

the user, together with input options. 

CALLING SEQUENCE: CALL GBHEAD 

REFERENCES : GBMAIN 
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routine name : GBINV 


DESCRIPTION: This routine computes the inverse of the modified form 

factor matrix stored either in core or in large blocks on disk 

CALLING SEQUENCE; CALL GDINV (A, B) 

A - Working block 1 
B - Working block 2 


REFERENCES; GBSCFA 


ROUTINE NAME; GBMAIN 

DESCRIPTION: This is the main controlling routine in the gray-body 

link that directs main control and defines the block 
size . 

CALLING SEQUENCE; CALL GBMAIN 
REFERENCES ; GBPROG 


ROUTINE NAME; GBP RE 

DESCRIPTION: This routine can be replaced by the user to intervene 

prior to the call to the invert routine. 

CALLING SEQUENCE: CAI.L GBP RE 

REFERENCES: OBMAIN 


ROUTINE NAME: GBPROG 

DESCRIPTION; This is a preprocessor-generated routine that defines 

the necessary commons and calls into the main gray-body 
logic. 

CAI.LLNG SEQUENCE; CALL GBPROG 
REFERENCES; TRASYS (root segment) 
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ROUTINE NAME 


GBRSI 


DESCRIPTION: This routine reads Rrav-body factor data from the 

Master restart tape and writes it to the rtSCRl file 
and . the master restart output tape. 

CALLING SEQUENCE- CALL GBRSI (TYPE, ICOMP) 

TYPE - type of gray-body factors on NRSI 
(solar, infrared) 

ICOMP - compute/no compute flag 

REFERENCE : GBMAIN 

FILES; NRSI, NRSO, NSCRl, NGBIR, NGBSO, NOUT 


ROUTINE NAME: GBRTI 

DESCRIPTION: This routine reads gray body factor data from the 

temporary restart tape and writes it to a scratch 
file and the temporary restart output file. 

CALLING SEQUENCE: CALL GBRTI (IGO, A) 

IGO - flag indicating stage reached in matrix 
inversion before restart 
A - storage array 


REFERENCE: GBINV 

FILES: NRTI, NRTO, NSCRl, NSCR2, NSCR3, NOUT 
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ROUTINE NAIffi 


GBSCFA 


DESCRIPTION: Given the set of form factors, this routine blocks 

the matrix and applies the proper factors to generate 
a positive definite matrix, and guarantees that the 
inverse exists. 

CALLING SEQUENCE: CALL GBSCFA 

REFERENCES ; GBMAIK 

FILES: Reads NFF, writes NSCRl 

NGBIR, NGBSO, NOUT, NRSO, NSCR3 


ROUTINE NAliE; I ACT 

DESCRIPTION: This routine returns the number of elements in an 

array defined in the array data block. The integer 
count of the array is stored in the zeroth cell of 
the array. 

CALLING SEQUENCE: Function lACT (lA) 

lA - First data word of array 

REFERENCES : PLDATA, RCDATA 
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ROUTIHE NAME 


INITF 




I 


DESCRIPTION: 


This routine is a series of entry points that allow 
the user to rewind and reuse data files as desired. 


CALLING SEQUENCES: CALL DIINIT 

CALL FFINIT 
CALL 

CALL AQINIT 
CALL QOINIT 


REFERENCE: ODPROG (Users Operations Data) 

FILES: NDI, NFF, NGBSO, NGBIR, NTQ, NOUT 


ROUTINE NAME: INIT28 

DESCRIPTION: This routine initializes the plot link and sets up 

the plot file. 

CALLING SEQUENCE: CALL INIT28 

REFERENCES: PLMAIN, NPMAIN, OPMAIN 
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ROUTINE NAME: INTOD 


DESCRIPTION; This routine sets up the NTITLE array used in page 

headings and is called each time the operations data 
are called after a link call. 

CALLING SEQUENCE: CALL INTOD 

REFERENCES; OPPROG (generated by preprocessor) 


ROUTINE NAME: JUPIDD 

DESCRIPTION: This routine sets up the planet parameters concerning 

the planet Jupiter. 

CALLING SEQUENCE: CALL JUPIDD 

KEY VARIABLES: PRAD - Planet radius 

SOL - Solar constant 
PALB - Planet albedo factor 
HDS - Planet darkside emissive power 
USS - Planet subsolar emissive power 
GRAV - Gravitational constant 

REFERENCES: ORBITI, ORBIT2 , DITTP 


ROUTINE NAJIE: LINE 

DESCRIPTION: This routine connects the points P]^ and P 2 defined by 

Pi = (Xi, Yi), P2 = (X2, Y2) 

CALLING SEQUENCE: CALL LINE (Xl, Yl, X2, Y2) 

XI, Yl - Coordinates of point 1 
X2, Y2 - Coordinates of point 2 

REFERENCES: NPFPLT, OPFPLT, PLGRID 
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ROUTINE NAME : LINEOP 


DESCRIPTION; This routine defines the intensity of the lines drawn 
by the plot routines. 

CALLING SEQUENCE; CALL LINEOP (DUM, INTEN) 

DUM - Dummy parameter 

INTEN - Variable defining intensity 

REFERENCES; NPSCAL, OPS CAL 


ROUTINE NAilE: LINES 

DESCRIPTION; This routine connects the arrays of X and Y with a line. 

CALLING SEQUENCE: CALL LINES (X, Y, N) 

X - Array of X 
Y - Array of Y 

N - Number of points in X and Y arrays 


REFERENCES: PLDRIV 


ROUTINE NAME: LIST 

DESCRIPTION: This routine allows the user to print and line-number 

data written on files NUSERl and NBCDOU by the FF, DI, 
RC and QO segments. 

CALLING SEQUENCE: CALL LIST (NF, lEOF) 

NF - name of file to be listed 
lEOF - number of end-of-files to be passed 
during list operation 

REFERENCE : ODPROG (Users Operations Data) 

FILES: NUSERl, NBCDOU, NOUT 
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BOUTINE NAME : MAP 


DESCRIPTION; This routine maps the plot object space into subject 
space . 

CALLING SEQUENCE: CALL NUP (XMN, YMIN, XMAX, YMAX, XMI, YMA, YMA) 

XMIN, YMIN, XMAX, YMAX - Corner points of object 

space 

XMI, XMA, YMC, YMA - New corner points of 

subject space 

REFERENCES; NPSCAL, NPINFO, OPSCAL, OPINFO, PLGRID 


ROUTINE NAME: MARSD 

DESCRIPTION: This routine defines the planet parameters concerning the 

planet Mars. 

CALLING SEQUENCE: CALL MARSD 


KEY VARIABLES: 


PRAD - Planet radius 

SOL - Solar Constant 

PALB " Planet albedo factor 

WDS - Planet darkside emissive power 

HSS - Planet subsolar emissive power 

GRAV - Gravitational constant 


REFERENCES: ORBITl, ORBIT2, DITTP 


I 


ROUTINE NME: ?1ERCUD 

DESCRIPTION: This routine sets up the planet parameters concerning 

the planet Mercury, 

CALLING SEQUENCE: CALL MERCUD 


KEY VARIABLES: 


PRAD - Planet radius 
SOL - Solar constant 
PALB - Planet albedo factor' 
IJDS — Darkside emissive power 
WSS - Subsolar emissive power 
GRAV “ Gravitational constant 


REFERENCES: ORBITl , ORBIT2, DITTP 
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ROUTINE NAME : MODAR 


DESCRIPTION: This routine is called by the user in order to change 

the value of node areas. 

CALLING SEQUENCE; See Users Manual, Appendix D. 

REFERENCES: ODPROG (Users Operations Data Block) 

ROUTINE NAME: MO DPR 

DESCRIPTION: This routine is called by the user in order to change 

the values of diffuse nodal surface properties, 

CALLING SEQUENCE: See Users Manual, Appendix D. 

REFERENCES; ODPROG (Users Operations Data) 

GHGPRP 

ROUTINE NAME : HODPRS 

DESCRIPTION: This routine is called by the user in order to change 

the values of specular surface properties. 

CALLING SEQITINCE: See Users Manual, Appendix D, 

REFERENCES; ODPROG (User Operations Data) 

CHGPRP 

PILES: NRAN (random access file) 

NOUT 
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ROUTINE NAME: 


MODSHD 


DESCRIPTION: This routine is called by the user in order to modify 

the SHADE/BSHADE flags for surfaces. 

GALLING SEQUENCE: See Users Manual, Appendix D. 

REFERENCES: ODPROG (Users Operations Data) 

FILES: NRAN, NOUT 


ROUTINE NAME : MODTR 

DESCRIPTION: This routine is called by the user in order to change 

the values of nodal surface ;ransraissivities . 

CALLING SEQUENCE; See Users Manual, Appendix D. 

REFERENCES: ODPROG (Users Operations Data) 

FILES : NRAN, NOUT 


ROUTINE NAME : MOOND 

DESCRIPTION: This routine defines the planet parameters concerning 

the moon. 


CALLING SEQUENCE: CALL MOOND 

KEY VARIABLES: TRAD - Planet radius 

SOL - Solar constant 
PALB - Albedo factor 
WDS - Darkside emissive power 
WSS - Subsolar emissive power 
GRAV - Gravitational constant 

REFERENCES: ORBITl, 0RBIT2, DITTP 
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ROUTINE NAME: KDATA 


DESCRIPTION: This is a user routine to preset the node plot optional 

parameters. 

CALLING SEQUENCE: CALL NDATA (NV, IW, SCL, ISELN, ITIT, NPHI, NPSI, 

NOMI, PHI, PS I, OMI) 

KEY VARIABLES: See users manual Appendix D. 

REFERENCES: User-defined call in the Operations Data block 


ROUTINE NAME : NDATAS 

DESCRIPTION: This routine is the short-form, user-callable routine 

to define the node plotter options the user wishes to use. 

CALLING SEQUENCE: CALL NDATAS (NV, IVU, .SCL) 

NV - Plot frame number (1 ^ NV ^ 6) 

IVU - Plot view type (X, Y, Z, 3-D, GEN) 

SU - Plot scale number 

REFERENCES: User-called in the Operations Data block. 


ROUTINE NAME : NDUPCK 

DESCRIPTION; This routine checks to see if duplicate node numbers 
are present in the node array. 

CALLING SEQUENCE: CALL NDUPCK 

REFERENCES: BUILDC 



vy 



ROUTINE NAME: 


HEFTD 


DESCRIPTION: This routine defines the planet parameters concerning 

the planet Neptune. 

CALLING SEQUENCE: CALL NEPTD 

KEY VARIABLES: PRAD - Planet radius 

SOL - Solar constand 
PALE - Planet albedo factor 
WDS - Darkside emissive power 
WSS - Planet subsolar emissive power 
GRAY - Gravitational constant 

REFERENCES: ORBITl, 0RBIT2, DITTP 


ROUTINE NAME: NODDAT 

DESCRIPTION: This routine is called by the user in order to print 

the current values of all nodal surface properties 
and their areas. 

CALLING SEQUENCE: See Users Manual, Appendix D. 

REFERENCES: ODPROG (Users Operations Data Block I 

FILES : NOUT I 
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ROUTINE NAME 


NPAXES 


DESCRIPTION: This routine defines the size and orientation of all 

characters to be small and horizonal, and draws and 
labels the CCS axes. 

CALLING SEQUENCE; CALL NPAXES 

KEY VARIABLES: VROT (3, 3) - Direction cosine matrix converting 

points in the plotter reference coordinate 
system to points in the desired view. 

REFERENCES : NPMAIN 


ROUTINE NAME; NPCONE 

DESCRIPTION; This routine* given the surface/node dimensions in the 
SCS generates the portion of the cone defined* converts 
it to the CCS, and plots the results. 

CALLING SEQUENCE; CALL NPCONE 

REFERENCES ; NPMAIN 


ROUTINE NAliE : NPCONV 

DESCRIPTION: This routine converts the plot control data from common 

and checks for errors. If no data are defined in common, 
it sets up the proper default parameters. 

CALLING SEQUENCE; CALL NPCONV (NV, NNP , NEND) 

NV - Plot view frame number 

NNP - Number of nodes to be selectively plotted 
NEND - End flag determining if more plots are needed 

REFERENCES; NPMAIN 


ROUTINE NAME: NPCYLQ 

DESCRIPTION: This routine, given the surface /node dimensions in the 

SCS, determines the proper plot calls to draw a cylinder. 

CALLING SEQUENCE: CALL NPCYLO 

REFERENCES; NPMAIN 
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ROUTINE NAME: NPPISC 


\ 


DESCRIPTION: This routine, given the surface/node dimensions in the 

SCS, determines the proper plot calls to draw a disk. 

CALLING SEQUENCE: CALL NPDISC 

REFERENCES: NPMAIN 


ROUTINE NAME: NPPOTL 

DESCRIPTION: This routine connects PI and P2 with a dotted line. 

CALLING SEQUENCE: CALL NPDOTL (XI, Yl, Zl, X2, Y2, Z2) 

XI, Yl, Zl - Coordinates of PI 
X2, Y2, Z2 - Coordinates of P2 

REFERENCES; NPCONE, NPCYLO, NPPARA, NPSPHE 


ROUTINE NAI-IE; HPFPLT 

DESCRIPTION: This routine directs the actions of the pen/beam 

plotter. I is either even or odd, or negative or 
positive. 

CALLING SEQUENCE: CALL NPFPLT (I, X, Y) 

I - Plotting designator 

- Even, [Draw to (X, Y) ] 

- Odd, [Position at (X, Y) ] 

- Positive, (generates new origin) 

- Negative, (keeps same origin) 

X, Y - Coordinates of new point 

REFERENCES: NPAXES, NPCONE, NPCYLO, NPDOTL, NPTPLT, NPDISC, NPPARA, 

NPRECT, NPSPHE, NPTEAP 
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ROUTINE NAME: 


NPINFO 


DESCRIPTION; This routine labels the plot. 

CALLING SEQUENCE; CALL NPINFO (NV, KEND, KCC) 

NV - View number 
KEND - = 0 Good plot 

^ 0 Error in plot 
KCC “ Internal view number 

REFERENCES ; NPMAIN 
ROUTINE NAiffi: NPMAIN 

DESCRIPTION: This routine is the main driving logic to control the 

node plot link. 

CALLING SEQUENCE; CALL NPMAIN 
REFERENCES: NPPROG (preprocessor-generated) 


ROUTINE NAME; NPPARA 

DESCRIPTION: This routine, given the node dimensions in the SCS, 

generates the plot calls necessary to draw a paraboloid. 

CALLING SEQUENCE; CALL NPPARA 

REFERENCES : NPMAIN 


ROUTINE NAME: NPPROG 

DESCRIPTION; This routine is generated by the preprocessor. It defines 
the necessary labeled common blocks and calls into the 
main driving logic for the node plotter. 

CALLING SEQUENCE; CALL NPMAIN 

REFERENCES; TRASYS (root segment) 
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ROUTINE NAME; NPRECT 

DESCRIPTION: This routine j given ‘‘he nodal dimensions, generates the 

plot calls to d . jw a rectangle, 

CALLING SEQUENCE; CALL NPRECT 

REFERENCES : NPMAIN 


ROUTINE NAME: NPRIDA 

DESCRIPTION: This routine is unique to the Univac version, and 

defines a system- labeled common block with a user name 
and location for labeling the plots. An initialization 
routine is called from this routine to initialize the 
plot. 

CALLING SEQUENCE; CALL NPRIDA (lARRAY) 

lARRAY (1) - Name (word I) 
lARBAY (2) - Name (word 2) 
lARRAY (3) - User mail number 
lARRAY (4) - Extension 
lARRAY (5) “ Blank 
lARRAY (6) - Project 

REFERENCES; User-callable routine from the Operations Data block 
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ROUTINE NAilE : NPRNT 


DESCRIPTION: This routine prints a sununary of node information as the 

configuration is generated by routines BUILDC/ADD. 

CALLING SEQUENCE: CALL NPSNT (NEMASS » NAME) 

NRMASS “ Array of node data 

NAME - Block coordinate system name 


REFERENCES : 


BUILDC 


ROUTINE NAME: NPROTA 

DESCRIPTION: This routine defines the transformation of the direction cosine 

matrix, given the desired view. 

CALLING SEQUENCE; CALL NPROTA 

KEY VARIABLES: KG Desired view number 

PH, PS, OM “ Angles necessary to arrive at desired view 

REFERENCES : NFMAIN 


ROUTINE NAME: NFS CAL 

DESCRIPTION: This routine determines the grid subject space, maps 

into the new space, and locates the pen/beam at the new 
origin. 

CALLING SEQUENCE: GALL NPSCAL (XS, YS, X, Y) 

XS - Scale factor 
YS “ Scale factor 
X - 0.0 
Y - 0.0 

REFERENCES: NPAXES 


ROUTINE NAME: NPSPHE 

DESCRIPTION: This routine, given the nodal dimensions, generates the 

plot calls necessary to draw a sphere or segment of a 
sphere. 

GALLING SEQUENCE: CALL NPSPHE 

REFERENCES : NPMAXN 

REV 1 
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ROUTINE NAME: 


NPTPLT 


DESCRIPTION; This routine transforms a point in 3-D space to X, Y coordinates 
in the subject space, checks to determine if they are in the 
allowable range, and moves the pen/beam to the coordinates of 
that point. 

CALLING SEQUENCE: CALL NPTPLT (X, XP, YP, X3, Y3, Z3) 

I - Pen-up, pen-down flag 

^XP, YP - Coordinates of transformed point 
X3, Y3, Z3 - Coordinates of input point 

REFERENCES; NPCONE, NPCYLO, NPDISC, NPPARA, NPRECT, NPSPHE, NPTRA:? NPMAIN | 


ROUTINE NAME: NPTRAP 

DESCRIPTION; This routine, given the nodal dimensions in the SCS generates 
the plot calls necessary to draw a trapezoid. 

CALLING SEQUENCE: CALL NPTRAP 

REFERENCE S ; NPMAIN 


ROUTINE NAME: NUMBER 

DESCRIPTION; This routine, given a value and a format, generates a number 
on the plot frame. 

CALLING SEQUENCE; CALL NUMBER (X, F) 

X - Floating or integer number to be output 
F - Output format 

REFERENCES: NPINFO, OPINFO, PLDRIV, PLGRID 
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ROUTINE NAME: 


ODATA 


DESCRIPTION: This routine is user-called to define orbit data plotter 

options. 

CALLING SEQUENCE: CALL ODATA (KV, SCL, SCLR, RPLN, TRUE, TIMER, ISELN, 

ITIT, NPHF, NPSI, NOMI, OMI, PSI, PHI) 

KEY VARIABLES: See Appendix D of users manual. 

REFEKENCES: User-called in the Operations Data block 


ROUTINE NAME: QPAXES 

DESCRIPTION: This routine defines the size and orientation of all 

characters to be small and horizonal, and draws and labels 
the CCS axis and sun line. 

CALLING SEQUENCE: CALL OPAXES 

REFERENCES: OPMAIN, OPPLAN 


ROUTINE NAME: OPCOMB 

DESCRIPTION; This routine does matrix multiplication of the A and B input 
matrices and returns the results in matrix C. 

galling SEQUENCE: CALI. OPCOMB (A, B, C) 

A - Input 3x3 matrix 
B - Input 3x3 matrix 
C - Resultant matrix (C = A*B) 

REFERENCES; OPLOC, 0PL0C2 


ROUTINE NAME : OF CONE 

DESCRIPTION: This routine, given the surface dimensions in the SCS, 

generates the portion of the cone defined, converts it to 
the CCS, and plots the results. 

CALLING SEQUENCE: CALL OPCONE 

REFERENCES ; OPMAIN 
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ROUTINE NAME; 


QPCONV 


DESCRIPTION: This routine converts the plot control data from common and 

checks for errors. If no data are defined in common, it sets 
the proper default parameters . 

CALLING SEQUENCE; CALL OPCONV (NV, NNP, KEND) 

NV - Plot view frame numher 

NNP “ Number of nodes to be selectively plotted 
KEND - End flag determining if more plots are needed 

REFERENCES : OPMAIN 


ROUTINE NAME; QPCYLQ 

DESCRIPTION: This routine, given the surface/nodal dimensions in the SCS, 

determines the proper plot calls to draw a cylinder. 

CALLING SEQUENCE: CALL OPCYLO 

REFERENCES; OPMAIN 


ROUTINE NAME: OP DISC 

DESCRIPTION: This routine, given the surface/nodal dimensions in the SCS- 

determines the proper plot calls to draw a disk. 

CALLING SEQUENCE: CALL OPDISC 

REFERENCES : OPMAIN 

ROUTINE NAME; OPDOTL 

DESCRIPTION; This routine connects PI and P2 with a dotted line. 

CALLING SEQUENCE: GALL OPDOTL (XI, Yl, Zl, X2, Y2, Z2) 

XI, Yl, Zl - Coordinates of PI 
X2, Y2, Z2 - Coordinates of P2 

REEERENCES: OPCONE, OPCYLO, OPPARA, OPSPHE 
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ROUTIHE NAME: 


OPEDOT 


DESCRIPTION: This routine generates the dotted lines used in the planet- and 

planet-shadow generation routines and connects points PI and P2. 

CALLING SEQUENCE: CALL OPEDOT (XI, Yl, Zl, X2, Y2, 22, A) 

XI, Yl, 21 - Coordinates of PI 
X2, Y2, 22 - Coordinates of P2 

A - Factor determining the length of the line increment 
REFERENCES: OPPLAN, OPSHAD 


ROUTINE NAME: OPFPLT 

DESCRIPTION: This routine directs the actions of the plot pen/beam, I may 

be even or odd, negative or positive. The value of I directs 
the plotter to the coordinates X, Y with the pen up or doxm. 

CALLING SEQUENCE: CALL OPFPLT (I, X, Y) 

I - Plot designator 

= Even, draw to (X, Y) 

= Odd position 

= Positive (generates new origin) 

= Negative (keep satne origin) 

REFERENCES: OPAXES, OPCONE , OPDOTL, OPTPLT, OPDISC, OPPASA, OPEECT, OPEDOT, 

OPSPHE, OPTRA? , OPVCS, OPCYLO 


ROUTINE NAME: OPINF O 

DESCRIPTION: This routine labels all plot alphabetic data on the plot frame. 

CALLING SEQUENCE: CALL OPINFO (NV, KEND) 

NV - View frame number 
KEND - Error flag 

REFERENCES ; OPMAIN 


ROUTINE NAME: OPLOC 

DESCRIPTION; This routine, given the oribital parameters, determines the proper 
direction cosine matrix to give the proper orientation, 

GALLING SEQUENCE t GALL OPLOC 

REFERENCES : OPMAIN 
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ROUTINE NAME; 0P LQC2 

DESCRIPTION: This routine locates the vehicle relative to the sun 

and planet according to parameters defined in subrou- 
tine DIDT2. 

CALLING SEQUENCE: CALL 0PL0C2 

REFERENCES OPLOG 


ROUTINE NAME: OPMAIN 

DESCRIPTION; This is the main driving logic to control the orbit 
plot link. 

CALLING SEQUENCE: CALL OPMAIN 

REFERENCES: OPPROG (preprocessor-generated) 


ROUTINE NAME; OPMAX 


DESCRIPTION: This routine determines the maximum radius needed to 

enclose any surface. This value is then used to 
determine the proper scale factor. 


CALLING SEQUENCE: SCL - OPMAX (NNP) 

NNP - Number of surfaces to be selectively 
plotted 


REFERENCES : OPMAIN, OPCONV 


ROUTINE NAME: OPPARA 

DESCRIPTION: This routine, given the surface dimensions in the SCS, 

generates the plot calls necessary to draw a paraboloid. 

CALLING SEQUENCE- CALL OPPARA 

REFERENCES : OPMAIN 
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ROUTINE HAHE: OPPLAN 

DESCRIPTION: This routine, given the desired planet radius for plotting 

(in inches), generates the logic necessary to draw the planet, 

CALLING SEQUENCE: GALL OPPLAN 

REFERENCES: OPMAIN 


ROUTINE NAME: OPPRNT 

DESCRIPTION: This routine generates a summary table of orbital parameters 

on the output file. 

CALLING SEQUENCE: GALL OPPRNT 

REFERENCES: OPMAIN 


ROUTINE NAME: OPPROG 

DESCRIPTION: This routine is generated by the preprocessor. It defines the 

necessary labeled common blocks and the calls into the main 
driving logic. 

CALLING SEQUENCE: GALL OPPROG 

REFERENCES: TRASYS (root segment) 


ROUTINE NAME: OPRECT 

DESCRIPTION: This routine, given the surface dimensions, generates the plot 

calls to draw a rectangle. 

CALLING SEQUENCE: GALL OPRECT 


REFERENCES: OillAlN 


ROUTINE NAME: OPROTA 


DESCRIPTION: This routine defines the transformed direction cosine 

matrix, given the desired view. 

CALLING SEQUENCE: CALL OPROTA 

KEY VARIABLES: KC - Desired view number 

PH, PS, OM - Angles necessary to arrive at desired view 

REFERENCES : OPMAIN 


ROUTINE NAME: OPSCAL 

DESCRIPTION: This routine determines the grid subject space, maps into 

the new space, and locates the pen/beam at the new origin. 

CALLING SEQUENCE: CALL OPSCAL (XS, YS, X, Y) 

XS Scale factor 
YS - Scale factor 
X = 0 
Y = 0 

REFERENCES: OPAXES, OPMAIN 


ROUTINE NAME : OPSHAD 

DESCRIPTION: This routine draws the planet shadow. 

CALLING SEQUENCE: CALL OPSHAD 

REFERENCES : OPtiAIN 
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ROUTINE NAME ; OPSPHE 


DESCRIPTION: This routine, given the surface dimensions, generates the 

plot calls necessary to draw a sphere or segment of a 
sphere . 

CALLING SEQUENCE: CALL OPSPHE 

REFERENCES: OPMAIN 


ROUTINE NAME: OPTIM E 

DESCRIPTION: This routine, given the orbit eccentricity, present orbital time 

and orbital period, computes the time anomaly 

CALLING SEQUENCE: CALL OPTIME (TIME, PER, ECC, TRU) 



TIME - 
PER “ 
ECC - 
TRU - 

Input orbital time 
Orbit period 
Orbit eccentricity 
Output true anomaly 

REFERENCES: 

OPLOC 


ROUTINE NAME: 

: OPTPLT 



DESCRIPTION: This routine transforms a point in 3~D space to X, Y coordinates 

in the subject space, checks to determine if the points are in 
the allowable range, and moves the pen/beam to the coordinates 
of that point. 

CALLING SEQUENCE: CALL OPTPLT (I, XP, YP, X3 , Y3, Z3) 

I - Pen-up, pen-down flag 

XP, YP - Coordinates of transformed point 

X3, Y3, Z3 - Coordinates of input point 

REFERENCES; OPGONE, OPCYLO, OPDISC, OPPARA, OPSPHE, OPTRAP, OPRECT 


ROUTINE NAME: OPTRAP 

DESCRIPTION: This routine, given the surface dimensions in the SCS, generates 

the plot calls necessary to draw a trapezoid. 

CALLING SEQUENCE: CALL OPTRAP 

REFERENCES: OPMAIN 
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ROUTINE NAME: OPTRNP 


DESCRIPTION; This routine, given a 3 x 3 matrix, computes the transpose. 
CALLING SEQUENCE: CALL OPTRNP (A, B) 


A - 3 X 3 matrix to be transposed 
B - 3 X 3 resultant matrix 


REFERENCES; 


OPDOTL, OPTPLT, 


OPLOC OPLOC2 


ROUTINE NAME; 0PTRS3 

DESCRIPTION; This routine changes a translation vector from one coordinate 
system to the corresponding translation vector in the new 
system. 

CALLING SEQUENCE: CALL 0PTRS3 (X, Y, Z, A, B, C, EX, RY, RZ, TRAN) 

X, Y, Z " New translation vector 
A, B, C - Translation vector in the old system 

RX, RY, RZ ~ Translation vector in the new system 
TRAN “ Direction cosines relating the old system to 

the new system 

REFERENCES: OPAXES, OPDOTL, OPTPLT, OPEDOT, OPVCS 

ROUTINE NAME; OFUNIT 

DESCRIPTION: This routine generates an identity matrix. 

CALLING SEQUENCE: CALL OPUNIT (A) 

A - Output unit matrix 
REFERENCES : OPLOC, 0PL0C2 

ROUTINE NAME: OPVCS 

DESCRIPTION; This routine labels the vehicle axes on the plot. 

CALLING SEQUENCE; CALL OPVCS 
REFERENCES : OPMAIN 
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ROUTINE NAME: ORBITl 


DESCRIPTION: This routine can be called by the user in the Operations Data 

block to define the orbit. 

CALLING SEQUENCE: CALL ORBITl (PLANAM, ALANI, APEI, OIN, TIME, HPI, HAl, 

SBA, SDE, STA, STD) 

KEY VARIABLES: See users manual, Appendix D 

REFERENCES: ODPROG 

ROUTINE NAME: OKBIT2 

DESCRIPTION: This is a user callable routine to define the desired orbit. 
calling SEQUENCE: CALL 0RBIT2 (PLANAM, GIG, BET, GIGS, BETS, TIME, HPI, HAI) 
KEY VARIABLES: See users manual. Appendix D 

REFERENCES : ODPROG 

ROUTINE NAME: ORIENT 

DESCRIPTION: This routine allows the user to define the vehicle orien- 

tation , 

CALLING SEQUENCE: CALL ORIENT (TYPE, IROTX, IROTY, IROTE, ROTX, ROTf, ROTZ) 

KEY VARIABLES; See users manual. Appendix D 
REFERENCES : ODPROG 

ROUTINE NAME: PAGE 

DESCRIPTION; This routine is called prior to every written statement to 

the output file to count the line printed and take care of all 
paging. 

CALLING SEQUENCE: II = PAGE (I) 

II = 0 New page was written 
1 No new page 

I - Number of lines to be printed 
REFERENCES; All output generating routines. 
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ROUTINE NA^ffi : PLCMBH 

DESCRIPTION; This routine combines incident and absorbed flux 

data according to correspondence data prior to data 
plotting. 

GALLING SEQUENCE:* GALL PLCMBN (SOL, ALB, PLA, NDCOM, NNDC, NUNIT) 

SOL ~ Solar flux/rate array 

ALB - Albedo flux/ rate array 

PLA - Planetary flux/rate array 

NOCOM - Combine/no combine flag 

NNDG “ Number of nodes after combining 

NUNIT - flux/ rate data storage file name (NDI, NTQ) 

REFEPLINCES : FNDFLP 

FILES: NTQ, NOUT, NGBIRR 


ROUTINE NAME: PLCFIT 

DESCRIPTION: This routine does .smooth-curve fitting. Given two arrays 

LO words long (X, Y) , it generates NO points bett^een each 
set of X, Y points and stores the points in arrays U and V. 

GALLING SEQUENCE: CALL PLCFIT (LO, X, Y, NO, U, V) 

LO - Number of input points 

X - Independent variable array 

Y - Dependent variable array 

NO - Number of divisions between each set of points 

U, V “ Output arrays 

REFERENCES : PLDCON 


ROUTINE NAME; PLDATA 

DESCRIPTION: This routine can be called by the user in the Operations Data 

block to define PLOT link options. 

calling SEQUENCE: CALL PLDAtA (IP, INS, IS, CRVF, TLX, TLY, Tl , T2 , 

XMPF, YMPF) 

KEY VARIABLES; See users manual, Appendix D 
REFERENCES; Operations Data block 
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ROUTINE NAME: 


PLDGON 


DESCRIPTION: This routine checks for discontinuities in the output plot 

and calls for curve fitting between all discontinuities. 

CALLING SEQUENCE: CALL PLDGON (NTIME, TIME, PARRAY, NODIV, NTOTj TIMP, 

PLOTP) 

NTIME - Number of orbit points 

TIME “ Array of times 

PARRAY ” Array of dependent data 

NODIV “ Number of divisions between each point 

NTDT “ Number of output points 

TIMP - Array of curve-fit points (independent) 

PLOTP - Array of curve-fit points (dependent) 

REFERENCES: PLDRiV 


ROUTINE name : PLDRIV 

DESCRIPTION: This routine is the main driving routine in the plot link. 

It decades the type of plot and calls the proper routines. 

CALLING SEQUENCE: CALL PLDRIV (RINDEP, S, DEPEND, RINDO, DEPO, IS, NDIV, 

NINDV) 

RINDEP - 
S 

DEPEND - 
RINDO - 
DEPO 
IS 

NDIV - 
NINDV - 

REFERENCES : PLMAIN 


ROUTINE NAME: PLGRID 

DESCRIPTION: This routine draws the plot grid and labels the complete frame. 

CALLING SEQUENCE: CALL PLGRID (mEN, XMAX, YMIN, YMAX, NODEN) 

XMIN, XMAX, YMIN, YMAX - Minimum and .aximum dimensions of 

plot frame 

NODEN - Node number of current frame 
REFERENCES : PLDRIV 


Independent variable array 
Temporary working array 
Dependent variable array 
Independent variable array computed 
Dependent variable array output 
Temporary array (same as S) 

Number of divisions beti^een curve-fit points 
Number of independent variable points allowed 
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EQUTINE NAME: ELLOAD 



DESCRIPTION: This routine determines the type of plot and finds the 

plot data if absorbed or incident flux data are to be 
plotted. 

galling SEQUENCE: CALL PLLOAD (DATA, TIME, ITIME, NTIME, NINDV) 

DATA - Working temporary array 

TIME, ITIME - Doublet array containing time and 
step number 

NTIliE " Number of time points 

NINDV - Number of independent-variable points 

REFERENCES : PLMAIN 


ROUTINE NAME: PLMAIH 

DESCRIPTION: This is the main driving routine in the plot segment. Its 

primary functions are to determine the size of the array and 
allocate array space. 

CALLING SEQUENCE: CALL PLMAIN 

REFERENCES: PLPROG (preprocessor-generated) 


ROUTINE NAME; PLOUT 

DESCRIPTION: This routine applies the dependent variable multiplier and 

converts fluxes to rates. 

CALLING SEQUENCE: CALL PLOUT (NV, Al, A2, A3, A4, ITYPE) 

NV - Number of arrays 

Al, A2, A3, A4 - Data arrays 
ITYPE - Flag indicating flux or rate 

^ 4HFLUX heat flux 
= 4HRATE heat rate 

REFERENCES : FNDFLP 


ROUTINE NAME: PLPROG 

DESCRIPTION : This routine is generated by the preprocessor and calls in the main 

processor routines. 

CALLING SEQUENCE: CALL PLPROG 

REFERENCES: TRASYS (root segment) 
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ROUTINE NAME: RLSCL 


DESCRIPTION: This routines given the maximum and minimum values, determines 

a scale using convenient units . 

CALLING SEQUENCE: CALL PLSCL (BMAX, BMIN, NSQ, AD, AL, S) 

BMAX, BMIN - Maximum and minimum values to be plotted 

NSQ “ Increment desired 

AU, AL - New upper and lower values 

S - Scale factor 

REFERENCES: PLGRID 


ROUTINE NAME: PLSYMB 

DESCRIPTION: This routine, given an array of Hollerith data, outputs the 

data to the plot frame and supplies the proper character 
terndnator . 

calling SEQUENCE: CALL PLS-XMB (ARRAY, N) 

ARRAY - Array of data to be printed 
N - Number of words in array 

REFERENCES; PLGRID 
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ROUTIIJE NAME 


PLUTOD 


DESCRIPTION: This routine defines the' planet parameters concerning 

the planet Pluto. 


CALLING SEQUENCE : CALL 

PLUTOD 


KEY VARIABLES : PRAD 

- Planet radius 


SOL 

- Solar constant 


PALB 

- Albedo factor 


WDS 

- Darks id a emissive 

power 

WSS 

- Subsolar emissive 

power 

GRAV 

- Gravitational constant 


REFERENCES: ORBITI, 0KBIT2, DITTP 


ROUTINE NAME : PRDUMP 

DESCRIPTION: This routine writes the node, area and property 

arrays to the RSO file. 

CALLING SEQUENCE: CALL PRDUMP (LABEL) 

LABEL - Flag indicating whether properties 

written out were defined per surface 
data block, as modified by "MOD" 
routines or are combined over corres- 
pondence data. 


REFERENCES: ODPRGG, CMMAIN 

FILES: NRSO, NOUT 


REV 1 
III -103 


ROUTIHE NAME : QOAVGR 

DESCRIPTION: This routine computes the integrated average of the 

absorbed heat rates that were input. 

CALLING SEQUENCE: CALL QOAVGR (NTIME, TIME, NDNIT) 

NTIME - Number of time points 

TIME ” Doublet array of time and step number 

NUNIT •“ Unit, number sorted fluxes heat rates are stored on 

REFERENCES: QOMAIN 


ROUTINE NAME : QOAVGS 

DESCRIPTION: This routine, given the final data, outputs the data according 

to the user requirements. 

CALLING SEQUENCE: CALL QOAVGS (NODNO, QAV, AREAT) 

NODNO “ Node number 
QAV - Averaged Q value 
AREAT - Area of node 


REFERENCES : QOSAVE 


ROUTINE NAME: QQCMBH 

DESCRIPTION: This routine generates the combining data read from the 

correspondence data. 

GALLING SEQUENCE: CALL QOCMBN (IGOMBL, IPIRSL) 

ICOMBL - Number of points in Combine Array 1 
IFIRSL “ Number of points in Combine Array 2 


REFERENCES: QOMAIN 
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ROUTINE NAME ; QOGOMB 


DESCRIPTION: Based on the combine arrays, this routine combines and stores 

the data. A final call to this routine causes the combined 
tables to be output 

CALLING SEQUENCE: CALL QOCOMB (DATAl, DATA2, WDIMI, ND1M2, NDIM, IFIRSL, 

ICOMBL, DATAs IRON, KTAB) 

DATAl, DATA2 - Working storage blocks 

NDIMI, ND1M2 ~ Length of data' .blocks 1 and 2 

NDIM - Number of nodes 

lEIRSL, IGOMBL - Length of combining arrays 

DATA “ Working storage array 

IROW “ Sequence number of current row 

KTAB - Table number 

REFERENCES : QOSAVE 


ROUTINE NAME : QQDATA 

DESCRIPTION: This routine is called by the user to define particular 

parameters for the absorbed heat rate output. 

CALLING SEQUENCE; CALL QODATA (lARRAY, NTMARY, QOTP, SUN, AMPF, FMPF, 

TMPF, TYPE, KCOR) 

KEY VARIABLES; See users manual, Appendix D for definition of variables 
REFERENCES : Operations Data block 


ROUTINE NAME; QOFLIP 

DESCRIPTION: This routine, given a matrix (NTIME x ICOL) stored on disk 
NSl, converts this to a matrix (ICOL x NTIME) on unit NS2 . 

CALLING SEQUENCE: CALL QOFLIP (DATA, NTIME, ICOL, NSl, NS2) 


DATA 

NTIME 

ICOL 

NSl, NS2 


Working storage area 
Number of rows 
Number of columns 
Input, output unit numbers 


REFERENCES : 


QQMAIN 



ROUTINE NAME: 


OOHEAD 


DESCRIPTION: This routine outputs default and user control parameters 

to the output file. 

CALLING SEQUENCE : CALL QOHEAD 

REFERENCES : QOMAIN 

ROUTINE NAZiE: QQMAIN 

DESCRIPTION: This routine is the main driving logic for the QOCAL link. 

CALLING SEQUENCE: GALL QOMAIN 

REFERENCES: QOPROG 


ROUTINE NAME: QOPROG 

DESCRIPTION: This routine is defined by the preprocessor. It defines 

all needed labeled common blocks and calls in the driving 
logic for the calculation of absorbed heat rates. 

KEY VARIABLES : CALL QOPROG 

REFERENCES: TRASYS (root segment) 


ROUTINE NAME: QQSAVE 

DESCRIPTION: This routine reads the uncombined data points andj from the 

input options directs the combining and output options. 

CALLING SEQUENCE: CALL QOSAVE (I GOMEL, IFIRSL, NTIME, TIME) 

ICOMBL, IFIRSL - Length of the combining arrays 

NTIME - Number of time points 

TIME ” Double array of step numbers and times 

REFERENCES: QOMAIN 
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ROUTINE MME: QOSBCD / 

DESCRIPTION: This routine is called to output /he subroutine Call cards. 

CALLING SEQUENCE-: CALL QOSBCD (PER, ITIME, lOCAB, AREAT, NOMO) 

PER - Orbit period / 

ITIME " Time array numb^ 

KTAB - Table nuniber rexerence number 
AREAT - Area / 

NODNO - Node number ’ 

REFERENCES : QOSAVE ^ 



ROUTINE NAME: QOTABS 

' ' ' / 

DESCRIPTION: This routine is calle/ to generate the final output tables. 

/' 

CALLING SEQUENCE: CALL QOTABS (KTAB, DATA, NTIME, ARE) 


KTAB - Output array number 
DATA ~ Arrays of Q data 
NTIME - Number of data points 
ARE “ Area of none 


REFERENCES : QOSAVE, QOCOMB 


ROUTINE NAME: QOTMES E 

DESCRIPTION: This routine is called to output the time array in final form. 

CALLING SEQUENCE: GALL QOTMES (itmE, TIME, NTIME) 

XTME - Time array number 
TIME - Array of times 
NTIME - Number of time points 

REFERENCES : QOSAVE 
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ROUTINE NAME; RBAREA 


DESCRIPTION; This is a function routine that calculates elemental areas 
for, unevenly distributed elements of unequal size. 

CALLING SEQUENCE: RBAREA (ILPj ALREj BMIN, BMAX, GMN, GMAX) 

ILP - Surface type 
ALPH - \ 

BMIN -/ 

BMAX - V Surface dimensions of current element 
GMEN “( 

GMAX -) 

REFERENCES ; RBEXPN, RBXPNI 


ROUTINE NAME; RBCAL 

DESCRIPTION; This routine contains the driving logic to determine the 

number and distribution of elements. on node IN, as well as on 
the image of node IN as seen in specular surface MIRROR. It 
also calcxilates the image factor from node IN to node JN. 

GALLING SEQUENCE: CALL REGAL (MIRROR) 

MIRROR - Surface sequence number of current specular 
surface 

KEY VARIABLES; IN, JN - Sequence numbers of current node pair 

FE, FA - IR and solar form factors from node IN to image of 
node JN in specular surface MIRROR 
V?E, WA - IR and solar shadowing factors 

RBVALI, RBVALS ~ IR and solar image factors between node IN 
and node JN 

KN “ Sequence number of current specular node 
SREFLI, SREFLS - Arrays of IR and solar specular reflectance 
values 

REFERENCES : RBPRE 


FILES: NOUT - System output file 



ROUTINE NAME: 


RBCALI 


DESCRIPTION: This routine is analogous to REGAL. It is called 

when the I - J order of the nodes was reversed due 
to area differences. 

CALLING SEQUENCE: See REGAL 

REFERENCES : See REGAL 

FILES: See REGAL 


ROUTINE NAME : RE DATA 

DESCRIPTION: This routine is called by the user in order to define 

parameters necessary to execute the REGAL segment. 

CALLING SEQUENCE: CALL RBDATA (NFIGFF, FFACC, FFACCS, FFRATL, 

FFPRT) 

NFIGFF - Configuration name for form factor 
access 

FFACC ^ 

FFACCS See subroutine FFDATA, Users Manual, 
FFRATL ^ Appendix D 
FFPRT J 

REFERENCES: ODPROG (Users Operations Data) 
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ROUTINE NAME: 


RBELEM 


LESCRIPTION: Glvan the number of elements in each direction, this 

routine calculates the elemental position and area 
vectors - 


CALLING SEQUENCE: CALL RBELEM (NB , NG, ILP , IC, DATA, POS, ARA, 

TRAN, RX, RY, RZ) 


NB, NG 

ILP 

IC 


DATA 
POS 
ARA 
TRAN 
RX, RY, 


- Number of elements in the beta and 
gamma directions 

- Surface type 

- Counter for total number of elements 
on node pair 

- Array of node dimensions 

- Array of elemental position vectors 

“ Array of elemental area vectors 

- Matrix of direction cosines 

RZ - Components of SCS origin position vector 
in the ICS, BCS , or CCS 


REFERENCES: RBCAL, RBEXPN , RBXPNI 


ROUTINE NAME: RBELSL 

DESCRIPTION: Given the total number of elements required on a node, 

this routine makes them as square as possible and deter- 
mines the number of elements in each direction. 

CALLING SEQUENCE; CALL RBELSL (NB, NG, ILP, DATA, NTOT) 

NB, NG - Number of elements in the beta and gamma 
directions 

ILP - Surface type 

DATA - Array of node dimensions 

NTOT - Number of elements required on a node 


REFERENCES: RBCAL, RBEXPN, RBCALI, RBXPNI 


ROUTINE NAME: 


EBEND 


DESCRIPTION: This routine provides for user intervention after the 

image factor calculation. It normally calls for a time 
accounting for the problem and end-files the image factor 
file on NFF. 

CALLING SEQUENCE: CALL RBEND ' 

REFERENCES : REMAIN 

FILES: NFF - Disk file for storing image factors. 


ROUTINE NAME: RBESUM 

DESCRIPTION: This routine provides an accounting, in CP seconds, of 

the time required to calculate image factors. 

CALLING SEQUENCE; CALL RBESUM 

REFERENCES ; RBEND 

FILES; NOUT - System output file 
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ROUTINE NAME 


RBEXPN 


DESCRIPTION; This routine expands nodes into subnodes, determines 

the number and distribution of elements on each subnode, 
and calculates more accurate image factors than would 
be possible on a nodal basis. 


CALLING SEQUENCE; CALL RBEXPN (POSI, ARAI, FE, FA, RATI, RATJ, JFLAG, 

JST, NEI, NEJ, NSUR, lADDR, MIRROR) 


POSI 
ARAI 
FE, FA 
RATI, RATJ - 

JFLAG 

JST 

NEI , NEJ 
NSUR 
lADDR 
MIRROR 


Array of elemental position vectors 

Array of elemental area vectors 

Resultant IR and solar form factors 

Ratio of maximum number of elements to the no. 

of elements indicated by aritimnetic averaging 

Flag indicating that maximum number of 

elements was exceeded 

Number of possible shadowers 

Number of elements on nodes I and J 

Actual number of shadowing surfaces 

Maximum time return sequence (not used) 

Surface sequence number of current 

specular surface 


REFERENCES: RBCAL 


ROUTINE NAME: RBGTST 

DESCRIPTION: This is an integer function routine that tests a given 

value, GT, to determine if it falls in the range GN£GT_<GX. 
If it does, the function value is 0 (zero)i if not, the 
function value is 1 (one) . 

CALLING SEQUENCE: RBGTST (GN, GX, GT) 

GN - Lower bound 
GX “ Upper bound 
GT - Value to be tested 

REFERENCES: RBSHD 
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ROUTINE NAME! SBIMAG 


DESCRIPTION: This routine images surfaces for shadowing purposes 

and stores the results on NSCR2 for later use. 

CALLING SEQUENCE: CALL RBIMAG (NST) 

NST - Total number of shadowing surfaces 

KEY VARIABLES: JS - Surface sequence number of current 

specular surface 

NSURFS “ NST plus number of images in specular 

surface JS 

IPS, IKS, PR, 

DSTR, DIMS, 

PSH, TSTR - Surface description parameters 
REFERENCES : REMAIN 

FILES: NSCR2 - Scratch file 


ROUTINE NAME : REMAIN 

DESCRIPTION; This routine contains the main driving logic for calcu- 
lating image factors. 

CALLING SEQUENCE: CALL REMAIN 

KEY VARIABLES: IN - Sequence number of "viewer" node 

JN - Sequence number of node whose image in KN 

is viewed by IN 

KN - Sequence number of specular node 

ISPEC - Array of specular surface sequence numbers 

REFERENCES: RBPROG (preprocessor-developed) 

FILES: NRAN - Random access file 

NSCR2 - Scratch file 
NOUT, NRTI, NRTO 
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ROUTINE NAME: RBNIMG 


DESCRIPTION: This routine images the current node JN in the current 

specular surface. 

CALLING SEQUENCE: CALL RBNIMG (JM, NOIM) 

JM - Surface sequence number of current specular 
surface 

NOIM ~ Flag to indicate if an image of the current 
node exists in JM 
=0, an image exists 
= 1, no image 


KEY VARIABLES: POSNJ, RXJ i 

RYJ, RZJ, DATAJ,> Nodal description parameters 
TRANJ ; 

REFERENCES : RBtlAIN 

FILES: NRAN “ Random access file 


ROUTINE NAME: RBOUT 

■ 

DESCRIPTION; This is a user routine that can be replaced to change the 
form of output for image factors. The normal calls are to 
RBPNCH (punches image factors) and RBPRNT (prints image 
factors) . 

CALLING SEQUENCE: CALL RBOUT 

REFERENCES; RBMAIN, RBRSI 


V 


) 
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ROUTINE NAME; 


RBPOSI 


DESCRIPTION; This routine transposes a vector in the ICS, BCS, or CCS 
into the SCS of the current specular surface, negates the 
Z component, and transforms the vector back into the ICS, 
BCS, or CCS. 

CALLING SEQUENCE: CALL RBPOSI (X, Y, Z, TRAN) 

X, Y, Z - Components of vector in the ICS, BCS, or CCS 
TRAN - Matrix of direction cosines 

KEY VARIABLES: A, B, C - Components of vector in the SCS 

REFERENCES: RBIMAG, EBNIMG 


ROUTINE NAME; RBPRE 

DESCRIPTION: This routine provides for user intervention prior to 

the calculation of an image factor. The normal call is 
to RBCAL, which calculates the image factor. 

CALLING SEQUENCE: CALL RBPRE (MIRROR) 

MIRROR - Surface sequence number of current specular 
surface 


REFERENCES: REMAIN 


ROUTINE NAME: RBPRNT 

DESCRIPTION: This routine prints the image factor between node IN and 

node JN. The user can override this routine to change 
the output format if he desires. 

CALLING SEQUENCE: CALL RBPRNT 

KEY VAEIIABLES: IN “ Sequence number of "viewer" node 

JN - Sequence number of "viewee" node 

NODE - Array of node numbers 

RBVALI, RBVALS - IR and solar image factors 

REFERENCES : RBOUT 

FILES: NOUT - System output file 
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ROUTINE NAME 


EBPSHD 


DESCRIPTION: This routine determines all possible shadowing surfaces 

between node I and the image of node J in specular sur- 
face MIRROR. 

CALLING SEQUENCE: CALL RBPSHD (RADJ, RADI, POSJ, POSI , NSURF, NSS, JST, 

IN, JN, MIRROR) 


RADJ, RADI - 
^OSJ, POSI - 
NSURP 
NSS 

JST 

IN, JN 

MIRROR 


Radii of sphere enclosing nodes i and J 
Position vectors of enclosing spheres 
Number of shadowing surfaces 
NSURE plus all surface images in specular 
surface MIRROR 

Number of possible shadowing surfaces 
found 

Sequence ninubers of surfaces containing 
nodes I and J 

Sequence number of current specular 
surface 


REFERENCES: RBCAL, RBEXPN, RBCALI,' RBXPNI 


ROUTINE NAME: RBRDIN 

DESCRIPTION: This routine sets up the shadowing data arrays, reads 

form factors from the NFF file written in the FF link 
and writes them on a scratch file for later use, and 
initialises the NFF file for storage of image factors. 


CALLING SEQUENCE; CALL RBRDIN 

KEY VARIABLES ; NSPEC 
NSURF 

Dl>e, TSTR, 
DSTR, IFS, IKS, 
PR, PSH 
ISPEC 

SREFLI, SEEFLS 


REFERENCES ; REMAIN 


” Number of active specular surfaces 
■“ Number of shadowing surfaces 

- Shadowing data arrays 

- Array of specular surface sequence 
numbers 

“ IR and solar specular reflectance 
arrays 


FILES : NFF - Disk file for storage of image factors 

NOUT - System output file 
NRAN - Random access file 
NSCRl - Scratch file 
NRSI, NRSO 
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ROTJTII'JE NAME: RBEDRQ 

DESCRIPTION; This routine initializes the image factor arrays and 
reads restart values if they are available. 

CALLING SEQUENCE: CALL EBRDRQ 

KEY VARIABLES: FFVALI , FFVALS — IR and solar form factors from the FP link. 

Also used for image factor restart values 
RBVALI, EBVALS - IR and solar image factor arrays 

REFERENCES: REMAIN 

FILES: NSCRl - Scratch file i 


ROUTINE NAME: RBROW 


DESCRIPTION; This routine writes image factor data on file NFF by rows. 

CALLING SEQUENCE: CALL RBROW 

KEY VARIABLES: IN - Image factor row number 

NODE - Array of node numbers 
RBVALI, RBVALS - IR and solar image factors 

REFERENCES; REMAIN 

FILES: NFF - Disk file for storing image factors 

KRSO 


ROUTINE NAME: RBRPSN 

DESCRIPTION: This routine calculates the minimum radius of a sphere 

that will encompass a given subnode and determines the 
position vector to the center of the sphere. 

CALLING SEQUENCE; GALL RBRPSN (RADN, POSN, ILK, DATA, BETA, GAMMA, DB, DG) 

RADN - Radius of encompassing sphere 
POSN - Position vector of sphere 
ILK - Surface type of node 

DATA - Surface dimension values 

BETA - Distance from edge of node to center of 
current subnode in the beta direction 
GAMMA - Distance from edge of node to center of 
current subnode in the gamma direction 
DB, DG - Dimensions of subnode in the beta and 
gamma directions 

REFERENCES: RBEXPN, RBXPNI 
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ROUTINE NAME: 


RERSI 


DESCRIPTION: This routine reads rows of image factors from the 

NRSI file and individual image factors from the NRSI 
file. 

CALLING SEQUENCE: CALL RBRSI (I COMP) 

IGOMP “ flag indicating whether or not all 
image factors were found on restart 
tapes. 


REFERENCES ; RBRDIN 

FILES; NRSI, NRTI, NRSO, NFF, NOUT, NRTO 


ROUTINE NAME : RBRSUM 

DESCRIPTION: This routine prints out the time requir: ’ to calculate 

one row of image factors, 

CALLING SEQUENCE; GALL RBRSUM 

REFERENCES: RBROT, RBRSI 

FILES: NOUT - System output file 


REV 1 
III-119 


ROUTINE NAJIE; 


RBSKD 


DESCRIPTION: This routine calculates the elemental shadowing between 

surfaces IN and the image of surface JN in specular surface 
MIRROR. An element is either completely shadowed or not 
shadowed at all. Shadowing is reduced by the transmissivity 
of semitransparent shadowing surfaces. 

CALLING SEQUENCE: CALL RBSHD (ILKI, RX, RY, RZ, POSI, WE, WA, JST, IN, 

JN, I, NSS, RS, MIRROR) 

KEY VARIABLES: ILKI - Surface type of IN 

RX, RY , RZ - Components of vector from element on node I 
to element on node J 

POSI - Array of elemental position vectors 
WE, WA - Elemental shadowing factors for IR and solar image 
factors 

JST “ Number of possible shadowing surfaces 

IN, JN - Surface sequence numbers 

I - Element sequence number 

NSS - Number of passible shadowers (including images) 

RS - Square of the magnitude of the eleraent-to-element 

connecting vector (RX^ + RY^ + RZ^) 

MIRROR Sequence number of current specular surface 

REFERENCES: RBCAL, RBEXPN, RBCALI 


ROUTINE NAME; RBTMCK 

DESCRIPTION: This routine calculates the time remaining in the run and 

compares it with the estimated time for the next calculation. 
If insufficient time remains, the run is aborted. 

CALLING SEQUENCE: CALL RBRMCK (lADDR, MAXLFT) 

lADDR - Not used 

MAXLFT - Estimated time required for next calculation 
KEY VARIABLES: M - Time remaining in run 

REFERENCES; RBCAL, RBEXPN, RBCALI, RBXPNI 
FILES: NOUT - System output file 
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ROUTIHE NAME; 


RBTRS3 



DESGRIPTION: This routine transforms points in an SCS to paints in the ICS, 

BCS, or CCS, 

CALLING SEQUENCE; GALL RBTRS3 (X, Y, Z, A, B, C, EX, RY, RZ, TRAN) 

X, Y, Z - Coordinates of point in the ICS, BCS, or CCS 

A, E, C - Coordinates of point in the SCS 

RX, RY, RZ - Components of the SCS origin position vector 
in the ICS, BCS, or CCS 
TRAN “ Matrix of direction cosines 

REFERENCES: RBELEM, RBEXPN, RBSHD, RBXPNI 


ROUTINE NAME: REVIEW 

DESCRIPTION: This routine calculates the unshadowed form factor between 

node IN and the image of node JN. It also calculates RATI 
and RATJ, which will later be compared with the user-input 
variable, FFRATL, to see if the nodes should be expanded 
Into subnodes for better accuracy. 

GALLING SEQUENCE: CALL REVIEW (NT I, NTJ, POSI, ARAI, NTOTI, NTOTJ, ARI, 

ARJ, NOT, FRACT, RATI, RATJ) 

NTI “ Initial number of elements on node IN 

NTJ - Total initial number of elements on node IN 

and the image of node JN 
POSI - Array of elemental position vectors 

ARAI - Array of elemental area vectors 

NTOTI, NTOTJ - Number of elements required on node IN 
and the image of node JN 

ARI, ARJ - Areas of node IN and the image of node JN 

NOT - Flag to Indicate whether node IN can "see" 

the image of node JN 
= 0 Can see 
= 1 Cannot see 

FRACT “ Unshadowed nod e-to- image form factor 
RATI, RATJ - Ratios of indicated maximum elements re- 
quired to NTOTI and NTOTJ 

REFERENCES: RBCAL, RBCALI 
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ROUTINE NAME: 


RBVNT 


DESCRIPTION: This routine calculates unshadowed form factors between 

subnodes and determines the number of elements required, 
based on a weighted-average criterion. 

CALLING SEQUENCE: CALL RBVHT (NTI, NTJ, POSI, ARAI, NTOTI, NTOTJ, ARI, 

ARJ, NOT, FRACT) 

NTI 

NTJ 


POSI 

ARAI 

NTOTI, NTOTJ 


ARI, ARJ 
NOT 


FRACT 


REFERENCES: RBEXPN , RBXPNI 


ROUTINE NAME : RBXPNI 

DESCRIPTION: This routine is analogous to subroutine RBEXPN. It is 

called when the I - J order of node pairs are reversed 
due to area differences. 

CALLING SEQUENCE: Reference RBEXPN 

REFERENCES : Reference RBEXPN 

FILES: Reference RBEXPN 


- Number of elements on the subnode of 
node IN 

- Total number of elements on the suhnode 
of node IN and the subnode of the 
image of node JN 

** Array of elemental position vectors 

*■ Array of elemental area vectors 

- Number of elements required on the 
subnodes of node IN and the image of 
node JN 

“ Areas of subnodes 

- Flag to indicate whether the subnode 
on node IN can "see" the image of node 
JN 

- 0 Can see 
= 1 Cannot 'lee 

- Unshadowed subnod e-to-subnode form 
factor 
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ROUTINE NAME; RCCMBN 


DESCRIPTION: This routine, given an array defining combinations, 

combines and calls the output routines. 

CALLING SEQUENCE; CALL RCCMBN (ICOMB, ICOMBL, SF, SPACNO, NUNIT) 

ICOMB - Array of combination data 
ICOMBL - Length of ICOMB array 
SF - Temporary array to store script F 
SPACNO - Array to store script F to space 
NUNIT - Unit containing gray^body matrix 

REFERENCES; RCMAIN 
FILES; NSCR3, NSCRl, NSCR2 


ROUTINE NAME; RCDATA 

DESCRIPTION: This routine is user-called in the Operations Data block 

and defines parameters for the RCCAL link. 

CALLING SEQUENCE: CALL RCDATA (NSGBIR, PNCH, FMIN, IRKN, RKSPC, NSPAC, 

SIG, AMPF, TAPE, RADI, NEFI, IPRIM, ISEC) 

KEY VARIABLES: See user's manual. Appendix D 

REFERENCES: Operations data block 
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ROUTINE NAME: RCEND 


DESCRIPTION: This routine can be replaced by the user to intervene 

just prior to the end of the RCCAL link. 

CALLING SEQUENCE; CALL RCEND 

REFERENCES: RCMAIN 

ROUTINE NAME : RCHEAD 

DESCRIPTION: This routine prints the control parameters on the output 

file. 

CALLING SEQUENCE; CALL RCHEAD 
REFERENCES: RCMAIN 

ROUTINE NAME : RCINIT 

DESCRIPTION: Combines nodes according to correspondence data, 

GALLING SEQUENCE: CALL RCINIT (NOCOM) 

NOCOM “ flag to indicate if any combining was 
performed 

REFERENCES : RCMAIN 

FILES : NG3IRR, NOUT 

ROUTINE NAME: RCMAIN 

DESCRIPTION; This is the main driving logic of the RCCAL link and 
directs the main logic flow, 

CALLING SEQUENCE: CALL RCMAIN 

REFERENCES: RCPROG (preprocessor generated) 
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ROUTINE NAME: 


RCPNCH 


DESCRIPTION; 


This routine defines the final output form of the 
RADKs in the form of cards, BCD tape, or binary tape. 


CALLING SEQUENCE: CALL RCPNCH (ICN, NI, NJ, SIG, SFA) 


ICN - Conductor number 
NI - Node i 

NJ - Node j 

SIG “ Stephan-Boltzmann constant 
SFA - Script F area factor 

REFERENCES: RCMA.IN, RCPSS2 

FILES : NOUT 


ROUTINE NAME: RCPRE 

DESCRIPTION: This routine can be replaced by the user to provide 

intervention prior to the RGCAL computation in the 
RCCAL link. 

CALLING SEQUENCE: CALL kCPEE 

REFERENCES : RCMAIN 


ROUTINE NAME : RCPROG 

DESCRIPTION: This routine is generated by the preprocessor. It 

defines all labeled common blocks required by the 
RCCAL link, as well as calls into the main logic. 

CALLING SEQUENCE: CALL RCPROG 

REFERENCES: TRASYS (root segment) 
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ROUTINE NAME ; RCF5S 1 


DESCRIPTION: This routine sorts rows o£ radiation conductors into 

decreasing order and sets up special identification 
flags . 

CALLING SEQUENCE: CALL RCPSSl (NND, SORT, ISORT, IFLG) 

NND - number of values in row radiation 
conductors 

SORT 

ISORT arrays used in sorting process 
IFLG 

REFERENCES : RCMAIN 

FILES: NSCR2, NSCR3, NOUT 


ROUTINE NAME : RCPSS2 

DESCRIPTION: This routine condenses the matrix of radiation conduc- 

tors by saving only significant conductors according 
to user criteria and punches out significant radiation 
conductors. 

CALLING SEQUENCE: CALL RCPSS2 (NND, IGN, ESUM, IFLG) 

NND “ Reference RCPSSl 

ICN - Radiation conductor number 

ESUM - Array of radiation conductor sums 

IFLG - Scratch array 

REFERENCES ; RCMAIN 

FILES: NSCR2, NSCR3, NOUT 
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ROUTINE NAME : RCTRAN 

DESCRIPTION: This routine transposes the gray body factor matrix. 

GALLING SEQUENCE: CALL RCTRAN (A,B,D,E, NB, NN, NR, NRECT, 

ITRC50, MAXBC) 


Scratch arrays 

NBj NR - Size flags 
NN - Number of nodes 
NRECT - Record count flag 
ITRC50 - Trace flag 

MAXBC - Amount of core available for transposing 
operation 

REFERENCES : RCMAIN 

FILES: NSCR2, NGBIR, NOUT, NSCRl 



ROUTINE NAME ; REMAIN 

DESCRIPTION: This routine initializes the variables used in the 

processor via three methods: 

1) Reads in the sequential file; 

2) Sets variables with predefined data; 

3) Reads in the random I/O file, 

CALLING SEQUENCE: CALL RDMAIN 

REFERENCES : RDPROG 

FILES: NSQNTL, NOUT, NRAN, NRSI, NRSO, WRTI 
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ROUTINE NAME: 


EEADHD 


DESCRIPTION: This routine reads the header records on the restart 

tapes. 

CALLING SEQUENCE: CALL READHD (NONIT, LABEL, ICALG) 

NUNIT - Limit to be read (NSRI, NRTI) 

LABEL - Type of data record (FF, DI, RB, 

CM, GB, SF). 

ICALC - Resume calculations flag 
REFERENCES: All restarfcable segments 

FILES : NOUT 


ROUTINE NAME : RETRVE 

DESCRIPTION: This routine retreives node number, area and surface 

property arrays from the NRARR file and places it in 
core. 

CALLING SEQUENCE: CALL RETRVE (IR) 

IR - Flag indicating arrays are combined or 
uncombined 

REFERENCES; Node combining routines 
FILES: NRARR, NOUT 


ROUTINE NAME : RIGPAC 

DESCRIPTION; This routine was written for the Univac computer to 
simulate the GDC random access package. Contains 
entry points NRITHM and READ^S . 

CALLING SEQUENCE: CALL RIOPAC (NUNIT, Fl^fA, NWDS, NR) 

NUNIT “ Unit number 

FWA - Address of first wojrd to read 
NWDS - Number of words to be read 
NR - Record number 


ROUTINE NAME: 


RKDATA 


DESCRIPTION: This routine is user-called in the Operations Data 

block and defines parameters for the RKGAL link. 

CALLING SEQUENCE: CALL RKDATA (NSGBIR, PNCH, FMIN, IRKN, RKSPC, 

NSPAC, SIG 3 AMPF, TAPE) 

KEY VARIABLES: See Users Manualj Appendix D 

f , 

REFERENCES: Operations Data Block 


ROUTINE NAME: RSRCCK 

DESCRIPTION: This routine determines if a parity error was 

encountered or a valid record was just read on a 
restart file. 

CALLING SEQUENCE: CALL RSRCCK (IFILE, lOK) 

IFILE - File name (NRSI, NRTI) 
lOK - Error/ no error flag 

REFERENCES: All restart tape read routines. 

FILES; NRSI, NRTI, NOUT 

ROUTINE NAME ; RSTON 

DESCRIPTION; This entry point in subroutine RESTRT sets up flags 
to resume reading the L:ji 3 ter restart tape, at the 
user's option. 

CALLING SE()UENCE; CALL RSTON (See Users Manual, Appendix D) . 

REFERENCES: ODPROG (Users Operations Data) 
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ROUTINE NAME: RSTOFF 


DESCRIPTION: This entry point into subroutine RESTRT sets up 

flags to interrupt the reading of the master restart 
file at the user's direction. 

CALLING SEQUENCE: CALL RSTOFF (reference Users Manual, Appendix D) . 

REFERENCES: ODPROG (Users Operations Data) 


ROUTINE NAME: SATURD 

DESCRIPTION: This routine sets up the planet parameters concerning 

the planet Saturn. 

CALLING SEQUENCE: CALL SATURD 


KEY VARIABLES: 


PRAD - Planet radius 

SOL - Solar constant 

PALB - Planet albedo factor 

WDS - Planet darks ide emissive power 

WSS " Planet subsolar emissive power 

GRAV - Gravitational constant 


REFERENCES : 


ORBITl, 0RBIT2, DITTP 


ROUTINE NAME: SETBEA 

DESCRIPTION; This routine positions a plot beam at the coordinates 
of a point In the subject mapping. 

CALLING SEQUENCE: CALL SETBEA (X, Y) 

X, Y - Coordinates of point 

REFERENCES; NPAXES, NPFPLT, NPSCAL, OPAXES, OPFPLT, OPSCAL, 

OPVCS , PLDRIV 
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ROUTINE NAME: 


SETFLG 


DESCRIPTION: This subroutine decodes the NMODIR (Z, ID) word to 

obtain FECMB and lAUTOC flags as needed for appro- 
priate node property access. 

CALLING SEQUENCE: CALL SETFLG 

REFERENCES : All segments 


ROUTINE NAME : SFELAV 

DESCRIPTION: This routine computes the position and area vectors 

for elements on a node given the nodal dimensions. 


CALLING SEQUENCE: CALL SFELAV (ILP, DATA, TRAN, RX, RY, RZ, NB, 

NG) 


ILP 

DATA 

TRAN 

RX, RY, RZ 
NB, NG 


“ Surface type 

- Array of surface dimensions 

- Direct cosine matrix for node 

- Position vector for node 

- Number of elements in the fccjo 
directions the node is to be 
divided into 


REFERENCES : SFMAIN 
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ROUTINE NAME: 


SFELEM 


DESCRIPTION: This routine, given the nodal dimensions and the 

total number of elements required, determines the 
number of elements in each direction. 


CALLING SEQUENCE; CALL SEELEM (ILP, DATA, TRAN, RX, RY, R2, 

NTOT) 


ILP 

DATA 

TRAN 

RX, RY, RZ 
NTOT 


- Node type 

- Array of node dimensions 

“ Direction cosine matrix for 
node 

- Position vector for node 

- Total number of elements 
required on node 


REFERENCES: SFMAIN 
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ROUTINE NAJIE: SFELMT 


k 


i 


i 






DESCRIPTION: This routine establishes the minimum number of elements 

on a node to determine if a view is possible. 

CALLING SEQUENCE: CALL SFELMT (NB, NG) 

NB, WG - Minimum number of elements required to provide 
a representative view from the node 


REFERENCES : SFMAIN 


ROUTINE NAME : SFELSL 

DESCRIPTION: This routine, given the total number of elements and the 

nodal dimensions, computes the number of elements in each 
direction to provide near-square elements. 

CALLING SEQUENCE; CALL SFELSL (NB, NG, ILP, DATA, NTOT) 

NB, NG “ Computed number of elements required in various 
directions 

ILP - Surface type 

DATA “ Array of surface dimensions 

NTOT - Total number of elements to be used on the node 

REFERENCES : SFELEM 


ROUTINE NAME; SFGTST 

DESCRIPTION: This integer function routine, given the allowable range on 

gamma concerning shadowing and given a gamma value, returns a 
0 if the gamma value falls within the allowable range or a 1 
if the value is outside the allowable range. 

CALLING SEQUENCE: SFGTST (GN, GX, GT) 

GN, GX - Minimum, maximum range 
GT “ Value to be checked 

REFERENCES: SFSHAD 


i 

■A 
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ROUTINE NAME: 


SFMAIN 


DESCRIPTION: This is the main driving routine in the SFCAL link 

and directs the logic flow necessary to generate a 
shadow tape. 

CALLING SEQUENCE; CALL SFMAIN 
REFERENCES: SFPROG (preprocessor-generated) 


ROUTINE NAME ; SFPAKA 

DESCRIPTION: This routine packs 9 solar waveband shadow-factor 

values into one word and generates 19 words each time 
the routine is called. These 19 words are then written 
to NPLS and NRSO (if requested) once every 10 calls. 
Entry point SFFLSH is called to write the last record 
and complete the writing sequence. 

CALLING SEQUENCE: CALL SFPAKA (TABSHA, ICNT) 

TABSHA - Unpacked shadow factor array 
ICNT - A counter that indicates when shadow 
data for 10 nodes have been packed 


REFERENCES : SFMAIN 


ROUTINE NAME ; SFPAKE 

DESCRIPTION: This routine is analogous to SFPAKA for the infrared 

waveband. 

CALLING SEQUENCE; CALL SFPAKE (TABSHE, ICNT) 

REFERENCES : SFMAIN 
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ROUTIM NAME: SFPL ' [ 

DESCRIPTION: This routine computes ‘ thh shadow factor to a planet 

element when Interpolation of the shadow tables is 
doubtful. 

CALLING SEQUENCE; GALL SPPL CASEE,ASFA, lEP, RES) 

AS EE- - Infrared waveband shadow factor 

ASFA “ Solar waveband shadow factor 

lEP - Planet element pointer | 

KBS ~ Node can be shaded flag 8 i 

REFERENCES : DICALP 

ROUTINE NAME: SFPROG 

DESCRIPTION; This routine is generated by the preprocessor into 
the SFCAL link, and provides the necessary labeled 
commons, 

CALLING SEQUENCE: CALL SFPROG 

REFERENCES; TRASYS (root segment) 

ROUTINE NAME; SFPRTR 

DESCRIPTION; This routine px-ints out shadoxj factor tables that 
are read in from the master restart tape. 

CALLING SEQUENCE: CALL SFPRTR 

REFERENCES: SFMAIN 

FILES; NSCRl, NOUT 

ROUTINE NAME : SEPSHS 

DESCRIPTION: This routine determines which surfaces could 

possibly shadow between node i and the sun. 
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CALLING SEQUENCE: CAIL SFPSHS (RADS, POS, iUNP, NSDRF, NSHAD, IN) 

RADS - Radius of sphere encompassing node i 
POS - Position vector of encompassing sphere 
SUNP ~ Position vector for sun 
NSURF - Number of shadowing surfaces 
NSHAD - Number of possible shadowers 
IN “ Sequence number of node i 

KEY VARIABLES: ISHAD - Array of possible shadowers 

REFERENCES: SFMAIN 


ROUTINE NAME: SFRDIN 

DESCRIPTION: This routine checks for restart information and sets the 

proper flags for the program. The shadowing surfaces are 
also set up in terms of the labeled commons. 

CALLING SEQUENCE: CALL SFRDIN 

REFERENCES: SFMAIN 


ROUTINE NAME: SFRDRQ 

DESCRIPTION: This routine reads in the restart data for one node at a 

time. 

CALLING SEQUENCE; CALL SFRDRQ (IFLG, TABSHA, IN) 

IFLG Data check flag 

5= 2HN0 No data found 
= 3HYES Data found 

TABSHA - Array of restart data for node i 
IN - Sequence number for node ± 


REFERENCES; SFMAIN 
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ROUTINE NAME: SFSHAD 

DESCRIPTION: This routine computes element-to-element shadowing between 

node i and the sun. 

CALLING SEQUENCE: CALL SFSHAD (RX, RY, RZ, ME, WA, IN, I, NSS, RS) 

RX, RY, RZ - Components of vector from element on 
node i to sun 

WE, WA “ Elemental shadowing factors for IR and 
solar 

IN - Surface sequence number 

I “ Element sequence number 

NSS “ Number of possible shadowers 

RS - Square of the magnitude of the element- 

to-sun connecting vector 

REFERENCES: SFMAIN 


ROUTINE NAME: SFTRS3 

DESCRIPTION: This routine transforms points in an SCS to points in the 

ICS, BCS, or CCS. 

CALLEJG SEQUENCE: CALL SFTRS3 (X, Y, 1, A, B, C, RX, RY, RZ, TRAN) 

X, Y, Z - Coordinates of a point in the ICS, BCS, 
or CCS 

A, B, C - Coordinates of the point in the SCS 
RX, RY, RZ - Components of the SCS origin position 
vector in the ICS, BCS, or CCS 
TRAN - Matrix of direction cosines 

REFERENCES; SFELAV, SFELEM 


ROUTINE NAME; SFUNCT 

DESCRIPTION: This function, given a clock angle and a cone angle, inter- 

polates the packed clock angle-cone angle array located in 
labeled common, computes the function value, A and stores the 
value in SFT. 

CALLING SEQUENCE: A = SFUNCT (CL, CO, SFT) 

CL - Clock angle 

CO - Cone angle 

SFT, A - Interpolated value 


REFERENCES: DICALS, DICALP 
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ROUTINE NAME: SHADE'! 


DESCRIPTION: This routine computes the shadow points, given the orbit 

definition. 

CALLING SEQUENCE: CALL SHADPT 

REFERENCES : OEBITl , 0RBIT2 

ROUTINE NAME: SINX 

DESCRIPTION: This routine computes the sine of an argument, in degrees. 

galling SEQUENCE: CALL SINX (A) 

A - Argument, in degrees 

ROUTINE NAME: SKFILE 

DESCRIPTION: This routine advances a unit N number of files. 

calling SEQUENCE: CALL SKFILE (NUNIT, N) 

NUNIT - Unit number 
N - Number of files 

REFERENCES: FINDST, STORE, QOMAIN , DITTP, PLLOAD 

ROUTINE NAME: SORTDL 

DESCRIPTION: This routine numerically sorts a doublet array within 

itself. 

calling SEQUENCE; CALL SORTDL (lA, NA) 

lA “ Doublet array 

NA - Total number of elements 

REFERENCES: QOMAIN, PLLOAD, RCPSSl 
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ROUTINE NA:1E: SORTS 


DESCRIPTION: This routine numerically sorts a singlet array. 

CALLING SEQUENCE: CALL SORTS (A, JJ) 

A - Singlet input array 

JJ - Number of elements in the array 

REFERENCE : NDUPCK 


ROUTINE NAME: SPIN 

DESCRIPTION: This routine is user-callable and allows the user to define 

the spin axis, spin rate, and start time. 

CALLING SEQUENCE: CALL SPIN (CLOG, CON, RAT, ANGLE, TIMS) 

CLOC - Clock angle locating spin axes 
CON - Cone angle locating spin axes 

RAT “ Spin rate 

ANGLE - True anomaly angle where spin begins 
TIMS - Start time of spin 

REFERENCES: User's Operations Data block 


ROUTINE NAME : STFAQ 

DESCRIPTION: This routine is user-called in the Operations Data block 

to generate a duplicate orbit point, given a true anomaly/ 
time and a step number to retrieve the data from. 

CALLING SEQUENCE: CALL STFAQ (ANGLE, TIM, NST) 

ANGLE - True anomaly the data are to be stored under 
TIM - Current orbital time data are to be stored under 
NST - Step number to retrieve data from 

REFERENCES: Operations Data block (user-called) 

FILES: NDI, NOUT, NSCRl, NTQ 
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ROUTINE NAME : STFFLG 

DESCRIPTION: This routine stuffs bits into NMODIR (Z, ID) for 

the FFCMB and lAUTOC flags. 

CALLING SEQUENCE: CALL STFFLG (NBITS) 

NBITS “ Bit location of FFCMB and lAUTOC 
flags in NMODIR (Z, ID) word. 

REFERENCES: All segments. 


ROUTINE NAME : STORE 

DESCRIPTION: This routine generates the header record on the 

required file. If the file has been rewound or read, 

STORE repositions the file to the proper point. 

CALLING SEQUENCE: CALL STORE (NUNIT, ISTEP, LABELl, LABEL2) 

NUNIT - Unit to write header on 
ISTEP - Step number 
LABELl - Identifier 
LABEL2 “ Identifier 

REFERENCES: STFAQ, FFRDIN, SFRDIN, DIEND, DIRDRQ, GBSCFA, AQMAIN, ADSURF, 

CMCMBN, CMRSI, GBAPRX, GBRSI, RBRDIN, • RBRSI 

ROUTINE NAME: SUND 

DESCRIPTION: This routine defines the parameters concerning the 

sun. 

CALLING SEQUENCE: CALL SUND 

KEY VARIABLES: PRAD - Planet radius 

RSUN “ Sun radius 
WSUN - Sun enmiissive power 
PALB - Sun albedo factor (0.0) 

REFERENCES: ORBITl, ORBIT2, DITTP 
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ROUTINE NAME : SURER 


DESCRIPTION: This routine sets up the spin rate, spin axis, param- 

eters and computes dawn and dusk times for the planet 
surface option, 

CALLING SEQUENCE: CALL SURER (PNAM, AL, SUI^, AT) 

PNAM - Planet name 

AL - Latitude on surface 

SUNL - Sun latitude 

AT - Atmospheric extinction coefficient 

REFERENCES: User-called 
FILES : NOUT 


ROUTINE NAME: SYMBOL 

DESCRIPTION: This routine generates Hollerith data on the plot frame. 

The sequence of characters is terminated by a ($.). 

CALLING SEQUENCE: CALL SYMBOL (A) 

A - Array of Hollerith data 

REFERENCES; NPMAIN, NPAXES, NPINFO, OPMAIN, OPAXES, OPINFO, OPVCS, 
PLGRID, PLSYMB 


:> 
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ROUTINE NAME: 


TPERR 


DESCRIPTION: This routine writes out a message when an error is 

encountered reading a data file. 

CALLING SEQUENCE: CALL TPERR (NUNIT, ISUB) 

NUNIT - Data file where error was encountered 
ISUB - Name of subroutine attempting to 
read file 


REFERENCES: All segments. 

FILES : NOUT 


{ 


ROUTINE NAME : TRANS F 

DESCRIPTION: This routine comhines the direction cosine matrices 

to generate a final C-S and S-C matrix. 


CALLING SEQUENCE: CALL TRANS F 

TRANBS) 

RMASS 

TRANCB 

BX, BY, BZ 
TRANCS 

TRANBS 


REFERENCES : BUILDG 


(RMASS, TRANCB, BX, BY, BZ, TRANCS, 

“ Array of surface data 

- Direction cosine matrix central 
to block system 

- Block coordinate component vector 

- Direction cosine matrix central 
to surface 

“ Direction cosine matrix block 
to surface 


ROUTINE NAME : TRASYS 

DESCRIPTION: This is the root segment defined by the preprocessor 

enabling the operations data to direct all logic flow. 


CALLING SEQUENCE: CALL TRASYS 
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ROUTINE NAME: 


TRNSP 


DESCRIPTION: This routine transposes matrix A and stores the result 

in matrix B . 

CALLING SEQUENCE; CALL TRANS P (A, B) 

A - Input matrix 
B - A^ 


REFERENCES : DILOC 

ROUTINE NAME : URANUD 

DESGRIP^iON: This routine defines the parameters necessary to 

orbit Uranus. 


CALLING SEQUENCE: GALL URANUD 


KEY VARIABLES : PRAD 

“ Planet radius 

SOL 

“ Solar constant 

PALB 

- Planet albedo factor 

WDS 

- Darkside emissive power 

WSS 

“ Subsolar emissive power 

GRAV 

- Gravitational constant 

REFERENCES: ORBITl, 

ORB IT 2, DITTP 


ROUTINE NAME : VENUS D 

DESCRIPTION; This routine defines the parameters necessary to 
orbit Venus. 


CALLING SEQUENCE: CALL VENUS D 


KEY VARIABLES: PRAD 

- Planet radius 


SOL 

- Solar constant 


PALB 

- Albedo factor 


WDS 

“ Darkside emissive 

Power 

WSS 

" Subsolar emissive 

power 

GRAV 

- Gravitational constant 

REFERENCES: ORBITl, 

0RBIT2, DITTP 



REV 1 
III-143 


ROUTINE NAME ; WRAPUP 


DESCRIPTION: This routine sets up the normal exit from the TRASYS 

processor and prints some information messages. 

CALLING SEQUENCE: CALL WRAPUP 

REFERENCES: TRASYS (root segment) 

FILES : NOUT 


ROUTINE NAME : WRITHD 

DESCRIPTION; This routine writes header records to the restart 
tapes . 

CALLING SEQUENCE: CALL WRITHD (NUNIT, LABEL) 

NUNIT - Name of file to write header record on 
LABEL - Type of information to follow header 

REFERENCES: All restart segments. 


ROUTINE NAME : WRITHT 

DESCRIPTION: This routine writes trailer records after informa- 

tion on the restart files. 

CALLING SEQUENCE: CALL WRITHT (NUNIT, LABEL) 

I REFERENCES: All restart segments. 


ROUTINE NAME ; ZNPMAX 

DESCRIPTION; This function routine is called from the node plotter 
(NPLOT) and scales the nodes to be plotted. 

CALLING SEQUENCE: ZNPMAX (NNP) 

NNP “ Number of surface/nodes to be selectively 
plotted 

REFERENCES: NPMAIN, NPCONV 
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FILE DEFINITIONS - PROCESSOR LIBRARY 


FILE NAME : BCDOL 

PROGRAM VARIABLE NAME: NBCDOU 

UNIT REFERENCE (UNIVAC/JSC) : 29 

PURPOSE: This file is the output tape for BCD data in the thermal 

analyzer input format. 

SEGMENT REFERENCES; QOCAL (WRITE) 

RKCAL (WRITE) 

RCCAL (WRITE) 


FILE NAME: DI 

PROGRAM VARIABLE NAME; NDI 

UNIT REFERENCE (UNIVAC/JSC): 10 

PURPOSE: This file is for storage of solar, planetary, and albedo 

direct Irradiation data, 

SEGMENT REFERENCES: DICAL (WRITE) 

DRCAL (WRITE) 

AQCAL (READ) 

PLOT (READ) 


FILE NAME; ^ 

PROGRAM VARIABLE NAME; NDIR 
UIvflT REFERENCE (UNIVAC/JSC); 22 

PURPOSE: This file is for storage of direct irradiation data input 

through the flux data block. 

SEGMENT REFERENCES: DICAL (READ) 



FILE NAME : FF 


PROGRAM VARIABLE NAME: NFF 

UNIT REFERENCE (UNIVAC/JSC) : 9 

PURPOSE: This file is used to store form factor data. 

SEGMENT REFERENCES: FFCAL (WRITE) 

GBCAL (READ) 

RBCAL (WRITE) 


FILE NAME : FFR 

PROGRAM VARIABLE NAME: NFFR 

Limit reference (UNIVAC/JSC): 21 

PURPOSE: This file is used to store form factor data input through 

the form factor data blocks. 


SEGMENT REFERENCES: FFCAL (READ) 


FILE NAME : GBIR 

PROGRAM VARIABLE NAME: NGBXR 

UNIT REFERENCE (UNIVAC/JSC) : 11 

PURPOSE; This file is for storage of infrared waveband gray-body 
factor data. 

SEGMENT REFERENCES: AQCAL (READ) 

GBCAL (WRITE) 

ODPROG - Subroutine GBAPRX (WRITE) 

RCCAL (READ) 

RKCAL (READ) 
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FILE NAME: GBIRR 


PROGRAM VARIABLE NAME: NGBIRR 

UNIT REFERENCE (UNI VAC/ JSC) : 23 

PURPOSE: This is the correspondence data storage file. 

SEGMENT REFERENCES: RCCAL (READ) 

RICCAL (READ) ■ 

QOCAL (READ) 


FILE NAME : GBSQ 

PROGRAM VARIABLE NAtlE: NGBSO 

UNIT REFERENCE (UNIVAC/JSC) : 12 

PURPOSE: This file is for storage of solar-waveband gray-body factor 

data. 

SEGMENT REFERENCES; AQCAL (READ) 

GBCAL (WRITE) 

ODPROG - Subroutine GBAPRX (^^^RITE) 

FILE NAME : OUTPUT 

PROGRAM VARIABLE NAME: NOUT 

UNIT REFERENCE (UNIVAC/JSC): 6 

PURPOSE; This is the print output file, 

SEGMENT REFERENCES; All 
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FILE NAFIE: PLS 


PROGEAM VARIABLE NAME: NPLS 

UNIT REFERENCE (UNIVAC/JSC) : 13 

PURPOSE; This is used to store the spacecraft/planet form factor 
matrix and shadow factor data. 

SEGMENT REFERENCES: DICAL (READ/WRITE) 

SFCAL (READ /WRITE) 


FILE NAME ; PLSR 

PROGRAM VARIABLE NAME: NPLSR 

UNIT REFERENCE (UNIVAC/JSC) ; 25 

PURPOSE: This file is used to store data from the shadow factor data 

hlock, 

SEGMENT REFERENCES; SFCAL (READ) 


FILE NAME : PUNCH 

PROGRAM VARIABLE NAME; NPUN 

UNIT REFERENCE (UNIVAC/JSC); 7 

PURPOSE: This file is the punch output file. 

SEGMENT REFERENCES: DICAL (WRITE) 

FFCAL (WRITE) 

RCCAL (WRITE) 

RKCAL (WRITE) 

QOCAL (WRITE) 
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"LE NAME: RIO 


i^IOGRAM VARIABLE NAME: NRAN 

UNIT REFERENCE (UNIVAC/JSC) : 8 

PURPOSE: This file is the primary random access file 3 which is 

used to store all node and surface description data. 

SEGMENT REFERENCES: DICAL (READ) NPLOT (READ) 

DRCAL (READ) OPLOT (READ) 

FFCAL (READ) RBCAL (READ) 

FILE NAME : RSI 

PROGRAM VARIABLE NAME: K7 31 

UNIT REFERENCE (UNIVAG/JSC) : 14 

PURPOSE; This file is the master restart input tape. 

SEGMENT REFERENCES: DICAL (.READ), FFCAL (READ), GBCAL (READ), 

SFGAL (READ), CMCAL (READ), RBCAL (READ), 
DRCAL (READ) 

FILE NAME: RSO 

PROGRAM VARIABLE NAME: NRSO 

UNIT REFERENCE (UNIVAC/JSC) : 15 

PITRPOSE: This file serves as the restart output tape. 

SEGMENT REFERENCES: DICAL (WRITE), FFCAL (WRITE), GBCAL (WRITE), 

SFCAL (WRITE), RBCAL (WRITE), DRCAL (WRITE) 
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FILE NAME : RTI 

PROGRAM VARIABLE NAME: NRTI 

UNIT REFERENCE (UNIVAC/JSG) : 17 

PURPOSE: This file is the temporary restart input tape. 

SEGMENT REFERENCES: DICAL (READ), FFGAL (READ), GBGAL (READ) 

RBCAL (READ), DRCAL (READ 


FILE NAME : RTO 

PROGRAM VARIABLE NAME: NRTO 

UNIT REFERENCE (UNIVAC/JSC) : 18 

PURPOSE: This file is the temporary restart output tape, 

SEGMENT REFERENCES: DICAL (WRITE), FFCAL (WRITE), GBGAL (WRITE) 
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FILE NAME: SQNIL 


PROGRAM VARIABLE NAME; NSQNTL 
UNIT REFERENCE (UKIVAC/JSC) ; 16 

PURPOSE: This file contains painters for the random access file, plus 

miscellaneous flags and quantities generated in the preprocessor 
for use by processor segments. 

SEGMENT REFERENCES; RDPROG (READ) 


FILE NAME: TAPEl 

PROGRAM VARIABLE NAME: NSCRl 

UNIT REFERENCE (UNIVAC/JSC) : 1 

PURPOSE: This file is scratch file 1. Scratch files are never used 

to pass information between segments. 

SEGMENT REFERENCES: GBCAL (READ/WRITE) 

RCCAL (READ/WRITE) 

RKCAL (READ/WRITE) 


FILE NAME: TAPE 2 

PROGRAM VARIABLE NAME: NSCR2 

LTIIT REFERENCE (UNIVAC/JSC) : 2 

PURPOSE: This file is scratch file 2 . 

SEGMENT REFERi'NCES; GBCAL (READ/WRITE) 

PLOT (READ/WRITE) 
QOCAL (READ/WRITE) 
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FILE NAME: TAPES 


PROGRAM VARIABLE NAME; NSCR3 

UNIT REFERENCE (UNIVAC/JSC) ; 3 

PURPOSE: This file is scratch file 3. 

SEGMENT REFERENCES; FFCAL (READ/WRITE) 

GBCAL (READ/WRITE) 
RCCAL (READ /WRITE) 


FILE NAME: ^ 

PROGRAM VARIABLE NAME: NTQ 

UNIT REFERENCE (UNIVAC/JSC) : 14 

PURPOSE; This file is used for storage of absorbed heat data, 

SEGMENT REFEP,ENCES: AQCAL (WRITE) 

QOCAL (READ) 

PLOT (READ) 


FILE NAME; TQR 

PROGRAM VARIABLE NAME: NTQR 

UNIT REFERENCE (UNIVAC/JSC) : 26 

PURPOSE; This file is a restart file for absorbed heat data. 
SEGMENT REFERENCES: Not currently used. Name and unit reserve 

FILE NAME: TRAJ 

PROGRAM VARIABLE NAME: NTRAJ 

UNIT r.,E^”ENCE (UNIVAC/JSC): 4 

PURPOSE: This file is used to input trajectory tape data. 

SEGMENT REFERENCES: ODPROG - Subroutine DITTP (READ) 



FILE NAME: USERl 


PROGRAM VARIABLE NAFIE: NUSERl 

UNIT REFERENCE (UNIVAC/JSC) : 19 

PURPOSE: This is a scratch file reserved for the user 

SEGMENT REFERENCES: User option 

FILE NAME: USER2 

PROGRAM VARIABLE NAME: NUSER2 

UNIT REFERENCE (UNIVAC/JSC) r 20 

PURPOSE: This is a scratch file reserved for the user 

SEGMENT REFERENCES: User option 
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D. VARIABLE DEFINITIONS - PROCESSOR LIBRARY 

LABELED COMMON /ALPK/ 

This common block contains an array of nodal absorptivities in the 
solar waveband . 

ALPH “ An array of solar absorptivities for active nodes 


LABELED COMMON /AQQDP/ 

This common block provides a storage area in the AQPROG segment 
for incident planetary fluxes that are read in from the NDI file. 

QDP - An array of incident planetary fluxes 
LABELED COMMON /AQQDR/ 

This common block provides a storage area in the AQPROG segment 
for incident albedo fluxes that are read in from the NDI file. 

QDR - An array of incident albedo fluxes 
LABELED COMMON /AQQDS/ 

This common block provides a storage area in the AQPROG segment 
for incident solar fluxes that are read in from the NDI file. 

QDS - An array of incident solar fluxes 


LABELED COMMON /AQTEMP/ 

This common block provides temporary storage in the AQPROG segment 
for the node array, as read from the NDI file, for verifying 
c '’vrespondence data with the active model node array, NODE. 

ITEMP “ Temporary array of node numbers 
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LABELED COMMON /AREA/ 

This common block contains an array of the active model nodal areas. 
AREA ~ An array of nodal areas 


LABELED CO>MON /ARE AT/ 

This common block provides a temporary working and storage area 
for use in the QOPROG segment when combining nodal areas. 

AREAT " An array of combined nodal areas 


LABELED COMMON /ARRAYS/ 


This common block is set up by the preprocessor and contains 
user-input arrays from the Array Data block. 


ADUMMY 
NAMEl 
NAMEl (N) 
NAME2 


- Integer count of the first array 

- First user-input array 

- Integer count of the second array 

- Second user-input array 


LABELED COMMON /BCSN/ 

This common block contains all the user-input block coordinate 
system (bCS) names and/or the default BCS name, ALLBLK. 
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LABELED COMMON /BLOCK/ 

This common block contains blocking information used by blank 
common to invert a matrix Xijhen calculating gray-body factors. 

NBLCK - Maximum number of blocks into x^hich a matrix may be 
divided (set by a data statement in GBMAIN) 

NBUP - Minimum number of blocks in an upper triangular matrix, 
as determined by full utilization of blank common 
NBLO “ Minimum number of blocks in a lower triangular matrix, as 
determined by full utilisation of blank common 


LABELED COMMON /CAL280/ 

This common block contains plotting information. 

XC - The X-coordinate of the point to be plotted in the plotter 
coordinate system 

YC - The Y-coordinate of the point to be plotted in the plotter 
coordinate system 

IP “ The index on a computed GO TO statement used in NPFPLT and 
OPFPLT 


LABELED COMMON /BFA/ 

This common block contains form factor blockage (shadow) factors 
in the solar waveband. 


LABELED COMMON /BFE/ 

This common block contains form factor blockage (shadow) factors 
in the infrared waveband. 


LABELED COMMON /BLKDIR/ 


This common block contains the block corrdinate system directory. 


NBLKDR 

<1. 

I) 

NBLKDR 

(2, 

I) 

NBLKDR 

(3, 

I) 

NBLKDR 

(4, 

I) 


- Block coordinate system name 

- Length of the block data X'/ritten on the random 
access record 

~ Random access record number 

- Flag to indicate whether or not the BCS trans- 
formation data have been applied 
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LABELED CO^iMON / C CONST/ 


This conimon block contaias a collection of constants and variables 
used throughout the processor. 


AQPRNT 

ARAD 

DELCT 

DIACC 

DIACGS 

DIPNOSH 

DIPNCH 

DLTLNE 

DTR 

ELPBEA 

FFACC 

FFACCS 

FFCMB 

FFDISF 

FFMIN 

FFNAC 


FFNOSH 

FFPNCH 

FFPRKT 

FFRATL 


FFZERO 

FOG 

GAUSS 

GBWBND 

lAI 


lALBFL 

lAQGBI 


- Print/no print flag for abs.:rbed flux data printout 

- Not used 

“ Not used 

- Element selection accuracy factor for node-to-planet 
form factors 

- Element selection accuracy factor for direct flux 
shadowing calculations 

“ Shadow/no shadow flag for direct flux calculations 

- Direct flux punch flag 

- Not used 

- Conversion factor for degrees to radians 

- Not used 

- Element selection accuracy factor for node-to-node form 
factor calculations 

“ Element selection accuracy factor for form-factor 
shadowing calculations 

- Form factor combined/not combined flag 

- Flag to indicate whether FF, DI or SF segments have 
been executed 

“ Minimum form factor value to be saved. All form factors 
smaller than FFMIN are set to zero 

- Flag to allow the input of form factors from an RSI 
file with a node array that does not match the node 
array that resides in core 

- Shadow/no shadow flag for form-factor calculations 

- Form-factor punch flag 

- Form-factor print flag 

- Maximum allowable ratio of the maximum number of elements 
indicated on a node pair divided by the number of elements 
indicated by arithmetic averaging. If this ratio exceeds 
FFRATL, the two nodes are temporarily subdivided 

- The initialization value for the form factor request 
matrices (defaults to -1.0) 

- Not used 

- Not used 

Waveband definition for gray-body calculations 

- Step number indicating where the IR gray-body factor 
matrix is to be obtained for flux calculations using the 
ORBGEN option 

- Albedo flux compute/stuff flag 

- Step number from which IR gray-body factors are to be 
obtained from absorbed-Q calculations 

- Step number from ti?hich solar gray-body factors are 
obtained for absorbed-Q calculations 
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lAQGBS 


to be 



lAQSDA 

lAQSD? 

lAQSDS 

IAS 


lAUTOC 

ICtfflL — rrF— &embj DJ -M Ions dfifi ned in correspondence data 

lEQFF - Flag pointing to equivalent form factors 

IFFSHO “ Flag indicating whether or not to compute form factors 

to shadower-only nodes. Also used to specify whether or 
not shadower-only nodes are to be plotted. 

IGBSFF Step number from which form factors are to be obtained 
for use in gray-body calculations 
ILLUMN - Not used 

IMESS - Starting address of the secondary MESS node array 

INCORE - Flag indicating whether combined or uncombined node, 

area, and properties arrays are in core 
INTMF - Not used 

lOVL “ Index on a computed GO TO statement set in ODPROG and 
used in TRASYS, This determines which segment is to 
be called 

IPRDMP - Flag used to specify the waiting of node, area, and 
properties arrays on the RSO file following calls to 
BUILDC, ADD, or MODXX routines 

IQOARY - Array of step numbers where absorbed-Q data are stored 
IQOCOR - Step number from which correspondence data are to be 
obtained 

IQOTAF - Initial array number for the output Qs (= IQOTME +1) 

IQOTME - Output time array number 

IPIiAFL - Planetary flux compute/stuff flag 

IRKCN “ Initial radiation conductor number 

IRKNGB “ Step number from which gray-body factors are to be 

obtained for use in radiation conductor calculations 
IRKNSP - Space node number 
IRSI - RSI file designator 

IRTI - RTI file designator 

ISFAC - Flag to specify whether or not direct flux shadow factors 
are to be written on the RSO file 

ISKIP - Flag used to skip data on the RSO file until the proper 
header file is found 
ISOLFL “ Solar flux compute/stuff flag 

ISPND - Starting address of the primary MESS node array 
ISTRT - Not used 

ITEIALL - Not used 


- Step number from which direct albedo fluxes are to be 
obtained for absorbed-Q calculations 

- Step number from which direct planetary fluxes are to 
be obtained for absorbed-Q calculations 

“ Step number from which direct solar fluxes are to be 
obtained for absorbed-Q calculations 

- Step number indicating where the solar gray-body factor 
matrix is to be obtained for flux calculations using the 
ORBGEN option 

- Flag indicating whether to include autocombining (polygon 
combining) in form factor combining operations 
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ITRCAO 

ITRCBQ^-■ 

ITRCCO 

ITRCDO 

ITRCIO 

ITRC20 

ITRC30 

ITRC40 

ITRC50 

ITRC60 

ITRG70 
ITRC80 
ITRC90 
MAXBC 
Iffi ' 
MFCO 
MFLUX 
MG 

HITS IN 

MNND 

MRSP 

MSRF 

NBGDSK 

NBLKLN 

NCONST 

NELN 

NERN 

NFFTYP 

NFIGCO 

NFIGFF 

NFIGGB 

NIBBLE 

NMESS 

NMIR 

NMODLS 

NN 

NNOD 

NNODC 

NNODU 

NRMOD 

NRSP 

NS 

NSFO 
NS FT 


- Trace flag for the AWPROG segment 
-=- Trace flag for the QOPROG segment 

- Trace flag for the RBPROG segment (not used) 

- Trace flag for the PLPROG segment 

- Trace flag for routines FINEST and STORE 
“ Trace flag for BUILDC 

- Trace flag for the FFPROG segment (not used) 

- Trace flag for the SFPROG segment (not used) 

- Trace flag for the NPPROG and RCPROG segments 

“ Trace flag used to print the original radiation conductors 

in the RCPROG segment 
“ Trace flag for the DIPROG setraent 

- Trace flag for the gray-body calculations 

- Trace flag for the RKPROG segment 

- Length of blank common 



- Not used 

- Not used 

- Not used 

- Flag to specify MITAS or SINDA output formats 

- Total number of nodes defined in the surface data block 

- Not used 

- Not used 

- Not used 

- Number of block coordinate systems in the model 

- Not used 

- Not used 

•' Effective radiation node number 

- Type flag for form factors (FF, RB, CM) 

- Configuration name for correspondence data access 

- Configuration name for form factor access 

- Configuration name for gray-body factor access 

- Not used 

- Number of MESS node pairs 

- Not used 

- Number of configurations created by BTJILDC/ADD sequences 
(maximum of 20 allowed) 

- Total number of nodes defined by BUILDC/ADD 

- Number of active nodes in the model 

- Number of nodes after combining 

~ Number of nodes before combining 

“ Flag defining the modification level of an RSO/RSI 
restart file 

- Not used 

“ Total number of surfaces defined by BUILDC/ADD 

- Configuration name to be used to identify the file when 
writing the shadow-factor output tape (SHADO) 

- Configuration name to Identify the desired file on the 
shadow-factor input tape (SHADI) 


NSPND - Number of MESS node pairs plus the number of special nodes 

NSSTEP - Sequence number of the current substep 

NSTEP - Sequence number o£ the current step 

NSTPDI - Step number reference for DI data in DRCAL | 

NSTPL - Step number indicating where planetary fluxes are to be 

obtained for a planet-oriented case using the ORBGEN option 
NSTSOL - Step number indicating where solar fluxes are to be 

o obtained for a sun-oriented case using the ORBGEN option 
NSURP “ Number of active shadowi'ng surfaces in the model 
PI - The constant 

QOAMPF - Area multiplying factor for the output Qs 
OOFMPF - Energy multiplying factor for the output Qs 
QOPNCH - Flag to punch the output Qs 
QORMPF - Not used 

QOTAPE - Flag to write the output Qs on a BCD tape 
QOTMPF - Time multiplying factor 

QOTYPE - Flag to outpi t the Q tables, the orbital average Q, or both 
RALB ” Multiplying factor for the absorbed albedo heat 

RFRAC - Significant radiation fraction used in the RCPROG segment 

RKAMPF - Area multiplying factor used in the RKPROG and RCPROG 
segments 

RKMIN - Minimum value of / that xifill result in a valid radia- 
tion conductor 

RKPNCH - Radiation conductor punch flag 

RKSP - Flag for calculating radiation conductors to space 

RKTAPE - Flag to write radiation conductors to the BCD tape 

RFIAN - Multiplying factor for the absorbed planetary heat 

RSOLAR - Multiplying factor for the absorbed solar heat 

RTD " Conversion factor for radians to degrees 

RTOL - Tolerance flag for form factor combining (see Users Manual, 
Appendix D, subroutine RCDATA) 

SAGS “ Not used 

SFPRNT - Print/no print flag for shadow factors on restart tape 
SIGMA - S tefan-Boltzman constant 

STRAGK - Not used 

THGHT - Not used 

TPIAM - Not used 

TRUANF - True anomaly of the final point in the orbit as defined 

using the ORBGEN option 

TRUANI - True anomaly of the first point in the orbit as defined 
using the ORBGEN option 
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LABELED COMMON /DIMS/ 

This common block contains an array of position vectors for 
shadowing surfaces. 

DIMS - Position vector locating the origin of the shadox^ing-sur- 
face coordinate system in the central coordinate system 


LABELED COMMON /DIRCT/ 

This common block contains the translation and rotation informa- 
tion relating a block coordinate system (BCS) to the central 
coordinate system (CCS) , as well as a directory of random-access 
record numbers for the surfaces in the BCS. 


DIRCT 

(1) “ 

X-component of 

the 

translation 

vector 

DIRCT 

(2) - 

Y-component of 

the 

translation 

vector 

DIRCT 

(3) - 

Z-component of 

the 

translation 

vector 

DIRCT 

(4) - 

Rotation about 

the 

CCS X-axis 


DIRCT 

(5) - 

Rotation about 

the 

CCS Y-axis 


DIRCT 

(6) - 

Rotation about 

the 

CCS Z-axis 


DIRCT 

(7) 

1 




DIRCT 

(8) - 

) Order of rotation 

about the CCS X-, Y 

DIRCT 

(9) _ 

1 respectively 




DIRCT 

(10) ■ 






Random-access record numbers 


LABELED COMMON /DISURl/ 


This common block is used to store surface data read in from the 
random-access file in the DIPROG segment 


COMMl 

CSHDI 

DATAI 

DTE I 

DUMI (16) 
ILKI 
INOD 
KSI 

NRMSSl (1) 

OLDAI 

POSNI 


- Five 6-character words describing the surface 

- Can-shade flag 

- Array of the five surface description parameters, 

ALPHA, BMIN, BMAX, GMIN, and GMAX 

- Date of the run 

- Dummy array 

- Identifier for type of surface 

- Can-be-shaded flag 

- Sequence number of surface 

- Surface ID number 

- Surface area 

- Position vector locating the center of the encompassing 
sphere in the CCS 
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PROP I 


RADI 

RXI 

RYI 

RZI 

SPRII 

SPRSI 

IRANI 


- Array of surface optical properties: solar absorptivity 

(ALPH) , Infrared emisslvity (EMISS), infrared trans- 
missivity (TIR) , and solar transmissivity (TSO) 

- Radius of the preshadowing sphere encompassing the 
surface 

- I Components of the position vector locating the origin of 

- ) the SCS in the CCS 

- Specular reflectance of the surface in the IR waveband 

- Specular reflectance of the surface in the solar waveband 

- Transformation Matrix relating the SCS to the CCS 


LABELED COMMON /DRQDP/ 

This common block provides a storage area in the DRCAL segment for 
incident planetary fluxes that are read in from the NDI file. 

QDP - An array of incident planetary fluxes 


LABELED COMMON /DRQDR/ 

This common block provides a storage area in the DRPROG segment 
for incident albedo fluxes that are read in from the NDI file. 

QDR - An array of incident albedo fluxes 
LABELED COM^^ON /DRQDS/ 

This common block provides a storage area in the DRPROG segment 
for incident solar fluxes that are read in from the NDI file. 

QDS - An array of incident solar fluxes 


LABELED COMMON /DRSHDC/ 

This common block is used to store the sequence numbers of pos- 
sible shadowing surfaces in the preshadowing calculations of 
segment DRCAL. 

ISHAD - An array of sequence numbers of possible shadowing sur- 
faces 
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LABELED COMMON /DRSUN/ 


This common block contains information used in calculating the 
incident solar flux on a node in the DEPROG segment, 

PACT - Unshadowed solar flux on a node 

Y “ Square of the distance from the node to the sun 

Z - Ratio of the solar constant to the node-to-sun distance 


LABELED COMMON /DRSURl/ 

This common block serves the same purpose in the DEPROG segment 
that DISURl serves in the DIPROG segment (see the DISURl descrip- 
tion) . 


LABELED COMMON /DRTRAN/ 

This common block contains the matrix of direction cosines neces- 
sary to transform vectors defined in the planet-oriented vehicle 
coordinate system to vectors in the user-defined vehicle coordi- 
nate system, 

PLDC - Transformation matrix 


LABELED COMMON /DRTRSH/ 

This common block contains miscellaneous variables used to cal- 
culate Incident fluxes in the DRPROG segment. 


IN - Sequence number of the current node 
NCHECK - Not used 

KELT - Total number of elements used on the node 
NSHAD - Number of possible shadowing surfaces 
NSHADR - Not used 

SFAVT - Temporary array of elemental surface area vectors 
SFPVT - Temporary array of elemental surface position vectors 
SHADS - Shadow factor 
SUNPVT " Temporary sun position vector 
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LABELED COMMON /DRVCTR/ 

This common block contains vector information for use in calcu- 
lating incident fluxes in the DEPROG segment. 

WEPT - Not used 

WEST - Optimum number of elements used on the node 
SDWPV - Sun position vector 

SFAV “ Array of elemental surface area vectors 
SPPV - Array of elemental surface position vectors 


LABELED COMMON /DSTORE/ 

This common block contains information used in storing and re- 
trieving data from units assigned to TRASYS. 

IDSTR (I, 1) - Unit Identifier (NUWIT) 

IDSTR (I, 2) - Number of files written to NUNIT 

IDSTR (Ij 3) - Yes/'.io flag indicating, whether or not the unit has 
been repositioned 


LABELED COMMON /DSTR/ 


This common block contains surface description data for shadowing 
surfaces . 


DSTR (1, 
DSTR (2, 
DSTR (3, 
DSTR (4, 
DSTR (5, 


I) - 'ALPHA 
I) - BMIN 
I) - BMAX 
I) - GMIN 
I) - GMAX 


LABELED COMMON /EMISS/ 

This common block contains IR emissivities . 

EMISS “ Array of IR emissivities for the active nodes 


LABELED COMMON /FA/ 

This common block contains one row of a symmetric matrix. 
FA - Area-form factor product, or area-script F product 
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LABELED COMMON /FFDATl/ 


This common block contains miscellaneous variables used in the 
FFCAL segment. 

FANS “ Unshadowed solar form factor 
FENS - Unshadowed IR form factor 

lABTME " Abort flag set when maximum number of elements per node 
is exceeded 

ICALTP - Flag indicating whether the form factors came from cards, 
tape, were calculated, or Xi?ere equivalenced 
lEOFFR - End-of-file f ag on unit NFFR 
IN - Sequence number of the current "looker" node 

JN “ Sequence number of the current "lookee" node 

KTAE - Not used 

NEX - Flag to indicate when nodes have been divided into sub 
nodes 

TIMEF - Cp time at beginning of a form factor calculation 

TIMER “ CP time at beginning of a row of form factor calculations 

TIMET - Time at the beginning of the form-factor calculations 
TPER - Not used 


LABELED COMMON /FFEQ/ 

This common block contains the master index for the random-access 
file (NSCR3) used in equivalenced form factors. 

INDXF - Master index 


LABELED COMMON /FFSHDC/ 

This common block serves the same purpose in the FFPROG segment 
as DRSHDC serves in the DRPROG segment (see the DI^HDC descrlp- 
C ion) . 


L.VBELED COMMON /FFSUMC/ 

This common block contains form-factor sums. 
SUM - An array of foriii- factor sums 
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LABELED COMMON /FFVALI/ 


This common block contains a row of IR area-form factor products. 
FFVALI - IR area-form factor products 


LABELED COMMON /FFVALS/ 

This common block contains a row of solar area-form factor prod- 
ucts. 

FFVALS - Solar area-form factor products 


LABELED COMMON /GBIR/ 

This common block contains a row of IR gray-body factors for use 
in the AQPROG segment. 

GBIR - IR gray-body factors 


LABELED COMMON /GBSO/ 

This common block contains a row of solar gray-body factors for 
use in the AQPROG segment. 

GBSO - Solar gray-body factors 


LABELED COMMON /TFS/ 

This common block contains an array of sequence numbers for 
shadowing surfaces. 

IFS - Shadowing-surface sequence numbers 


LABELED COMMON /IKS/ 

This common block contains an array of surface-type identifiers 
for shadowing surfaces. 

IKS - Type ot shadowing surface 
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LABELED COMMON /INDX/ 

This common block contains the master index for the random-access 
file NRAN. 

INDX - Master index 


LABELED COMMON /INDXN/ 

This common block contains an array of random-access .record num- 
bers for active nodes, 

INDXN - Array of random-access record numbers for active nodes 


LABELED COMMON /INDXS/ 

This common block contains an array of random-access record num- 
bers for active surfaces. 

INDXS - Array of random-access record numbers for active surfaces 


LABELED COMMON /ISHAD/ 

This common block serves the same purpose in the DIPROG segment 
as DRSHDC serves in the DRPROG segment (see the DRSHDC descrip- 
tion) . 


LABELED COMMON /ISPEC/ 

This common block contains an array of active specular-surface 
sequence numbers. 

ISPEC - Array of sequence numbers of active specular surfaces 


LABELED COMtION /ISPN/ 

This common block contains an array of primary MESS nodes and/or 
special nodes for use in the RCPROG segment. 

ISPN - Array of primary MESS nodes and/or special nodes 
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EEPRODUCIBILiTY OF THE 
PEIGDsfAL PAGE pr'iOR 

LABELED COMMON /ISTPDR/ 

This common block contains a directory of user-assigned step 
numbers. 

ISTPDR - Array of user-assigned step numbers 


LABELED COMMON /JSURF/ 

This common block contains an array of all surface numbers defined 
by BUILDC/ADD. 

JSURF - Array of surface numbers 


LABELED COMMON /LNGSEG/ 


LAQSEG 

- 

Length 

of 

AQ segment 




LCMCOM 

- 

Length 

of 

dynamic core 

for 

CM 

segment 

LCMSEG 

- 

Length 

of 

CM segment 




LDISEG 

- 

Length 

of 

DI segment 




LDRSEG 

- 

Length of 

DR :>egment 




LFFSEG 

- 

Length 

of 

FF segment 




LGBCOM 

- 

Length 

of 

dynamic core 

for 

GB 

segment 

LGBSEG 

- 

Length 

of 

GB segment 




LMFSEG 

- 

Not used 





LNPSEG 

- 

Length 

of 

NP segment 




LODSEG 

- 

Length 

of 

OD segment 




LOPSEG 

- 

Length 

of 

OP segment 




LPLCOM 

- 

Length 

of 

dynamic core 

for 

PL 

segment 

LPLSEG 

- 

Length 

of 

PL segment 




LQOCOM 

- 

Length 

of 

dynamic core 

for 

QO 

segment 

LQOSEG 

- 

Length 

of 

QO segment 




LRBSEG 

- 

Length 

of 

RB segment 




LRCCOM 

- 

Length 

of 

dynamic core 

for 

RC 

segment 

LRCSEG 

- 

Length 

of 

RC segment 




LRDSEG 

- 

Length 

of 

RD segment 




LSFSEG 

- 

Length 

of 

SF segment 




MAXFL 

- 

Maximum field length av 

ailable 


NCURFL 

- 

Current 

: field length 
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LABELED COMMON /MNP/ 


This common block contains a list of selected nodes to be plotted 
in the NPPROG segment, 

MNP - Array of selected node numbers to be plotted 


LABELED COMMON /MSND/ 

This common block contains an array of secondary MESS node num- 
bers '!'or use in the RCPROG segment. 

MSND - Array of secondary MESS node numbers 


LABELED COMMON /MSP/ 

This common block contains a list of selected surfaces to be 
plotted in the OPPROG segment. 

MSP - Array of selected surface numbers to be plotted 
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LABELED CO^^MON /NCONST/ 


This conunon block contains miscellaneous variables for use in Che 
NPPROG segment. 


ITITLE - Title to be written on each plot frame 
KC - Flag indicating which view is to be plotted 
KERR - Scaling error flag 
LINPL -* Litia plot conttol flAg 

NOMI Values to indicate the order of rotation of Che plotter 
NPHI - coordinate system through the angles OMI, PHI, and PSI, 
NPSI respectively. 

OMI - Angle of rotation of the plotter coordinate system (PCS) 
about the X-axis 

PHI - Angle of rotation of the PCS about the Z-axis 

PSI - Angle of rotation of the PCS about the Y-axis 

SCL - Scale factor for plotting 

VROT - Transformation matrix of direction cosines resulting 

from the rotations PHI, PSI, and OMI. Transforms vectors 
in the PCS to the user-defined view 


LABELED COMMON /NDS/ 

This common block is used in the RCPROG segment for temporarily 
storing node numbers. 

NODET - Temporary node array 


LABELED COMMON /NODE/ 

11) is common block contains an array of the active node numbers, 
NODE - Array of active node numbers 


LABELED COMMON /NOROLO/ 

Tills common block contains information pertinent to partitioning 
a lower triangular matrix for use in the GBPROG segment. 

NLO (i, - Number of rows in block I 

NLO (2, 1) - Number of elements in block I 
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LABELED COMMON /NOROUP/ 

This cotranon block contains information pertinent to partitioning 
an upper triangular matrix for use in the GBPROG segment. 

NUP (1, I) - Number of rows in block I 
NUP (2, I) - Number of elements in block I 


LABELED COMMON /NPMASS/ 

This common block is used to store surface data read in from the 
random-access file in the NPPROG segment. 

CC - Transformation matrix relating the SCS to the CCS 
DATA - Array of the five surface description parameters, ALPHA, 
BMIN, BMAX, GMIN, and GMAX 
ILK - Identifier for type of surface 
ElMASS - Dummy array of 11 xrords 
RX -j 

RY "I Components of the position vector locating the origin of 
RZ “) the SCS in the CCS 


LABELED COMMON /NSPEC/ 

This common block contains the number of active specular surfaces 
NSPEC - Number of active specular surfaces 


LABELED COMMON /OCONST/ 

This common block serves the same purpose in the OPPROG segment 
as NCONST serves in the NPPROG segment (see the NCONST descrip- 
tion), except for SCL and the addition of txtfo variables. 

RPLN - Planet radius in plot frame dimensions 

SCL - Maximum spacecraft dimension, measured from the CCS origin 
in plot frame dimensions 
SCLR - Orbit radius in plot frame dimensions 
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LABELED COMMON /ODTEMP/ 


This common block provides a scratch array available to the pro- 
grammer within any given segment. 

ODTEMP - Scratch array dimensioned to the maximum of 100 and the 
total number of nodes input in the Surface Data block 


LABELED COMMON /OPMASS/ 

TJiis common block serves the same purpose in the OPPROG segment 
as NPMASS serves in the NPPROG segment (see the NPMASS descrip- 
tion) 


LABELED COMMON /OPTRAN/ 

This common block contains variables used in the OPPROG segment. 

S - Transformation matrix to locate the solar position vector 
in the plotter coordinate system 
CTR - Transformation matrix to transform vectors in the CCS to 
vectors in the orbit coordinate system (OCS) 

XV -j 

YV -/ Components of the vehicle position vector in the OCS 
ZV “) 
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LABELED COMMON /ORBIT/ 


This labeled common contains all variables associated with orbit 
definition and other orbit-associated parameters. 

ALAN - Longitude of the ascending node 

APER - Argument of perifocus 

ASUN - Not used 

ATMT - Atmospheric transmissivity 

BETA “ Angle from the Zp-axis of the orbit coordinate system to 
the sun vector (vertex at planet center) 

BETAS - Same as BETA, except star vector 

CIGMA - Angle, in the orbital plane, from periapsis to the pro- 
jection of the solar vector in the direction of the sp 
spacecraft's motion 

CIGMAS “ Same as CIGMA, except star vector 
CLOCK - Clock angle to the spin vector 
COKE - Cone angle to the spin vector 

DAWN - Sunrise time 

DUSK - Sunset time 

DWP - Subsolar planet emissive power, less the darkside 

emissive power 

ECC - Orbit eccentricity 

GRAV - Planet gravitational constant 

HA - Altitude of apoapsis 

HP - Altitude of periapsis 

ICALFL - Not used 

INSHAD - Flag indicating whether the point is in or out of the 
planet shadow 

lORBIT - Flag for type of orbit 
* 1 Planet 
= 2 Sun 
= 3 Star 
= 4 Tape 

IROTX - Order of performing rotation ROTX (1, 2, or 3) 

IROTY “ Order of performing rotation ROTY (1, 2, or 3) 

IROTZ - Order of performing rotation ROTZ (1, 2, or 3) 

ISFT - Flag directing the use of shadow factor data 

ISKPSO - Not used 

NSPPF - Step number for storing the spacecraft-to-planet form 
factors 

OINC - Orbit inclination 

ORNT - Transformation matrix from the vehicle coordinate system 
to the central coordinate system 
PALB - Planet albedo (solar reflectivity) 

PERIOD - Orbital period 

PLGL - Clock angle-to-planet position vector (in the CCS) 

PLTYPE - Spacecraft-planet form factor read/store flag 
PNAME - Name of planet being orbited 


REV 1 
III-174 



PRAD 

PSD 

RATE 

ROTX 

ROTY 

ROTZ 

RSUN 

RTHET 

SHADIN 

SHAOUT 

SOL 

SOLO 

SPINT 

STRRA 

STRDEC 

SUNCO 

SUNCL 

SUNDEC 

SUNPVO 

SUNRA 

STRDEC 

STRRA 

TIMEPR 

TIMES T 

TIMSP 

TRUEAN 

WDS 

WSS 

WSUN 


Planet radius 

- Planet-sun distance (set to l.OE + 15 for planetary orbits) 

- Spin rate 

- Rotation about the X-axis from the VCS into the CCS 

- Rotation about the Y-axis from the VCS into the CCS 

- Rotation about the Z-axis from the VCS into the CCS 

- Radius of sun 

- Geocentric altitude 

- True anomaly when entering the planet shadow 

- True anomaly when leaving the planet shadoii? 

- Solar "constant" 

- Extra - atmospheric solar constant 

- Transformation matrix from the CCS at zero spin time to 
the CCS at the current spin time 

- Right ascens'iion of star 

- Declination of star 

- Cone angle to sun vector (in the CCS) 

- Clock angle to sun vector (in the CCS) 

- Declination of the sun 

- Sun vector at noon 

- Right ascension of the sun 

- Declination of star 

- Right ascension of star 

- Current problem time 

- Time of periapsis passage 

- Time at which spinning begins 

- Orbit true anomaly 

- Planet's darkside emissive power 

- Subsolar emissive pother of the planet 

- Emissive power of the sun 
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LABELED COMMON /PLANET/ 


This common block contains variables used to compute the planet/ 
spacecraft geometry relationship and the planet element break- 
do\m . 


ALB 

ALG 

APEAPL - 
BETPM - 


E 


POSP 

RADP 

RADPJ - 


Product of the planet radius and BETPM 
AREAPL/ALB 

Area of the portion of the planet visible from the 
spacecraft 

Angular measurement from the subspacecraft point on the 
planet to the planet horizon, with the center of the 
planet as the vertex 
(PRAD + H)*PRAD 
where ; 

PRAD = Planet radius 
H = Altitude of spacecraft 

Components of a vector from the CCS origin to the center 
of the planet 

Radius of a preshadowing sphere associated with a planet 
element 

Arra> containing radii of the preshadowing spheres asso- 
ciated with each planet element 
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LABELED COMMON /PLOT/ 


This common block contains arrays of infcrmation used in the 
NPPROG and OP’^ROG segments. 

lOPNNP - An array that is functionally analogous to NPNNP for use 
in the OPPROG segment 
lOPNV - Not used 

lOPTIT - An array that is functionally analogous to NPTIT for 
use in the OPPROG segment 

lOPVU - An array of viexjs to be plotted in the OPPROG segment 

NPNNP - An array of starting locations in labeled common ARRAYS 

for user-input arrays of selected nodes to be plotted 
in the NPPROG segment 

NPTIT - An array of starting locations in labeled common ARRAYS 
for user-input title arrays that are to be written on 
plot frames in the NPPROG segment 
NPVU “ An array of views to be plotted in the NPPROG segment 

OPROT - An array of user-input Euler angles for general views 

in the OPPROG segment 

OPRPLK “ An array of the desired plot sizes of the planet radius 
for use in the OPPROG segment 

OPSCL “ An array specifying the desired plot size for maximum 
spacecraft dimensions, as measured from the CCS origin 
for use in the OPPROG segment 

OPSCLR - An array specifying the desired plot sizes of the orbit 
radius for use in the OPPROG segment 
OPTIMP - An array of present times in the orbit used in conjunc- 
tion with OPTIMS to calculate true anomalies in the 
OPPROG segment 

OPTIMS - An array of perigee passage times for use in the OPPROG 
segment 

OPTRUE - An array of true anomalies for use In the OPPROG segment 
ZNPROT - An array of user-input Euler angles for general views 
in the NPPROG segment 

ZNPSCL - An array of scale factors for use in the NPPROG segment 


LABELED COMMON /PLOTTR/ 

This labeled common contains the variables, flags, and Hollerith 
title data for controlling data plot operations, 

IPLNA - Array of node numbers for selective data plotting 

IPLSN - Array of step numbers for selective data plotting 

IPLUNT - Plot flag for type of data 

PLCRVF - Yes-no flag for curve fitting data plots 

PLLABX - X-axis plot label array 

PLLABY - Y-axis plot label array 

PLTITl - Plot label title 1 

PLTIT2 - Plot label title 2 

PLXMPF - Plot multiplying factor for X-axis 

PLYMPF - Plot multiplying factor for. Y-axis 
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labeled common /PR/ 


A 2xNS array containing the IR and solar transmissivities asso- 
ciated with each surface. 

PR (1, N) - IR Transmissivity of surface N 
PR ( 2 , N) - Solar transmissivity of surface N 


LABELED COMMON /PSH/ 

A 4xNS array containing data associated with the preshadowing 
sphere for each surface, 

PSH (1, N) - Radius of the preshadowing sphere for surface N 
PSH (2, N) 

PSH (3, N) - > Components of the position vector to the center of 
PSH (A, N) -) the preshadowing sphere for surface N 


LABELED COMMON /QAP/ 

This labeled common contains an array of absorbed planetary heat 
rates for each node. 

QAP (N) - Absorbed planetary heat rate (energy/unit time) for 
node N 


LABELED COMMON /QAR/ 

This labeled common contains n array of absorbed reflected 
(planetary albedo) heat rates for each node, 

QAR (N) - Absorbed albedo heat rate (energy/unit time) for node N 


LABELED COMMON /QAS/ 

This labeled common contains an array of absorbed solar heat rates 
for each node. 

QAS (N) - Absorbed solar heat rate (energy /unit time) for node N 
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LABELED COMMON /QAVERG/ 

This labeled common contains an array of orbital average absorbed 
total heat rates for each node. 

QAVERG (N) - Average absorbed heat rate for node N 


LABELED COMMON /QDP/ 

This labeled common contains an array of direct planetary heat 
fluxes for each node. 

QDP (K) - Direct planetary heat flux (energy/unit time/unit area) 
for node N 


LABELED COMMON /QDR/ 

This labeled common contains an array of reflected (planetary 
albedo heat fluxes for each node. 


LABELED COMMON /QDS/ 

This labeled common contains an array of solar heat fluxes for 
each node. 


labeled coition /QOCMB/ 

This labeled commoi. contains an array used for working storage of 
correspondence dat a when combining absorbed-Q data for output. 

ICOMB - a. An array of node numbers identical with the list on 

the right side of an equal sign in correspondence data 
b. On the second pass of the combining logic 5 node num- 
bers duplicated in the correspondence data are set 
negative 
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I^BELED COMMON /QOFRST/ 


This labeled common contains an array of pointers that refer to 
node numbers duplicated in the correspondence data. 

IFRST - An array of sequence numbers pointing to node numbers on 
the left side of equal signs in the correspondence data 
that are duplicated by node numbers on the right of any 
equal sign 


LABELED COMMON /QONODT/ 

This labeled common contains a temporary-node-number array. 

NODET - A temporary-node-number array used in node combining oper- 
ations. Combine operations begin in a manner identical 
to those in the NODE (uncorabined) array and end as a com- 
bined node array 


LABELED COMMON /RBDATl/ 

This common block serves the same purpose in the RBPROG segment 
that FFDAT?. serves in the FFPROG segment (see the FFDATl descrip- 
tion) . 


LABELED COMMON /RBFFVI/ 

This common block contains a row of IR area-form factor products 
that were read in from the NSCRl file or the NFFR file. 

FFVALI - IR area-form factor products 


LABELED COMMON /RBFFVS/ 

This common block contains a row of solar area-form factor pro- 
ducts that were read in from the NSCRl file or the NFFR file. 

FFVALS - Solar area-form factor products 
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LABELED COMMON /RBSHDC/ 


This conmon block serves the same purpose in the RBPROG segment 
as DRSHDC serves in the DEPROG segment (see the DRSHDC descrip- 
tion) . 


LABELED COMMONS /RBSURl/ and /RBSUR2/ 

These common blocks serve the same purpose in the RBPROG segment 
as DISURl serves in the DIPROG segment (see the DISURl descrip- 
tion) . 


LABELED COMMON /RBVALI/ 

This common block contains a row of IR area-image factor products, 
RBVALI - IR area- image factor products 


LABELED COMMON /RBVALS/ 

This common block contains a row of solar area-image factor 
products. 

RBVALS - Solar area-image factor products 


LABELED COMMON /RKCMB/ 

This common block serves the same purpose in the RKPROG segment 
as QOCMB serves in the QOPROG segment (see the QOCMB descrip- 
tion) , 


LABELED COMMON /RKFRST/ 

This common block serves the same purpose in the RKPROG segment 
as QOERST serves in the QOPROG segment (see the QDFRST descrip- 
tion) . 
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LABELED COMMON /RKNODT/ 

This common block serves the same purpose in the EKPROG segment as 
QONODT serves in the QOPROG segment (see the QONODT description). 


LABELED COMMONS /RMASSl/ and /RMASS2/ 

These common blocks serve the same purpose in the FFPROG segment 
as DISURl serves in the DiPROG segment (see the DISURI descrip- 
tion) . 


L/iBELED COMMON /RMASSM/ 

This common block serves the sdme purpose in the RBPROG segment 
as DISURl serves in the DIPROG segment (see the DISURl descrip- 
tion) . 
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LABELED COMMON /RNDMRP/ 


This common block contains the random-access record numbers and 
record lengths that are passed from the preprocessor to the pro- 
cessor to enable the processor to access preprocessor-generated 
data from the random access file. These variables are used in 
the RDPROG segment to set up processor common blocks. 


NLAD - The record length of the user array directory 

array 

NLRIO - The maximum number of records that can exist on 

the random I/O file 

NRAND - The random I/O record number of the user array 

name directory array 

NRAPD - The random I/O record number of the user array 

position directory array 

NRAV, NLAV — The random I/O record number and the record 

length of the user array value array. Used to 
set up common ARRAYS 

NRBCSD, NLECSD - The random I/O record number and the record 
length of the BCS directory 

NRBCSN, NLBCSN - The random I/O record number and the record 

length of the BCS name array. Used to set up 
common BCSN 

NRBCSRj NLBCSR - The random I/O record number and the record 

length of the BCS index. Used to set up common 
BLKDIR 


NRCQD, NLCQD 
NRCQV, NLCQV 

NR IN, NLIN 
NRIS, NLIS 
NRSD, NLSD 

NRT, NLT 

NRTD, NLTD 
NRUQD, NLUQD 
NRUQV, NLUQV 


- The random l/o record number and the record 
length of the control constants directory array 

" The random I/O record number and the record 
length of the control constants value array. 
Used to set up common CCONST 

- The random I/O record number and the record 
length of the node index 

- The random I/O record number and the record 
length of the surface index 

- The random I/O record number and the record 
length of the step directory. Used to set up 
common ISTPDR 

- The random I/O record number and the record 
length for the title array. Used to set up 
common TITLE 

- The random I/O record number and the record 
length of the combined directory 

- The random I/O record number and the record 
length of the user constant directory array 

- The random I/O record number and the record 
length of the user constant value array. Used 
to set up common UCONST 
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LABELED COMMON /SF/ 


This labeled common contains a storage array for script— F (gray- 
body factor) data, 

SF — An array used to store one row of script— F values during node- 
combining operations in the RKPROG segment 


LABELED .COMMON /SPS/ 

This labeled common contains a storage array for script-F data. 

SFS - An array used to store one row of script-F values during 
node- combining operations in the RCPROG segment 


LABELED COMMON /SFSHDC/ 

This labeled common used in the SFPROG segment in same manner as 
labeled common DRSHDC in DRPROG. 


LABELED COMMON /SFQDP/ 
Not used 


LABELED COMMON /SFQDR/ 
Not used 


LABELED COMMON /SFQDS/ 
Not used 


LABELED COMMON /SFSURl/ 

This labeled common contains the same variable names and is used 
identically with labeled common DISURl. 
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LABELED COMMON /SFVECC/ 


This labeled common contains variables used to define surface and 
planet elements used in shadow-factor tape calculations- 

NEPT - Number of elements on the planet 
NEST - Number of elements on the node 

PLAV - Planet- element area vector array (3 components in the CCS 
for each element; 

PLPV - Planet-element position vector array 
SFAV - Surface— element area vector array 
SFPV - Surface-element position vector array 
SDWPV - Sun position vector 


LABELED COMMON /SPACE/ 

This labeled common contains an array of script-F (gray-body) 
factors from each node to space. 


LABELED COMMON /SPACNO/ 

This labeled common contains radiation conductor values from each 
corrbined node to space, and is used in the RKPROG segment. 


LABELED COMMON /SPCNO/ 

This labeled common, used in RCPROG, is exactly analogous to 
SPACNO. 


LABELED COMMON /SREFLI/ 

This common block contains an array of surface IR specular re- 
flectivities used in the DRPROG and RBPROG segments, 

SREFLI - Array of IR specular reflectivities for active specular 
surfaces 



labeled common /s reels/ 


This common Block contains an array of surface solar specular re- 
flectivities used in the DEPROG and RBPROG segments. 

S REELS - Array of solar specular reflectivities for active specular 
surfaces 


LABELED COmON /SRIR/ 

This common block contains an array of nodal IR specular reflec- 
tivities used in the DKPRGG and the RBPROG segments. 

SRIR - Array of IR specular reflectivities for active nodes 


LABELED COMMON /SRSO/ 

This common block contains an array of nodal solar specular re- 
flectivities used in the DRPROG and the RBPROG segments, 

SRSO - Array of solar specular reflectivities for active nodes 


LABELED COMMON /SUN/ 

This labeled common is identical to labeled common DRSUN. 


LABELED COMMON /TAPE/ 

This labeled common contains all variable names for the processor 
library files (see Section III-C) . 
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LABELED COMMON /TITLE/ 

This labeled common contains title information. 


TITLE 

- 

Problem title input by the user 

NTITLE 


Title identifying the segment in the segment 
header printouts 

DTE 

- 

Sun data 

IHSTEP 

- 

Step number in hollerith format 

IPAGE 

- 

Page number counter 

LINE 

- 

Line number counter 

MLINE 

- 

Maximum lines allowed per page 

MODELN 

“ 

Configuration name 

NJOB 

- 

Job identification number 

TME 


Time of day 
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LABELED COMMON /TRANS/ 


This labeled common contains a transformation matrix to transform 
planet-oriented vectors to the user-defined VCS . 

PLDC - Transformation matrix used to transform vectors from a 
planet-centered coordinate system to the VCS 


LABELED COMMON /TRASH/ 

This labeled common contains miscellaneous vector arrays and 
variables used in direct flux calculations. 

IN - Current-node sequence number 

NCHECK - Check flag to eliminate unnecessary calls to the pre- 
shadowing routine 

NELT - Temporary storage address for the number of elements on 
a node during element optimization operations 
NSHAD - Number of possible shadowing surfaces associated with a 
node 

NSHADR - Counter used in conjunction with NCHiICk! to eliminate 
unnecessary calls to pre-shadowing routine 
PLAVT - Temporary storage address for planet-element area vectors 

PLPVT - Temporary storage address for planet-element position vectors 

SFAVT - Temporary storage address for node-element area vectors 

SFPVT - Temporary storage address for node-element position vectors 

SUNPVT “ Temporary storage address for the sun position vector 


LABELED COMMON /TRIR/ 

This common block contains an array of nodal IR transmissivities. 
TRIR - Array of IR transmissivities for active nodes 


LABELED COMMON /TRSO/ 

This common block contains an array of nodal solar specular 
transmissivities , 

TRSO - Array of solar transmissivities for active nodes 
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LABELED COMMON /TSTR/ 


Tliis common block contains an array of transformation matrices for 
shadowing surfaces, 

TSTR - Array of transformation matrices relating the SCS coor- 
dinates of shadowing surfaces to the CCS 


LABELED COMMON /UCONST/ 

This array contains user-defined constants that were input in the 
Quantities Data block. 


LABELED COMON /VARBL/ 

This common block contains variables used to calculate albedo and 
planetary fluxes in the DIPROG segment. 

ALBF - Product of the solar constant and form factor from an in- 
solated planetary element to the sun 
PLNF - Emissive power of a planetary element 


LABET.ED COMMON /VECTOR/ 

This labeled common contains the same variable names, used in the 
same way as those in labeled common SFVECC. 


LABELED COMMON /XSPACE/ 

This labeled common contains a scratch array used to store one 
row of the inverted gray-body factor matrix. 


BLANK COMMON 

Blank common in the TRASYS processor is used as a scratch storage 
area in central memory. 


E, EMTRY POINTS PROCESSOR LIBRARY 


EHTRY POmr SUBROUTINE 
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BUILDC 
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INITF 
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PFINIT 

INITF 
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INITF 
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